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ABSTRACT 


Three  new  methods  for  the  detection  and  interception  of  frequency-hopped 
waveforms  are  presented.  The  first  method  extends  the  optimal,  fixed-block 
detection  method  based  on  the  likelihood  ratio  to  a  sequential  one  based  on  the 
Sequential  Probability  Ratio  Test  (SPRT).  The  second  method  is  structured 
around  a  compressive  receiver  and  is  highly  efficient  yet  easily  implemented. 
The  third  method  is  based  on  the  new  concept  of  Amplitude  Distribution  Func¬ 
tion  (ADF)  and  results  in  a  detector  that  is  an  extension  of  the  radiometer. 

The  first  method  presents  a  detector  structured  to  make  a  decision  sequen¬ 
tially,  that  is,  as  each  data  element  is  collected.  Initially,  a  purely  sequential 
test  is  derived  and  shown  to  require  fewer  data  for  a  decision.  A  truncated 
sequential  method  is  also  derived  and  shown  to  reduce  the  data  needed  for  a 
decision  while  operating  under  poor  signal-to-noise  ratios  (SNRs).  A  detailed 
performance  analysis  is  presented  along  with  numerical  and  Monte  Carlo  anal¬ 
yses  of  the  detectors. 

The  second  method  assumes  stationary,  colored  Gaussian  interference  and 
presents  a  detailed  model  of  the  compressive  receiver.  A  locally  optimal  de¬ 
tector  is  developed  via  the  likelihood  ratio  theory  and  yields  a  reference  to 
which  previous  ad  hoc  schemes  are  compared.  A  simplified,  suboptimal  scheme 
is  developed  that  trades  off  duty  cycle  for  performance,  and  a  technique  for 
estimating  hop  frequency  is  developed.  The  performance  of  the  optimal  and 
suboptimal  detectors  is  quantified.  For  the  suboptimal  scheme,  the  trade-off 
with  duty  cycle  is  studied.  The  reliability  of  the  hop  frequency  estimator  is 
bounded  and  traded  off  against  duty  cycle. 

In  the  third  method,  a  precise  definition  of  the  ADF  is  given,  from  which 
follows  a  convolutional  relationship  between  the  ADFs  of  signal  and  additive 
noise.  A  technique  is  given  for  deconvolving  the  ADF,  with  which  signal  and 
noise  components  can  be  separated.  A  detection  statistic  based  directly  on 
this  deconvolution  technique  is  defined  and  statistically  characterized,  yielding 
a  framework  on  which  to  synthesize  a  detector.  The  detector’s  performance  is 
analyzed  and  compared  with  the  radiometer. 
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CHAPTER  1 

INTRODUCTION  AND 
BACKGROUND 

1.1  MOTIVATION 

The  goal  of  the  intercept  receiver  is  to  detect,  identify,  and  geolocate  hostile 
electromagnetic  (EM)  sources  and  use  this  information  to  counter  with  Elec¬ 
tronic  Counter  Measures  (ECM)  and  Electronic  Counter  Counter  Measures 
(ECCM).  For  military  communications,  the  receiver  is  expected  to  intercept 
any  one  of  a  set  of  target  communications.  The  interception  task  is  hindered 
by  a  dense  EM  environment  that  consists  of  other  communication  signals,  both 
friendly  and  hostile,  possibly  jamming  or  masking  signals  in  addition  to  ever¬ 
present  noise  interference.  With  the  advent  of  frequency-hopped  and  other 
spread-spectrum  communications,  the  search  bandwidth  that  assures  a  rea¬ 
sonable  probability  of  intercept  is  greatly  increased,  thereby  aggravating  the 
problem  of  interference,  because  greater  numbers  and  types  of  interferences 
obscure  th e  *•"»-«>*  signal  Th"  'P^reased  complexity  of  the  interception  prob¬ 
lem  motivates  the  search  for  new  methods  of  detection  and  interception  of 
frequency-hopped  waveforms. 

1.2  TARGET  SIGNALS 


Military  and  other  secure  communications  use  spread-spectrum  signaling 
involving  some  variety  of  modulation  whose  purpose  is  to  add  ambiguity  or 
“randomness”  to  the  waveform  as  a  measure  against  unintended  detection  or 
interception.  The  usual  procedure  for  randomizing  the  waveform  is  pseudo¬ 
random  variation  of  transmission  times  (time  hopping  or  TH),  phases  (direct 
sequence  or  DS),  or  frequencies  (frequency  hopping  or  FH).  This  work  concen¬ 
trates  solely  on  the  interception  of  FH  waveforms  that  have  form 


N* 

s(0  =  53^(0 


(1.1) 


where 

n(t)  equals  \fl. S'  sin (w*,<  +  6i)  for  i7\  <<<(»+  1)7*; 


is  a  family  of  known  frequencies  within  the  spread-spectrum  band¬ 
width; 
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{k,}  are  integer- valued,  independent,  uniformly  distributed,  random  vari¬ 
ables  ranging  inclusively  between  1  and  A’; 

{$,}  are  continuous,  independent,  uniformly  distributed,  random  variables 
ranging  between  0  and  2tr  that  represent  carrier  phase; 

S'  is  a  real  constant  denoting  the  average  signal  energy; 

T),  is  a  real  constant  denoting  the  epoch,  or  time  duration,  of  each  hop; 

Nh  is  a  positive  integer  denoting  the  number  of  hops  during  message 

transmission. 

This  general  model  of  frequency-hopped  waveforms  includes  a  large  number  of 
modulations  such  as  frequency  shift  keying  (FSK)  and  minimum  shift  keying 
(\1SK).  Some  important  modulations  not  included  are  those  whose  carrier  phase 
is  correlated  from  hop  to  hop,  for  example,  continuous  phase  FSK  (CPFSK). 
Even  for  these  cases,  these  results  apply  but  may  not  be  optimal. 

1.3  INTERCEPT  RECEIVER  FUNCTIONS 


An  intercept  receiver  extracts,  for  further  processing,  a  small  number  of 
candidate  signals  from  the  plethora  of  signals  in  a  communication  band  of 
interest.  The  initial  processing  steps  that  discard  signals  are  called  Pruning 
Functions.  After  pruning,  secondary  processing,  known  as  feature  extraction, 
yields  information  aiding  in  emitter  identification  and  countering.  After  feature 
extraction,  further  processing  could  yield  the  actual  information  embedded  in 
•he  communication  signal,  but  this  is  peripheral  to  the  primary  function  of  the 
intercept  receiver,  namely  ECM  and  weapons  support,  and  consequently  will 
not  be  explored  here. 

1.3.1  Pruning 

Given  the  frenetic  activity  in  most  communication  bands,  the  intercept  re¬ 
ceiver  must,  early  in  its  processing,  choose  out  of  all  candidate  signals  within 
the  band  a  small  number  of  potential  target  signals.  Pruning  Functions  achieve 
this  by  eliminating  all  but  the  most  promising  prospects  for  processing.  Prun¬ 
ing  Functions  fall  into  four  categories:  Initial  Detection,  Direction  Finding,  Fre¬ 
quency  Estimation,  and  Time  of  Arrival  Measurement.  As  is  evident  from  their 
names,  Initial  Detection  separates  potential  target  signals  from  background 
noise,  Direction  Finding  classifies  and  possibly  eliminates  signals  by  direction 
of  origin,  Frequency  Estimation  censors  signals  based  on  a  measurement  of  car¬ 
rier  frequency,  while  Time  of  Arrival  Measurement  differentiates  between  pulsed 
signals  based  on  their  arrival  times.  It  is  useful  conceptually  to  consider  the 
Pruning  Functions  as  independent  processes  applied  separately,  but  in  a  practi¬ 
cal  system  these  functions  are  usually  highly  coupled  in  that  a  single  processing 
step  may  accomplish  two  or  more  Pruning  Functions. 

Initial  Detection  separates  the  candidate  signal  from  the  background  noise, 
usually  in  the  form  of  a  threshold  operation  applied  on  a  test  statistic  derived 
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from  the  received  waveform.  It  may  be  a  dedicated  operation,  meaning  that 
the  only  information  gleaned  is  the  presence  or  nonpresence  of  a  communication 
waveform,  or  may  yield  other  information  such  as  the  time  of  arrival  (time  do¬ 
main)  or  the  dominant  frequency  (frequency  domain),  or,  for  feature  detectors, 
the  hop  rate. 

The  signal  direction  can  be  found  by  using  a  scanning  narrowbeam  antenna 
whose  scan  position  at  the  time  of  detection  determines  direction.  Difference 
in  signal  phase  from  two  separate  antennas  can  yield  the  angle  from  which  a 
candidate  signal  is  emitted;  however,  this  interferometry  technique  suffers  from 
a  vulnerability  to  coherent  interference.  In  a  similar  fashion,  the  amplitude 
difference  between  the  same  signal  received  from  different  antennas  or  antenna- 
patterns  can  determine  emission  angle  but  the  technique  is  even  more  vulnerable 
to  interferences,  coherent  or  noncoherent.  The  difference  between  the  times  of 
arrival  of  different  receivers  is  an  alternative  way  to  determine  angle.  This 
method  seems  more  directed  to  pulse  signals,  such  as  radar,  in  which  arrival 
time  is  a  relatively  simple  quantity  to  measure,  but  it  could  be  applied  to  more 
complicated  communicat ion  signals,  by  cross  correlation  for  instance. 

By  estimating  the  dominant  frequency  of  a  candidate  signal,  narrowband 
interferences  such  as  other  non-spread-sp-'ctrum  communication  signals  can  be 
identified  and  rejected.  Additionally,  the  current  hop  frequency  of  the  target 
signal  can  be  determined  and  subsequently  used  to  narrowband  jam  the  current 
hop  band.  The  estimated  hop  frequency  also  can  be  used  for  identification 
processing  or  information  extraction. 

Time  of  Arrival  Measurement  can  not  only  determine  emission  angle,  as 
previously  mentioned,  but  also  can  prune.  This  use  is  primarily  useful  for 
pulse  radar  signals,  in  which  times  of  arrival  can  determine  pulse  repetition 
rates  and  hence  associate  the  intercepted  pulse  with  the  emitting  radar.  For 
frequency-hopped  communications,  the  arrival  times  of  individual  hop  intervals 
determine  the  emitters'  hop  rate  and  hence  discriminate  between  target  signal 
and  interferers. 


1.3.2  Feature  Extraction 


Feature  Extraction  is  the  measurement  of  characterizing  features  of  the  com¬ 
munication  waveforms.  Features  such  as  hop  rate,  hop  frequency,  modulation 
type,  and  bit  periods  serve  as  examples.  Feature  Extraction  overlaps  the  Prun¬ 
ing  Functions  in  that  operations  such  as  center  frequency  estimation  and  time- 
of-arriva!  measurements  yield  useful  features.  The  loose  distinction  is  one  of 
precision  and  purpose.  Measurements  made  while  pruning  are  coarse  and  serve 
only  to  decimate  what  otherwise  would  be  an  unmanageable  number  of  can¬ 
didate  signals,  while  extracted  features  are  of  sufficient  accuracy  to  serve  the 
intercept  purposes  of  jamming  and  identification. 


1.4  EXISTING  METHODS 


All  aspects  and  functions  of  the  intercept  receiver  were  described.  However, 
this  work  concentrates  solely  on  initial  detection  and  feature  extraction  and  on 
feature  extraction;  only  hop  frequency  estimation  will  be  explored. 
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Figure  1.1.  Wideband  Energy  Detector 


1.4.1  Wideband  Energy  Detector 

The  Wideband  Energy  Detector1,2  is  the  simplest  to  implement  of  all  exist¬ 
ing  detection  schemes.  Also  called  a  radiometer,  it  is  a  device  for  measuring  the 
energy  of  a  signal  over  a  prescribed  time  and  bandwidth.  A  typical  radiome¬ 
ter  (Figure  1.1)  consists  of  a  bandpass  filter  followed  by  a  square  law  device 
arid  an  integrator.  The  bandpass  filter,  of  bandwidth  W,  selects  the  frequency 
band  ever  which  the  energy  measurement  is  made.  The  square  law  device  and 
integrator  calculate  the  energy  of  this  bandpassed  signal  over  the  interval  of 
length  T.  The  radiometer  is  used  to  detect  spread-spectrum  signals  by  setting 
the  center  frequency  and  bandwidth  of  the  bandpass  filter  so  that  the  filter  will 
pass  the  expected  spread-spectrum  signal.  A  decision  is  made  by  comparing 
the  output  of  the  radiometer  to  a  threshold. 

As  can  be  guessed,  the  wideband  r  diometric  detector  is  most  efficient  when 
the  bandwidth  exactly  matches  the  spread-spectrum  bandwidth  A'/T*  and  the 
integration  period  matches  the  transmission  time  N/,  1\  of  the  spread-spectrum 
signal.  Under  these  conditions  and  for  large  time  bandwidth  (TW  >  1000) 
products,  the  p-  rformance  of  the  wideband  radiometer  is  described  below  for 
the  case  of  white-noise  interference  with  single-sided  spectral  density  No  as 


and 

d  =  Q~l{PF)-Q-l{PD)  (13) 

where  Q~l(  )  is  the  inverse  of  the  complementary  Gaussian  probability  dis¬ 
tribution  and  S'Th/No  is  the  required  signal-to-noise  ratio  for  detection  with 
probability  Pd  and  false-alarm  probability  Pf . 

In  addition  to  being  the  easiest  detector  to  implement,  the  wideband  en¬ 
ergy  detector  assumes  the  least  known  about  the  spread-spectrum  signal.  For 
optimal  detection  of  a  given  spread-spectrum  waveform,  only  the  bandwidth 
and  message  duration  need  to  be  known.  However,  an  approximate  knowledge 

*M.K.  Simon,  J.K.  Omur*,  R.A.  Scholtz,  and  B.K.  Levitt,  Spread  Spectrum  Communica- 
lions,  vol.  Ill,  Computer  Sciences  Press,  Rockville,  Md.,  1985. 

3D.  Torrieri,  Principles  of  Secure  Communication  Systems.  Artech  House,  Dedham,  Mass., 
1985. 
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of  these  parameters  degrades  performance  only  slightly.  Because  the  perfor¬ 
mance  of  the  wideband  energy  detector  is  invariant  to  the  details  of  the  spread- 
spectrum  waveform,  is  equally  effective  in  the  detection  of  either  FH,  TH,  or 
DS  waveforms.  It  is  also  useful  as  a  lower  bound  on  the  performance  of  other 
detectors  designed  around  the  particulars  of  a  given  spread-spectrum  waveform. 


1.4.2  Optimal  Channelized  Detector 


The  Optimal  Channelized  Detector2,3  uses  a  more  precise  knowledge  of  a 
spread-spectrum  waveform  to  achieve  performance  gains  over  the  Wideband 
Energy  Detector  and  other  detector  configurations.  In  the  context  of  FH  wave¬ 
forms,  the  message  duration  and  the  period  and  phase  of  the  pulse  epochs  are 
assumed  to  be  known.  Also  assumed  to  be  known  are  the  exact  frequencies 
of  the  pulsed  sinusoids  that  constitute  the  individual  “hops”  of  the  spread- 
spectrum  waveform.  The  signal  amplitude  relative  to  the  background  noise 
is  also  assumed  to  be  known.  Not  known  are  the  relative  phase  between  the 
individual  hops  and,  of  course,  the  pseudo-random  code  that  produces  the  hops. 

With  these  assumptions  and  that  of  white-noise  interference,  detection  the¬ 
ory  yields  the  detector  shown  in  Figure  1.2.  This  detector  consists  of  individual 
filters  matched  in  time  and  frequency  to  each  of  the  possible  pulsed  sinusoids 
component  to  the  FH  waveform.  The  envelopes  of  the  matched  filter  outputs 
are  “emphasized”  by  normalizing  by  expected  noise  energy  and  applying  the 
zero-order  modified  Bessel  function  of  the  first  kind  I0.  The  emphasized  filter 
outputs  are  summed  to  yield  a  likelihood  function  over  a  single  epoch.  These 
individual  likelihood  functions  for  each  epoch  of  the  message  are  multiplied  to 
yield  the  overall  likelihood  function,  from  which  a  decision  can  be  made  via  a 
thresi,  rid  comparison. 

The  generalized  performance  expression  for  the  optimal  multichannel  de¬ 
tector  cannot  be  obtained  due  to  an  inability  to  specify  the  output  probability 
distribution  functions.  When  the  number  of  hops  Nh  is  large  (e.g.,  7Vh  >  100),  it 
is  possible  to  closely  approximate  the  true  answer  by  using  Gaussian  statistics. 
This  analysis  gives  S'Th/No  needed  for  a  given  P?a  and  Pd  as 


S'Th 

No 


i-  K  +  Ke 


(1.4) 


where  d  is  given  before. 

Unfortunately,  however,  the  Optimal  Channelized  Detector  is  only  of  aca¬ 
demic  interest  because  of  its  implementation  complexity  and  its  sensitivity  to 
the  FII  waveform  parameters.  It  is  useful  primarily  for  establishing  an  upper 
bound  to  the  performance  of  other  more  implementable  and  robust  detectors. 

1.4.3  Suboptimal  Channelized  Detectors 

Because  of  the  implementation  complexity  of  the  Optimal  Channelized  De¬ 
tector,  Suboptimal  Channelized  Detectors  are  considered.  These  are  several 

3D.G.  Woodring,  “Performance  of  Optimum  and  Suboptimum  Detectors  for  Spread  Spec- 
trum  Waveform**” ,  Naval  Research  Laboratory,  Washington,  D.C.,  Technical  Report  No. 
8432,  December  1980. 
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Figure  1.2:  Optimal  Channelized  Detector 

classes  of  detectors  that  are  simplifications  in  various  ways  of  the  Optimal 
Channelized  Detector.  These  simplifications  are  listed  below. 

The  combinations  of  a  matched  filter  followed  by  an  envelope  detector  are 
replaced  by  narrowband  radiometers  of  bandwidth  l/Th-  This  simplifica¬ 
tion  degrades  the  resulting  performance  over  that  of  the  predicted  perfor¬ 
mance  of  the  Optimal  Channelized  Detector.  However,  this  performance 
difference  would  be  less  in  practice,  because  the  Optimal  Channelized  De¬ 
tector  would  not  achieve  the  optimum  performance  due  to  the  inevitable 
mismatch  between  the  actual  frequencies  received  and  the  frequencies  for 
which  the  filters  are  matched.  Doppler  shift  and  transmitter  waveform 
diversity  are  the  likely  culprits  of  this  mismatch. 

The  emphasizing  function  is  linearized.  This  is  mainly  an  implementation  sim¬ 
plification.  However,  the  optimality  of  the  Optimal  Channelized  Detector 
depends  on  a  priori  knowledge  of  the  amplitude  of  the  FH  waveform,  which 
is  a  parameter  of  the  emphasis.  Thus  the  loss  of  optimality  in  practice  is 
lessened  because  of  the  uncertainty  of  the  FH  waveform  amplitude. 

Decisions  are  made  at  the  channel  level  and  are  then  combined  to  form  a 
statistic  upon  which  the  final  decision  is  based.  This  type  of  detector  is 
appealing  if  a  frequency  estimate  of  the  detected  is  also  desired. 

Instead  of  having  a  filter  for  each  FH  frequency,  the  entire  spread-spectrum 
bandwidth  is  subdivided  into  coarse  subbands.  The  subbands  most  likely 
to  contain  the  current  hop  are  selected  for  application  of  any  of  the  above 
channelised  detector  schemes. 
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(Varying  with  time) 


Decision 


Figure  1.3:  Autocorrelation  Detector 

1.4.4  Autocorrelation  Detector 

An  Autocorrelation  Detector4,5  is  composed  of  a  bank  of  autocorrelators, 
each  operating  on  a  subband  of  the  total  spread-spectrum  bandwidth  (Figure 
1.3).  The  autocorrelators  estimate  the  autocorrelation  y(r)  of  their  bandlimited 
input  over  the  time  period  T.  The  power  of  each  correlation  is  sampled  yielding 
Wk  which  are  weighted  with  a*  and  summed  to  derive  the  decision  statistic  Y . 
This  is  summarized  by  the  equation 

G-i  Hr 

Y  =  Y,akWk  <  Threshold-  (1.5) 

*=i  Ho 

There  are  three  issues  concerned  with  the  design  of  the  autocorrelation  de¬ 
tector.  The  first  is  the  coarseness  of  the  individual  subbands  relative  to  the 
total  spread-spectrum  bandwidth.  It  seems  intuitively  appealing  to  assume 
that  performance  would  improve  by  reducing  the  width  of  this  subband  up  to 
the  limits  of  the  FH  frequency  spacing;  however,  no  analytical  or  numerical 
results  confirm  this  conjecture.  The  second  design  issue  is  the  time  interval 
over  which  we  estimate  the  autocorrelation  function.  The  third  issue  involves 
the  weights  used  in  the  computation  of  the  decision  statistic. 


1.5  HOP  FREQUENCY  ESTIMATORS 


As  previously  described,  feature  detectors  focus  on  a  particular  time-domain 
feature  of  the  spread-spectrum  waveform.  In  this  work,  we  develop  a  detector 
based  on  the  feature  hop  frequency.  These  devices  can  do  the  job  of  initial 
detection,  but  we  focus  on  their  estimation  performance.  Two  noteworthy  es¬ 
timators  in  the  literature  fit  this  billing.  The  first6  is  a  maximum  likelihood 
estimator  with  a  structure  similar  to  the  optimal  detector  of  Figure  1.2,  except 
that,  instead  of  summing  the  outputs  of  each  channel,  it  selects  the  channel 

«A.  Polydoros  and  J.K.  Holmes,  "Autocorrelation  Techniques  for  Wideband  Detection  of 
FH/DS  Waveforms  in  Random  Tone  Interference”,  MILCOM  ’83  Conference  Proceedings, 
Boston,  Mass.,  pp.  781-785,  October  17-20,  1983. 

SA.  Polydoros  and  K.T.  Woo,  “LPI  Detection  of  Frequency-Hopping  Signals  using  Auto¬ 
correlation  Techniques”,  IEEE  J.  Select.  Areas  Commun.,  vol.  SAC-3,  no.  5,  September 
1985. 

SN.C.  Beaulieu,  W.L.  Hopkins,  and  P.J.  McLane,  “Interception  of  Frequency  Hopped  Spread 
Spectrum  Signals",  to  be  published  in  IEEE  Trans.  Common. 
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Figure  1.4:  Maximum  Likelihood  Hop  Frequency  Estimator 


with  the  output  of  maximum  magnitude.  The  hop  frequency  corresponding  to 
that  channel  is  declared  the  estimate  (Figure  1.4).  The  second  estimator  of 
note7  is  based  on  the  first  one  but  has  reduced  complexity.  In  this  method, 
wideband  radiometers  cover  the  spread-spectrum  bandwidth  in  order  to  select 
a  small  number  of  subbands  that  most  likely  contain  the  particular  hop.  These 
subbands  are  further  processed  into  fine  bands,  enabling  the  ultimate  selection 
of  the  band  with  the  current  hop. 


1.6  NEW  METHODS  FOR  DETECTION 
AND  HOP  ESTIMATION 

1.6.1  Sequential  Detector 


In  this  work,  a  new  detector  is  developed  based  on  the  ideas  of  sequential 
detection.  It  is  essentially  like  the  optimal  channelized  detector  but,  instead  of 
basing  its  decision  on  accumulated  energy  on  a  predetermined  large  number  of 
hop  dwells,  it  decides,  after  each  hop  dwell,  on  the  presence  or  nonpresence  of 
a  frequency-hopped  waveform.  Because  the  detector  is  based  on  the  Sequential 
Probability  Ratio  Test  (SPRT),  the  test  is  optimal  in  the  sense  that  no  other 

7W.W.  Short  and  R.D.  Chapman,  “Adaptively  Configured  Channelized  Receiver  for  Fre- 
qucncy  Hopped  Signal  Detection  and  Tracking”,  IEEE  International  Conjerence  on  Com¬ 
munications,  1965,  pp.  832-838,  1985. 
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sequential  test  will  make  a  decision  in  less  time  on  average  than  the  SPRT. 
However,  optimality  occurs  only  if  the  signal  is  present  at  a  predetermined 
SNR.  For  smaller  SNRs,  the  SPRT  actually  can  perform  worse  than  a  test 
based  on  a  Fixed  Sample  Size  (FSS).  This  shortcoming  is  averted  by  mixing 
the  SPRT  with  the  FSS  test  to  create  the  Truncated  Sequential  Test  (TST). 
From  these  results  is  designed  an  optimal  test  whose  worst-case  average  decision 
time  is  minimal.  Netted  are  three  new  detectors  that  exploit  the  advantages  of 
sequential  detection:  the  pure  SPRT,  the  TST,  and  the  optimal  TST. 

With  each  of  these  techniques,  the  number  of  samples  needed  for  a  reliable 
decision  is  dramatically  reduced.  One  way  this  detector  performance  gain  can 
yield  a  performance  gain  in  the  interceptor  itself  is  through  decreased  duty 
cycle.  A  typical  interceptor  might  scan  a  particular  direction  in  order  to  de¬ 
termine  the  presence  or  absence  of  communications.  Because  of  the  gains  due 
to  the  sequential  detector,  the  scan  time  is  significantly  reduced.  Another  way 
to  take  advantage  of  the  performance  gain  is  to  add  robustness  to  the  signal 
parameters.  It  was  already  mentioned  that  the  decision  time,  of  the  sequential 
tests  is  dependent  on  the  SNR.  By  hypothesizing  that  the  signal  rests  in  a  band 
of  SNRs,  a  sequential  test  can  be  designed  that  still  outperforms  tests  based  on 
a  fixed  sample  time. 

1.6.2  Compressive- Receiver-Based  Detector  and  Hop 
Estimator 

The  compressive  receiver,  which  simultaneously  estimates  frequency  com¬ 
ponents  over  a  wide,  predetermined  band,  has  promise  as  an  interceptor  with 
both  the  simplicity  of  a  wideband  device  and  the  performance  of  the  channelized 
device.  The  use  of  the  compressive  receiver  for  interception  is  a  largely  unex¬ 
plored  area  with  all  previous  results  being  superficial  and  ad  hoc.  By  contrast, 
two  different  detectors  and  a  hop  frequency  estimator  sire  developed  using  an 
optimal  likelihood  function  approach.  The  first,  the  locally  optimal  detector, 
is  a  detector  with  structure  similar  to  the  channelized  detector  but  operating 
on  the  output  of  the  compressive  receiver.  It  is  locally  optimal,  meaning  that 
for  signals  with  low  SNR  it  gives  the  greatest  probability  of  detection  for  a 
given  probability  of  false  alarm.  Because  the  locally  optimal  detector  has  an 
unwieldy  structure,  it  defeats  the  motivation  to  use  a  compressive  receiver:  sim¬ 
plicity  and  high  performance.  Therefore,  a  time-multiplexed  detector  is  used 
that,  at  the  expense  of  duty  cycle,  can  achieve  performance  as  close  to  optimal 
as  desired.  Both  the  locally  optimal  and  the  time-multiplexed  detector  have 
hop  frequency  estimator  versions.  By  choosing  as  the  hop  frequency  estimate 
the  hop  frequency  corresponding  to  the  detector  channel  with  maximum  out¬ 
put,  a  hop  frequency  estimator  is  formed.  In  conclusion,  two  detectors  and 
a  hop  frequency  estimator  are  developed  with  performance  comparable  to  the 
channelized  devices  but  with  the  simplicity  of  the  broadband  devices. 

1.6.3  Detector  Based  on  the  Amplitude  Distribution 
Function 


A  new  idea  for  detection  is  developed  based  on  the  Amplitude  Distribu¬ 
tion  Function  (ADF).  The  ADF  is  precisely  defined  as  a  function  from  which, 
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through  a  sequence  of  lemmas  and  theorems,  two  results  are  established.  One 
result  is  that  the  ADF  is  roughly  the  average  distribution  of  a  stochastic  process, 
and  the  second  is  that,  for  signal  plus  noise,  the  resulting  ADF  is  the  convolution 
of  the  ADFs  of  signal  and  noise  individually.  The  convolutional  relationship  for 
signal  plus  noise  motivates  the  construction  of  statistical  transform,  called  the 
deconvolution  statistic,  that  converges  to  something  that  is  arbitrarily  close  to 
the  signal  ADF  and  hence  has  potential  for  separating  signal  from  noise  even 
for  low  signal  levels.  How  close  the  deconvolved  ADF  matches  the  signal  ADF 
depends  on  the  proper  choice  of  the  kernel  of  the  deconvolution  statistic.  An  op¬ 
timal  detector  is  presented  that  directly  observes  samples  of  the  deconvolution 
statistic,  yielding  a  test  statistic  of  quadratic  form.  The  ADF-based  detector  is 
a  robust  device  that  is  a  generalization  of  the  radiometer  and  quite  immune  to 
the  details  of  spread-spectrum  modulation. 


1.7  DOCUMENT  ORGANIZATION 


This  work  is  partitioned  into  five  chapters.  The  first  chapter  presents  the 
problem,  precisely  defines  the  type  of  frequency-hopped  waveforms  under  con¬ 
sideration,  and  describes  the  functions  of  the  intercept  receiver  from  a  system 
viewpoint.  It  also  briefly  describes  existing  interception  methods  in  the  cate¬ 
gories  of  initial  detection  and  hop  frequency  estimation  and  then  contrasts  them 
with  the  new  methods  developed. 

The  second  chapter  describes  in  detail  the  new  sequential  detection  meth¬ 
ods.  Within  this  chapter,  the  likelihood  function  for  a  single  epoch  is  developed 
and  asymptotically  analyzed  for  a  large  number  of  hop  frequencies.  Based  on 
this  analysis,  the  synthesis  of  the  FSS  test,  SPRT,  and  TST  are  developed.  Per¬ 
formance  equations  are  presented  along  with  the  results  of  numerical  and  Monte 
Carlo  analyses.  The  optimal  TST  is  described  and  the  asymptotic  efficiencies, 
which  capture  the  low-SNR  test  behavior,  are  presented.  Finally,  conclusions 
are  drawn. 

The  third  chapter  gives  a  detailed  description  of  the  detector  and  hop  fre¬ 
quency  estimator  based  on  a  compressive  receiver.  It  does  this  first  by  precisely 
defining  the  signal  and  compressive-receiver  models  and  then  using  them  to 
develop  equations  for  the  output  signal  component  and  to  characterize  statisti¬ 
cally  the  noise  at  the  compressive- receiver  output.  With  the  detection  problem 
translated  to  the  output  of  the  compressive  receiver,  likelihood  ratio  theory  is 
applied  for  the  low-SNR  case  to  create  the  locally  optimal  detector.  A  sim¬ 
plified  detector,  the  time-multiplexed  detector,  is  also  presented  along  with  a 
hop  frequency  estimator.  All  detectors  are  performance  analyzed  and  numerical 
results  given.  Finally,  conclusions  are  drawn. 

The  fourth  chapter  introduces  the  ADF-based  detector  and  proceeds  with 
an  exposition  of  the  mathematical  tools  developed  for  the  ADF,  which  consist  of 
a  sequence  of  theorems  and  lemmas  culminating  in  a  convolutional  relationship 
between  the  ADFs  of  signal  and  noise.  The  deconvolution  statistic  is  intro¬ 
duced  along  with  family  of  kernels  to  be  used  in  the  statistic.  The  large-time 
statistical  character  of  the  deconvolution  statistic  is  shown  to  be  the  basis  of 
the  ADF-based  detector.  Synthesis  and  performance  analyses  of  the  detector 
are  presented  and  conclusions  are  drawn. 
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The  fifth  chapter  summarizes  the  previous  chapters  and  highlights  the  im¬ 
portant  points.  It  then  suggests  possible  extensions  to  be  investigated  in  the 
future  and  finally  concludes  the  document. 
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CHAPTER  2 

SEQUENTIAL  DETECTION 
METHODS 


2.1  BACKGROUND  AND 
INTRODUCTION 


The  first  task  in  the  interception  of  spread-spectrum  communications  is  the 
detection  of  the  waveform.  This  is  a  prelude  to  other  interception  processes, 
such  as  feature  detection,  channel  tracking,  and  message  extraction.  As  a  new 
development  toward  the  detection  problem,  this  chapter  applies  and  extends 
previously  published  results  in  sequential  detection  to  the  problem  of  the  op¬ 
timal  detection  of  noncoherent  frequency-hopped  (FH)  waveforms.  By  using 
likelihood  function  methods,  the  problem  was  solved8  for  an  FH  waveform  with 
a  known  signal-to-noise  ratio  (SNR)  and  epochs  with  known  starting  times  and 
durations.  However,  in  that  approach,  the  decision  was  based  on  a  data  seg¬ 
ment  of  fixed  size.  Here  a  sequential  approach  is  taken,  meaning  that  whenever 
a  new  data  element  is  collected,  a  decision  about  the  presence  or  nonpresence  of 
an  FH  waveform  is  attempted.  If  no  decision  is  reached,  another  data  element 
is  collected. 

The  sequential  approach  to  detection  has  a  rich  history.  For  the  binary  hy¬ 
pothesis  problem  with  discrete-time  independent  identically  distributed  (i.i.d.) 
data,  Wald9  has  derived  the  optimal  sequential  test.  This  test  is  optimal  in 
the  sense  that  no  other  test  can  reach  a  decision  of  the  same  Neyman-Pearson 
reliability  within  a  shorter  average  time.  This  result  has  been  extended  to  con¬ 
tinuous  time  data. 10,11  Others  have  suggested  tests  that  must  make  a  decision 
within  a  prescribed  time.  These  are  the  “truncated”  tests.12,13,14  Truncation 
is  desirable  not  only  for  implementation  reasons,  but  also  for  improving  the 
performance  of  a  sequential  test  when  the  input  statistics  differ  from  those  as¬ 
sumed  in  designing  the  test.  In  particular,  Tantaratana  and  Poor13  derive  a 

8J.D.  Edell,  “Wideband,  Noncoherent,  Frequency- Hopped  Waveforms  and  their  Hybrids  in 
Low-Probability-of-Intercept  Communications”,  Naval  Research  Laboratory,  Washington, 
D.C.,  Technical  Report  No.  8025,  November  1976. 

9 A.  Wald,  Sequential  Analysis ,  Wiley,  New  York,  1947. 

10N.  Shiryayev,  Optimal  Stopping  Rulcs>  Springer- Verlag,  New  York,  1977. 
nD.A.  Darling  and  A.J.F.  Siegert,  “The  First  Passage  Problem  for  a  Continuous  Markov 
Process”,  Ann.  Math.  Stat.,  vol.  24,  pp.  624-639,  1953. 

12T.W.  Anderson,  “A  Modification  of  the  Sequential  Probability  Ratio  Test  to  Reduce  the 
Sample  Size”,  Ann.  Math.  Stat .,  vol.  31,  pp  165-197,  1960. 

13S.  Tantaratana  and  H.V.  Poor,  “Asymptotic  Efficiencies  of  Truncated  Sequential  Tests”, 
IEEE  Trans.  Inform.  Theory ,  vol.  IT-28,  no.  6,  pp.  911-923,  November  1982. 

1 4 S.  Tantaratana  and  J.B.  Thomas,  “Truncated  Sequential  Probability  Ratio  Test”,  Inform. 
5ci.,  vol.  13,  pp.  283-300,  1977. 
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truncated  sequential  test  for  i.i.d.  Gaussian  data  with  an  unknown  mean,  which 
is  the  foundation  of  the  results  in  this  chapter. 

Development  of  the  sequential  test  is  begun  by  defining  the  observations 
model  for  a  composite  hypothesis  problem.  Specifically,  given  the  observation 
y(t),  the  problem  is  one  of  choosing  between  H o,  which  is  the  hypothesis  that 
an  FH  waveform  is  not  present,  and  //y,  which  is  the  hypothesis  that  an  FH 
waveform  is  present  with  an  SNR  7'  where  0  <  7'.  Exactly,  the  model  is 


H0  :  y(t)  =  n(t) 

vp  rsim 

//y  :  y(t)  =  s(t)  +  n(t) 


0  <  7' 


(2.1) 


where  s(t)  is  given  by  (1.1)  and  n(l)  is  white  Gaussian  noise  with  two-sided 
spectral  density  4^.  The  hypothesized  SNR  7'  is  related  to  the  other  model 
parameters  by  y'  =  S'Th/Np. 

Because  a  reliable  test  cannot  be  devised  for  an  FH  waveform  with  an  arbi¬ 
trarily  small  SNR,  the  preceding  composite  hypothesis  problem  is  simplified  to 
a  binary  hypothesis  problem:  H 0  versus  where  7  is  specified  as  the  smallest 
SNR  that  is  to  be  accurately  detected.  The  quantity  7  =  STh/Np  with  S  being 
the  corresponding  signal  energy.  The  relative  SNR  r  =  \Jy'  fy  is  also  used. 

Using  the  above  observations  model,  the  design  of  a  sequential  test  for  the 
detection  of  FH  signals  is  approached  as  follows.  An  asymptotically  optimal  test 
is  derived  by  applying  the  likelihood  function  theory  to  the  simplified  binary 
hypothesis  problem  Hp  versus  H-y.  The  parameters  of  this  test  are  specified 
to  ensure  a  maximum  probability  of  detection  for  a  given  probability  of  false 
alarm.  This  binary  hypothesis  test  is  then  applied  to  the  more  general  composite 
hypothesis  problem  with  a  resulting  degradation  in  detection  time  that  is  shown 
to  be  controllable  by  properly  truncating  the  test  procedure. 

The  derivation  of  the  asymptotically  optimal  test  begins  with  the  derivation 
of  the  likelihood  function  for  a  single-epoch  observation,  which  is  appropriately 
called  the  Single-Epoch  Likelihood  Function  (SELF).  By  invoking  the  central 
limit  theorem,  Gaussian  densities  are  found  that  are  asymptotic  to  the  actual 
SELF  densities  as  the  number  of  frequencies  becomes  large.  In  determining 
these  densities,  the  SELF’s  means  and  variances  will  be  explicitly  computed 
under  each  hypothesis.  By  next  considering  individual  SELFs  as  the  observa¬ 
tions,  the  problem  will  be  reduced  to  a  binary  hypothesis  problem  with  Gaus¬ 
sian  i.i.d.  observations.  This  simplification  is  justified  because  each  epoch  of 
FH  waveform  has  independent  statistics  and  because  the  SELF’s  statistics  do 
not  depend  on  the  particular  hop  frequency.  Using  these  equivalent  observa¬ 
tions  and  their  asymptotic  densities,  the  Asymptotic  Log-Likelihood  Function 
(ALLF)  is  derived.  The  ALLF  is  then  used  to  synthesize  tests  for  the  binary 
hypothesis  problem.  This  procedure  requires  extending  the  previously  pub¬ 
lished  sequential  tests  to  the  cases  of  data  with  variances  that  depend  on  the 
hypothesis.  Applying  these  results,  a  Fixed-Sample  Size  (FSS)  test,  a  Sequen¬ 
tial  Probability  Ratio  Test  (SPRT),  and  a  Truncated  Sequential  Test  (TST)  are 
designed. 

Each  of  the  three  tests  is  analyzed  by  approximating  the  test  statistic  by  a 
Wiener  process  and  then  employing  the  classical  theory  of  diffusion.11,12  This 
analysis  is  more  general  because  it  yields  the  performance  of  each  test  to  the 
composite  hypothesis  problem  rather  than  just  the  binary  hypothesis  problem 
on  which  the  tests  are  based.  This  analysis  yields  the  average  decision  time  of 
each  test  as  a  function  of  the  input  SNR,  as  well  as  the  operating  characteristic 
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of  each  test.  From  these  results  comes  an  optimal  test,  whose  worst-case  av¬ 
erage  decision  time  is  minimal.  Finally,  a  computer  simulation  confirms  these 
analytical  results. 

To  further  extend  these  results  to  the  case  of  a  test  that  was  synthesized 
under  the  expectation  of  detecting  an  FH  waveform  with  extremely  small  SNR, 
an  asymptotic  analysis  of  a  different  sort  is  undertaken.  This  analysis  shows 
how  the  above  tests  perform  for  the  composite  hypothesis  problem  as  the  min¬ 
imum  reliably  detectable  SNR  of  the  FH  waveform  becomes  increasingly  small. 
Numerical  results  for  this  case  are  given,  but  a  corresponding  computer  simu¬ 
lation  is  not  possible  due  to  the  rate  of  increase  of  the  number  of  computations 
required  as  the  SNR  diminishes. 


2.2  LIKELIHOOD  FUNCTION:  ONE 
EPOCH 


The  statistical  test  for  the  composite  hypothesis  problem  is  defined  by  find¬ 
ing  an  asymptotically  optimal  test  for  a  binary  hypothesis  problem  and  applying 
that  test  to  the  composite  case  and  accepting  the  resulting  degradation.  This 
simplified  binary  problem  consists  of  the  two  hypotheses  Ho,  where  no  signal  is 
present,  and  Hy,  where  a  signal  is  present  with  SNR  7.  For  this  binary  hypoth¬ 
esis  problem,  Appendix  2. A  contains  a  derivation  of  the  SELF,  which  is  the 
likelihood  function  A,  of  the  ith-epoch  observation  y(t)  for  *7/,  <  t  <  (i  +  1)7), . 
The  SELF  is  expressed  as 


A. (2/) 

=  MAi(y/*)] 

(2.2) 

k  =  0  '  ' 

(2.3) 

where  /0  is  the  zeroth  order  modified  Bess^.  function  of  the  first  kind  and 

2  /•('+  l)Th 

Hk 

—  /  y(t)  cos  a/*  t  dt 

jN0Th  J%Th  ’ 

(2.4) 

Qk 

-  ru-^r-  /  y(<)sin  dt. 

vNo Th  J iTs 

Because  of  the  statistical  independence  of  their  respective  observations,  the 
likelihood  function  of  the  n-epoch  observation  is  then  n«  =  i  Aj,  i .e. ,  the  product 
of  these  individual  SELFs. 

The  SELF  is  nicely  modeled  as  the  configuration  of  well-known  devices,  as 
indicated  in  Figure  2.1.  That  is,  the  SELF  is  channelized  so  that  each  channel 
has  a  matched  filter  that  is  tuned  to  a  particular  hop  frequency  and  whose 
output  is  envelope  detected  and  emphasized  by  a  Bessel  function  non-linearity. 
The  output  of  each  channel,  after  scaling  by  C~y  /  K ,  is  summed  to  produce  the 
SELF. 
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Matched  Envelope  Emphasizing 

Filter  Detector  Function 


Figure  2.1:  Block  Diagram  of  Single-Epoch  Likelihood  Function 


2.3  ASYMPTOTIC  LOG-LIKELIHOOD 
FUNCTION 

The  Asymptotic  Log-Likelihood  Function  (ALLF)  is  asymptotic  to  the  n-epoch 
likelihood  function,  11"=  1  as  the  number  of  FH  channels  becomes  large.  The 

critical  idea  behind  the  derivation  of  the  ALLF  is  the  application  of  the  central 
limit  theorem  to  yield  asymptotic  densities  for  the  SELF  from  which,  using  an 
n-epoch  collection  of  SELFs  as  an  equivalent  observation  set,  the  ALLF  will  be 
determined. 

The  SELF  (2.3)  was  computed  assuming  a  binary  hypothesis  problem,  i.e., 
Ho  is  the  hypothesis  that  no  FH  waveform  is  present,  while  //7  is  the  hypothesis 
that  an  FH  waveform  exists  with  a  known  SNR  7.  The  following  analysis  will 
assume  that  an  FH  waveform,  if  present,  will  have  a  SNR  7'  or  equivalently 
an  average  signal  energy  S'  that  is  not  necessarily  equal  to  the  average  signal 
energy  5  assumed  known  in  the  binary  case.  This  generalization  is  not  necessary 
for  deriving  the  ALLF  but  will  be  needed  to  analyze  the  performance  of  the 
ALLF  in  the  composite  hypothesis  problem. 

Proceeding  with  the  derivation  of  the  ALLF,  the  central  limit  theorem  is 
applied  to  the  SELF  to  obtain  an  asymptotic  density  under  all  hypotheses, 

0  <  7'.  The  central  limit  theorem  is  justified  here  because  the  SELF’s  output 
is  the  sum  of  many  channels  whose  statistics  will  be  shown  to  be  nearly  inde¬ 
pendent  and  nearly  identical.  It  will  be  shown  that  the  degree  of  dependence 
between  channels  is  determined  by  the  amount  of  isolation  between  channels, 
which  is  perfect  for  minimally  spaced  channels  as  is  the  case  assumed  here.  It 
also  will  be  shown  that  the  channel  means  and  variances,  while  different  for  the 
signal-present  and  signal-absent  cases,  are  of  a  commensurate  magnitude. 
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2.3.1  Matched  Filter  Output  Statistics 


Because  the  central  limit  theorem  requires  only  the  mean  and  variance  of 
each  channel,  only  the  statistics  of  the  matched  filter  output  need  be  deter¬ 
mined  exactly  since  the  SELF’S  mean  and  variances  can  be  determined  from 
these  statistics  alone.  Assuming  that  the  signal  present  is  in  the  ifcth  channel, 
then  the  matched  filter  output  in  the  /th  channel  can  be  found  from  (2.4)  as 


P,  ~ 


v/2Vsin  0  +  i/i 
v\ 


for  /  =  k 
for  /  /  k 


(2.5) 


f  \/2y  cos  0  +  £/  for  l  =  k 
\  b  for  l  ^  k 


(2.6) 


where 

a  2  /('+|)T‘ 

v\  -  ~r xFZF  n(t)cosu,tdt 

v  AoT\ 

2  /('+UA 

it  =  ~nrfW  /  n(f)sinuq<  dt. 

JiTt, 


(2.7) 


The  matched  filter  outputs  for  the  no-signal-present  hypothesis  //0  are  the 
special  case  of  the  above  expressions  for  7'  =  0.  Two  assumptions  were  made  in 
determining  these  approximate  expressions  for  the  matched  filter  outputs.  The 
first  assumption  is  that  a/* 7\  is  large  and  equivalent  to  requiring  a  large  number 
of  carrier  cycles  over  a  single  epoch.  The  second  assumption  of  orthogonally 
spaced  channels  [i.e.  (w*  ~u>i)Th/2w  is  an  integer]  implies,  in  essence,  that  the 
channels  are  isolated  from  one  another.  Another  condition  implying  channel 
isolation  is  wide  spacing  between  the  channels  [i.e.  (w*  —  uq)T/,  is  large].  In  a 
practical  implementation,  smooth  window  functions  also  could  have  been  used 
in  the  matched  filter  implementation  to  achieve  the  channel  isolation  assumed 
here. 

Simplified  expressions  for  the  matched  filter  outputs  are  represented  by  (2.5) 
and  (2.6).  The  statistical  nature  of  their  noise  components,  {iq}  and  {£/},  is 
determined  next.  From  (2.7),  it  follows  that  the  random  variables  {iq},{£;} 
are  Gaussian  with  zero  mean  and  unity  variance.  Under  the  isolated  channel 
assumption,  it  is  easy  to  show  that 


f[t'mi'n]  =  0  for  m  ^  n  1  <  m,  n  <  K 

£["m£r.]  —  0  for  all  m,n  \  <  m,n  <  K  (2.8) 

£[i|m£n]  =  0  former*  l<m,n<K. 


Thus  {iq},  {£(}  are  mutually  independent,  since  they  are  Gaussian.  These  rela¬ 
tions  also  determine  the  joint  density  of  tq  and  as 

(2.9) 

The  equations  (2.5)  and  (2.6)  and  (2.8),  along  with  the  joint  density  of  iq  and 
(2.9),  constitute  a  complete  statistical  description  of  the  matched  filter  outputs 
{P(}  and  {<?(}. 
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2.3.2  SELF  Moments 


The  statistics  of  {Py}  and  {(?/}  were  found  in  order  to  determine  the  mean 
and  variance  of  the  SELF  (2.3).  The  SELF  moments  are  needed  to  apply  the 
central  limit  theorem  and  thus  ultimately  produce  the  ALLF.  A  few  conditions 
for  the  application  of  the  central  limit  theorem  will  be  established  now.  First, 
since  the  random  variables  {P|}  and  {(?/}  are  mutually  independent,  each  chan¬ 
nel  output  of  the  SELF  is  also  independent.  Furthermore,  the  channel  outputs 
are  all  identically  distributed  except  for  the  output  of  the  channel  with  the 
signal  present.  This  particular  channel  output  will  be  shown  to  have  a  variance 
comparable  to  that  of  the  other  channel  outputs  and  thus  the  central  limit  the¬ 
orem  still  applies  and  with  it  we  get  a  density  asymptotic  to  the  actual  SELF 
density. 

To  continue,  we  need  explicit  expressions  for  the  mean  and  variance  of  the 
SELF.  Assuming  a  signal  is  present  with  a  relative  SNR  of  r  =  then 

the  matched  filter  outputs  of  the  channel  containing  the  signal  are  by  (2.5)  and 

(2.6) 


P,  =  \/Jy  sin  0  +  ui 
Qi  =  x/Iy7  cosO  +  £i . 


(2.10) 


If  fir  and  <x2  are  defined  to  be  the  mean  and  variance  for  this  channel  output, 
then  (2.9)  implies 


Pr  =  £  [/0  (v'FSV'f  +  4?) 


(2.11) 


=  ~  f  '  (2. 


12) 


With  the  rectangular-to-polar  conversion,  Pi  =  pcosip,  Qi  =  psiniji,  and  apply¬ 
ing  the  identity 

1  /2* 

/„(«)=—/  earo">  cb  (2.13) 

2*  Jo 

the  integral  becomes 

^r  =  c-y  /  Pi0(V^ip)  i0{VWp)e-^dP  (2.14) 

JO 


=  e1 1 


o  (2n/^7) 


(2.15) 


This  integral  was  evaluated  as  (1)  in  Section  13.31  of  Watson.15  The  variance 
is  now  evaluated  as  follows: 


*?  +  /*?  =  e[n Kv/W/TTo?) 


(2.16) 


=  ^  j"  j"  l20  (y/Tiy/f’*  +  (2.17) 


15G.N.  Watson,  A  Treatise  on  tie  Theory  of  Bessel  Functions,  Cambridge  University  Press, 
New  York,  1980. 
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which  becomes  with  rectangular-to-polar  conversion,  and  applying  (2.13), 

<x?  +  p2  =  e-7'  f  pi]  (y/Iyp)  I0  {VWp)  e-'rdp.  (2.18) 

Jo 

This  integral  is  evaluated  by  applying  (16)  from  Section  11.41  of  Watson15, 
which  states 

l  Jq  lo  (Va2  +b2  -2a6cos'tf)  d<t>  =  70(a)/0(6).  (2.19) 

Application  of  this  formula  and  an  interchange  of  integrations  reduces  the  in¬ 
tegral  (2.18)  to  a  simpler  integral  solved  as  (1)  in  Section  13.31  of  Watson.15 
The  net  result  is 


a2 +  p2r  = 


e27  r 


jT  e-2aco«0/(i  (4v/9ysin  d<t>.  (2.20) 


Summarizing,  for  a  signal  in  channel  /  with  a  SNR  y1 ,  the  channel  moments  are 


fir  =  C  Io(2ry) 


J2  =  e27  ijf  e-27co^/o(4r7sinf)^-/2(2r7)  . 


(2.21) 


where  r  =  \Jy‘ /y. 

The  above  calculations  give  expressions  for  the  channel  moments  for  a  chan¬ 
nel  with  a  signal  present.  The  moments  for  the  case  of  a  channel  without  a  signal 
present  are  special  cases  of  the  above  with  r  =  0  and  are  thus  denoted  by  po 
and  <7q.  From  (2.21)  and  the  Bessel  function  identity  (2.13),  they  are 


Po  =  e7 

a2  =  e27  [/o(27)  -  1] . 


(2.22) 


Likewise,  moments  for  the  /th  channel  (whenever  it  contains  a  signal  with 
strength  7)  correspond  to  the  above  moments  with  r  =  1  and  are  thus  de¬ 
noted  by  pi  and  crj. 

As  previously  mentioned,  the  application  of  the  central  limit  theorem  de¬ 
pends  on  the  various  channel  means  and  variances  having  commensurate  ampli¬ 
tudes.  The  relative  amplitudes  between  the  moments  are  computed  from  (2.21) 
and  (2.22)  and  are 


—  =  ^o(2r7) 

Vo 


(2.23) 


-j  =  1  +  2r27  +  2r272  +  C)(y3). 


(2.24) 


Hence,  for  small  assumed  SNRs  (7  <  1),  the  mean  and  variance  of  the  channel 
with  a  signal  present  and  the  mean  and  variance  of  the  channels  without  a 
signal  present  are  within  a  factor  of  three  of  each  other. 
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Expressions  for  the  mean  and  variance  of  the  SELF  are  now  immediate, 
since  the  SELF  is  the  sum  of  ail  K  channels  scaled  by  C*1  /  K .  The  expressions 
are 


Mr  =  ^[(A  -  l)/i0  +  /ir]  (2.25) 

A 

K  =  ^[(A'-1K’  +  ^].  (2.26) 

Here,  Mr  is  the  mean  of  the  SELF  when  a  signal  of  strength  7'  is  present  and 
Af0  and  My  are  written  for  the  special  cases  of  Afr,  when  r  =  0  and  r  =  1, 
respectively.  The  variances  Vr ,  Vo,  and  Vi  are  defined  similarly. 

2.3.3  Derivation  of  the  ALLF 


With  the  first  two  moments  of  the  SELF  determined,  the  central  limit  the¬ 
orem  gives  approximating  densities  to  the  SELF,  A,,  under  the  composite  hy¬ 
pothesis  problem.  These  densities  are 


H0  :  A, 

versus 

Hy  :  A, 


(A.  —  Mn)a 
2  Hi 


(2.27) 

for  0  <  7'  <  7 


which  give  a  simplified  statistical  characterization  of  the  SELF.  That  is,  the 
SELF  outputs,  { A; } ,  are  Gaussian  i.i.d.  variables  whose  means  and  variances 
depend  on  the  hypothesis. 

As  was  the  procedure  in  deriving  the  SELF,  the  Asymptotic  Log-Likelihood 
Function  (ALLF)  is  designed  using  the  simpler  binary  hypothesis  problem.  For 
a  single-epoch,  likelihood  function  theory  and  (2.27)  imply  a  log-likelihood  func¬ 
tion  of 

i,(A,)  =  c2A,2  +  c1A,  +c0  (2.28) 

where 


Independence  between  observations  over  different  epochs  implies  that  the  ALLF 
up  to  time  n  is 

n 

=  (2.32) 

1=1 

Now  that  ALLF  has  been  found,  its  mean  and  variance  will  be  computed  as  a 
prelude  to  investigating  its  performance  in  the  composite  hypothesis  problem. 
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2.3.4  Moments  of  the  ALLF 


For  the  analysis  that  follows,  it  is  useful  to  derive  the  moments  of  L,  from 
which  the  ALLF  moments  follow  trivially  from  (2.32),  starting  with  the  mean 

Mr  =  E[Li(\i)]  (2.33) 

=  c2E(A?)  +  ci£(Ai)  +  c0  (2.34) 

=  C2  (Mr  +Vr)+ClMr  +  C0  (2.35) 


which  expands  in  terms  of  the  SELF  moments  to 

Mr  =  \  In  y  +  (2V,  Vo)-1  [(Mr  -  M0)2  V,  -  (Mr  -  M,)2  V0  +  (V,  -  V0)  Vr]  . 

(2.36) 

Now  to  compute  the  variance  of  L,. 

Vr  —  V CLT  [Li  (A,)]  (2.37) 

which  upon  substitution  of  (2.28)  yields 

Vr  =  Var  [c2A2  + 1,  Ai  +  c0]  (2.38) 


=  Var  [(c2M2  +  c,Mr  +  c0)  +  (2c2Mr  +  c,)v  +  c2t>2]  (2.39) 

where  v  =  A,  —  Mr.  Proceeding, 

Vr  =  Var[(2c2Mr +c,)u  +  c2u2]  (2.40) 


=  (2c2Mr  +  c:)2Vr  +  2c2Vr2 


which  simplifies  to 


Vr 


YL(L 

2  VVo 


(2.41) 


(2.42) 


The  special  cases,  r  =  1  and  r  =  0,  of  the  moments  of  L,  are  respectively 
written  as  M\  and  Ma  for  the  means  and  as  V,  and  Vo  for  the  variances. 


2.3.5  Summary 


A  log-likelihood  function  for  the  binary  hypothesis  problem,  designated  the 
ALLF,  has  been  derived  that  is  asymptotic  to  the  true  log-likelihood  function 
as  the  number  of  channels  becomes  large.  The  ALLF  was  found  with  the  help  of 
the  likelihood  function  theory  by  considering  an  n-epoch  collection  of  SELFs  as 
a  set  of  i.i.d.  observations  assumed  Gaussian  by  the  central  limit  theorem.  The 
Gaussian  assumption  was  justified  by  showing  that  each  SELF  was  the  sum  of 
nearly  independent  and  nearly  identical  random  variables.  Various  means  and 
variances  were  also  derived  that  will  prove  useful  in  future  discussions.  The 
ALLF  now  will  be  used  to  design  an  FSS  test,  an  SPRT,  and  a  TST. 
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2.4  TEST  DESIGN 


The  results  above  reduced  the  problem  of  detecting  an  FH  waveform  to  that 
of  discriminating  between  two  sets  of  Gaussian  i.i.d.  data  with  different  means 
and  variances.  A  Fixed  Sample  Size  (FSS)  test,  a  Sequential  Probability  Ratio 
Test  (SPRT),  and  a  Truncated  Sequential  Test  (TST)  based  on  this  simplified 
model  will  be  discussed. 

2.4.1  FSS  Test  Design 


As  the  name  suggests,  an  FSS  test  consists  of  comparing  a  test  statistic  T t, 
based  on  a  fixed  number  of  observations  L,  to  a  threshold  r.  Then,  if  the  test 
statistic  is  greater  than  r,  hypothesis  H\  is  chosen,  while  a  test  statistic  less 
than  r  indicates  hypothesis  Hq.  Symbolically  this  is 

J  r  J  H„.  <2  «> 


In  our  case,  the  test  statistic  is  the  L-epoch  ALLF  and  the  test  parameters  L 
and  r  are  specified  to  correspond  to  prescribed  false  alarm  PF  and  detection 
Pp  probabilities.  To  determine  L  and  r,  the  density  of  the  Ti  is  needed  for 
each  hypothesis.  Although  this  density  equals  the  non-central  \2  density,  an 
approximate  Gaussian  density,  derived  via  the  central  limit  theorem,  is  used 
instead  to  yield  simplified  expressions  for  the  test  parameters.  These  densities 
are 


versus 


H0: 

n  ~ 

{T,.-L  Mo)3 

l  /•>  •Sr.u: 

v,2*tVoC 

Hi  : 

n  ~ 

(TY-LA.,)3 
*  &  3LV. 

/S  Av,  c 

(2.44) 


From  these  densities,  Pp  and  Pf  can  be  computed  in  terms  of  L  and  t  to  yield 
Pf  =  1  —  $  Pd  —  1  —  1 )  where  4>_1  is  the  inverse  c'  the 

distribution  function  of  a  zero-mean,  unity-variance  Gaussian  random  variable. 
These  are  solved  simultaneously  to  arrive  at 


PD)-Vj*-1(l- pF)]7 
______ 


(2.45) 


T 


L* 

(Mi  -  Mo) 


V*Mi<t>-l(  1  -  PF)  -  vtMoQ-'i  1  -  Pd)]  . 


(2.46) 


2.4.2  SPRT  Design 

Wald’s  sequential  probability  ratio  test  (SPRT)  now  can  be  defined  as  a  test 
with  test  statistic  Tn,  based  on  n  observations  and  two  thresholds  a  and  b.  The 
SPRT  works  as  follows.  Upon  the  nth  observation,  if  Tn  is  greater  than  a,  then 
hypothesis  Hi  is  chosen.  If  Tn  is  less  than  b,  then  hypothesis  Ho  is  chosen.  If, 
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instead,  Tn  is  between  a  and  6,  the  test  statistic  is  updated  to  include  n  +  1 
observations  and  the  process  is  iterated.  Symbolically  this  test  is  described  as 

(  >a  =>  Hi 

for  each  n,  Tn  <  <  6  H0  (2.47) 

(_  €  (a,  b)  =>  take  another  sample. 

The  threshold  values  a  and  b  are  assigned  to  give  the  desired  Neyman-Pearson 
probability  of  detection  Po  and  probability  of  false  alarm  Pp.  Relationships  be¬ 
tween  the  thresholds  and  these  probabilities  are  given  by  Wald’s  approximations9 


(2-48) 


2.4.3  TST  Design 


Truncated  Sequential  Test  (TST)  is  a  hybrid  of  the  above  two  tests.  Specif¬ 
ically,  TST  follows  the  rules  of  a  sequential  test  with  test  statistic  Tn  and  with 
thresholds  a  and  6,  but  has  the  added  feature  of  forcing  a  decision  at  time  L 
(if  no  decision  has  been  made)  by  comparing  the  test  statistic  to  a  threshold  r. 
Symbolically, 

(  >  a  =>  Hi 

for  each  n  <  L,  T„  <  <  6  =>  Ho 

l  €(a,6)  =>  take  another  sample 

V  (2.50) 

bul  for  ti  =  L,  Tl{*'t  % 

Two  relations  secure  the  specification  of  the  TST  parameters  a,  6,  L,  and  r. 
If  Pp  and  Pp  are  the  actual  Neyman-Pearson  probabilities  for  the  TST,  then 
from  Tantaratana  and  Poor13 

Pp  <  P£ss  +  P^PRT  (2.51) 


(1  -P’d)  <  (1  -P£ss)  +  (l-PZPRT)  (2.52) 

where  Pfss  is  the  probability  of  false  alarm  for  the  TST,  if  L  —  oo,  and  PpPRT 
is  the  false-alarm  probability  for  the  TST,  if  a  —  —  b  =  oo.  PpSS  and  PpPRT  are 
defined  similarly.  Thus  the  errors  of  the  TST  can  be  viewed  as  a  mixture  of  the 
errors  of  an  FSS  test  with  parameters  L  and  r  and  an  SPRT  with  parameters  a 
and  6.  These  inequalities  can  be  verified  by  viewing  the  ALLF  T„  as  a  discrete 
stochastic  process  with  time  index  n  and  enumerating  its  sample  paths.  For 
instance,  a  sample  path  leading  to  a  false  alarm  must  either  cross  threshold  a 
before  threshold  6  and  before  time  L  or  be  greater  than  threshold  r  at  time  L. 
Since  these  events  also  correspond  to  false  alarms  in  either  the  FSS  test  part 
or  the  SPRT  part  of  the  TST,  the  inequality  (2.51)  must  follow. 

The  above  inequalities  can  be  used  to  specify  a  TST,  whose  actual  error 
probabilities  Pp  and  1  —  Pp  are  less  than  any  specified  error  probabilities  Pp 
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and  1  -  Pd  Thus,  the  TST  can  be  designed  by  partitioning  the  bounding 
errors  (1  —  Pp)  and  Pp  among  the  SPRT  and  FSS  test  parts  of  the  TST  and 
then  using  the  appropriate  equation  to  compute  the  parameters  L,  r,  a,  and  b 
for  TST13.  Specifically,  this  partitioning  is  quantified  with  the  introduction  of 
two  constants,  0  <  C\  <  1  and  0  <  Ci  <  1,  which  are  defined  as  TST  mixture 
constants,  then 


Pfss  =  CiPF 

(2.53) 

pSPRT  _  (1  -  Ci) PF 

(2.54) 

5 

I 

£ 

II 

<0 

V) 

■fc? 

1 

(2.55) 

(1  -Psdprt)  =  (1  —  C2XI  —  Pd) 

(2.56) 

for  the  error  probabilities  of  the  FSS  test  and  SPRT  parts  of  the  TST.  From  the 
above  inequalities  and  (2.45),  (2.46),  (2.48),  and  (2.49),  the  TST  parameters 
are  determined  as  follows: 

-  P£ss)  -  V02<&-1(1  -  Pfss ) 

l  _  i - - - L- 

(Mi -Mo) 

(2.57) 

-  p£ss)\ 

(2.58) 

/ pSPHT \ 
a  ~  *n  pSPRT  J 

(2.59) 

L  ,  fl-PSDPRT\ 

6  ~  "  (l  -  PfPRT  )  ' 

(2.60) 

Note  that  (2.51)  and  (2.52)  guarantee  that  the  actual  detection 

errors 

£ 

VI 

£ 

(2.61) 

1 

a 

IA 

1 

(2.62) 

The  mixture  constants  Ci  and  C2  reflect  proportions  of  the  FSS  test  and  SPRT 
parts  of  the  TST  since,  if  C\  =  C2  =  1,  a  pure  FSS  test  is  defined,  and  if 
Ci  =  C2  =  0,  a  pure  SPRT  is  defined.  Criteria  for  choosing  the  mixture 
constants  will  be  discussed  in  Section  2.6. 


2.5  PERFORMANCE  OF  TESTS 


The  problem  addressed  by  the  preceding  tests  (the  FSS  test,  the  SPRT, 
and  the  TST)  is  the  detection  of  an  FH  waveform.  The  detection  of  the  FH 
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waveform  is  a  prelude  to  other  interception  processes,  such  as  feature  detection, 
channel  tracking,  and  message  extraction.  Here  the  performance  of  the  tests 
in  detecting  an  FH  waveform  with  variable  amplitude  and  in  the  presence  of 
white  Gaussian  noise  is  quantified. 

The  three  tests  were  designed  under  the  assumption  of  binary  hypotheses. 
These  hypotheses  are  Hq  (FH  waveform  is  not  present)  and  H-,  (FH  waveform 
is  present  and  has  SNR  7).  Of  concern  here  is  the  performance  of  the  three 
tests  when  the  actual  SNR  7'  of  the  FH  waveform  is  more  generally  0  <  Y  <  7. 
Two  parameters  characterize  a  test’s  performance  for  a  particular  7'.  The 
first,  denoted  by  E(N/r,  7),  is  the  Average  Sample  Number  (ASN)  defined  as 
the  average  of  the  number  of  samples  needed  to  reach  a  decision.  The  second 
parameter,  denoted  by  P0(r,y),  is  the  Operating  Characteristic  (OC)  defined 
as  the  probability  of  declaring  the  absence  of  an  FH  waveform.  Both  the  ASN 
and  OC  are  defined  as  functions  of  relative  SNR  r  and  the  assumed  SNR  7. 

2.5.1  Analysis  of  FSS  Test 


For  the  FSS  test,  the  ASN  is  obviously  L,  while  the  OC  can  be  determined 
by  approximating  the  ALLF  at  time  L  by  a  Gaussian  random  variable  with  the 
same  moments.  This  central-limit-theorem  argument  produces 

n(r--r)=*(rmL)  (2-63> 

for  the  OC. 


2.5.2  Analysis  of  SPRT 


For  the  SPRT,  the  analysis  is  more  difficult  but  can  be  approached  as  a 
diffusion  problem.  Here  we  approximate  the  test  statistic  by  a  Wiener  process. 
Specifically,  if  T(t)  is  a  Wiener  process  with  variance  function  Vr<  and  mean 
function  Mrt,  then  the  ALLF,  Tn,  converges  weakly  to  T(t)  at  integer  times 
<  =  n,  provided  n  is  sufficiently  large.  This  last  restriction  is  needed  to  ensure 
that  Tn  has  an  approximate  Gaussian  density  as  implied  by  the  central  limit 
theorem.  In  terms  of  the  approximating  Wiener  process  T(t),  the  problem  of 
finding  the  OC  function  is  now  the  problem  of  finding  the  probability  that  T(i ) 
will  “touch”  the  lower  threshold  6  before  the  upper  threshold  a.  Likewise,  the 
problem  of  finding  the  ASN  is  now  the  problem  of  finding  the  average  time 
that  T(t)  first  “touches”  either  threshold  (a  or  b).  This  time  is  also  called  the 
average  stopping  time.  Expressions  for  these  quantities  are  given1112  as 


Po(r,j)  =  < 


S(N/r,  7)  =  { 


e~2hifc  - 1 

_  g  — 2a- 


a  —  b 


Mr  tO 
Mr  =0 


(2.64) 


afoir, 7)  +  t[l  -  Po(r,7)] 
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VV 


Mr  to 


Mr-  0 


.  (2.65) 
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2.5.3  Analysis  of  TST 


The  diffusion  analysis  technique  also  applies  to  the  TST  but  is  more  in¬ 
volved.  The  ASN  is  by  Darling  and  Siegert11 

OO 

E(N/r,7)  =  A  £(- 1)"  ^Bn{e~k'L  -  1)  (2.66) 

nr  1  n 


where 


A 

1! 

'cT 

1  ^ 

to 

(2.67) 

Bn 

M,b  .  nira  Mra  .  nirb 

=  €  vr  sin - -  -e  v'  sin - - 

a  —  6  a  —  b 

(2.68) 

kn 

Ml  Vrn2ir2 

2Vr  +  2(a  -  6)2 ' 

(2.69) 

The  OC  function  is  defined  by  Anderson12  as 

E°°  f  ^  ~  -  2[na  -  (n  -  1)6]  \ 

V  vm  ) 


(n~  t 

_ g2  T7^n((|-a)<J) 


/  r  -  LMr  -  2n(a  -  6)\ 

V  v/IVr  ) 


2[n6  —  (n  —  l)a]  —  r  +  LM 


VZvr 

2n(6  —  a)  —  r  +  LM 


(2.70) 


Equations  (2.63)  through  (2.70)  represent  a  complete  characterization  of  the 
performance  of  the  FSS  test,  the  SPRT,  and  the  TST. 

The  fact  that  the  diffusion  technique  yields  accurate  expressions  for  the 
ASN  and  OC  functions  will  not  be  proved  here  but  will  be  verified  below  by 
computer  simulation. 


2.5.4  Numerical  Results 


The  FSS  test,  the  SPRT,  and  the  TST  were  simulated  by  the  computer  to 
verify  the  assumptions  of  the  analysis  and  as  an  independent  measure  of  the 
relative  performance  of  the  three  tests.  The  simulated  detector  consisted  of  512 
channels  and  each  test  was  synthesized  to  ensure  a  probability  of  false  alarm 
Pf  of  no  more  than  1%  and  a  probability  of  detection  Pd  of  at  least  99%. 


I 
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Here  the  fairly  relaxed  probability  of  false  alarm  of  1%  was  chosen  in  order 
to  limit  the  number  of  data  needed  for  a  decision.  Under  these  specifications, 
the  simulation  was  run  until  1,000  decisions  were  reached  for  each  of  11  SNRs 
evenly  spaced  between  0  and  7.  The  decisions  that  no  FH  waveform  was  present 
were  averaged  to  estimate  the  OC,  while  the  number  of  observations  taken  to 
reach  a  decision  was  averaged  to  estimate  the  ASN.  Additionally,  the  standard 
deviation  of  ASN  average  was  measured  to  indicate  the  ASN  estimation  error. 

Figure  2.2  and  Figure  2.3  are,  respectively,  the  ASNs  by  simulation  and 
by  theory  when  the  assumed  SNR  7=1,  while  Figure  2.4  and  Figure  2.5 
are  the  corresponding  curves  for  (7  =  0.3).  As  predicted,  the  ASN  is  greatly 
reduced,  by  about  57%,  for  the  SPRT  in  the  regions  around  7'  =  0  and  7'  =  7. 
These  curves  exemplify  a  general  property  of  the  SPRT:  to  perform  very  well 
when  the  observation  statistics  are  close  to  those  assumed,  but  to  exhibit  a 
degraded  performance,  often  to  the  point  of  being  worse  than  the  FSS  test, 
when  the  observation  statistics  are  different.  In  our  context,  this  degradation 
is  evidenced  by  a  large  ASN  for  the  SPRT,  when  the  actual  SNR  7'  is  midway 
between  the  two  assumed  values  0  and  7.  The  TST  reduces  the  ASN  in  this 
mid-SNR  region,  as  shown  by  the  figures,  but  it  does  so  at  the  expense  of 
performance  in  the  regions  around  7'  =  0  and  7'  =  7.  Despite  this  performance 
loss,  truncation  is  necessary  for  implementation  reasons.  It  also  will  be  shown 
that  the  TST  has  the  desirable  property  of  having  a  higher  detection  probability 
than  the  SPRT  at  small  SNRs  and  that,  through  optimization  of  the  mixture 
constants,  the  TST  can  regain  much  of  what  it  lost  in  ASN  around  7'  =  0  and 
7'  =  7- 

Focusing  on  the  OCs  (Figures  2.6  and  2.7  for  7  =  1,  Figures  2.8  and  2.9 
for  7  =  0.3) ,it  is  obvious  that  the  FSS  test  has  slightly  higher  probability  of 
detection  for  small  SNRs  while  the  SPRT  has  degraded  performance  in  this 
region.  Notice  that  these  test  performances  are  reversed  for  SNRs  close  to  7. 
The  OCs  also  show  that  the  TST’s  actual  detection  errors,  Pp  and  1  —  Pd,  are 
within  79%  of  their  specified  bounds,  Pp  and  1  —  Pp- 

Throughout  the  analysis,  various  simplifying  approximations  were  made 
whose  accuracies  were  hard  to  quantify,  especially  the  Wiener  process  approx¬ 
imations  to  the  ALLF.  Thus  the  computer  simulation  was  compared  quanti¬ 
tatively  to  results  predicted  by  theory  as  a  validation  of  assumptions  made. 
Table  2.1  for  7  =  1  and  Table  2.2  for  7  =  0.3  show  how  well  the  simulation 
of  the  three  tests  corresponds  to  the  analysis.  The  quantity  AASN  is  the  nor¬ 
malized  difference  between  the  theoretical  ASN  and  the  simulation  ASN,  where 
the  normalizing  factor  is  the  estimated  standard  deviation  of  the  average  used 
to  estimate  the  ASN.  The  AASN  values  show  a  good  correspondence  between 
theory  and  simulation,  since  they  are  within  two  standard  deviations  86%  of 
the  time.  The  quantity  AOC  is  the  normalized  difference  between  the  theo¬ 
retical  OC  and  the  simulation  OC.  Here  the  normalizing  factor  is  the  standard 
deviation  of  the  OC  average,  assuming  that  the  theoretical  OC  value  is  correct. 
In  other  words,  the  normalizing  factor  for  a  theoretical  OC  of  Po(y')  and  1,000 
simulation  runs  is  <r oc  =  ~  /VtOI/IOOO.  Here  again,  in  Tables  2.1 

and  2.2,  a  good  correspondence  between  theory  and  simulation  is  apparent. 

The  purpose  of  the  computer  simulation  was  to  validate  the  assumptions 
made  in  the  specification  and  analysis  of  the  three  tests:  the  FSS  test,  the 
SPRT,  and  the  TST.  The  accuracy  with  which  the  analysis  predicts  quantities 
measured  by  simulation,  as  shown  above,  substantiates  the  assumptions  made. 
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Figure  2.2:  ASN  from  Simulation  versus  SNR,  7=1 
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Figure  2.3:  ASN  from  Theory  versus  SNR,  7=1 
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Figure  2.4.  ASN  from  Simulation  versus  SNR,  y  =  0.3 
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Figure  2.5:  ASN  from  Theory  versus  SNR,  7  =  0.3 
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Figure  2.6:  OC  from  Simulation,  7  =  1 
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Figure  2.7:  OC  from  Theory,  7  =  1 
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Number  of  Channels  =  512 


Figure  2.8:  OC  from  Simulation,  7  =  0.3 
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Figure  2.9:  OC  from  Theory,  7  =  0.3 
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Table  2.1:  Comparison  Between  Theory  and  Simulation  for  7=1 


■1 

FSS 

test 

SPRT 

TST 

II 

A  ASN 

A  OC 

A  ASN 

A  OC 

A  ASN 

A  OC 

0.00 

0.00 

■iWSM 

-1.19 

-0.03 

-1.47 

0.03 

0.10 

0.00 

1.11 

-0.48 

mWFm 

-0.76 

2.03 

0.20 

0.45 

-1.05 

■WSEM 

-1.84 

1.02 

0.30 

0.00 

0.18 

1.16 

-0.09 

1.19 

-0.18 

0.40 

0.00 

-0  25 

1.09 

0.89 

0.20 

-0.40 

0.50 

I  Bum 

0.73 

0.60 

0.83 

-0.64 

1.53 

0.60 

TiTM 

0.22 

0.23 

0.02 

-0.14 

-0.15 

0.70 

0.00 

ItilKlrdh 

-1.47 

1.64 

-2.35 

3.20 

0.80 

0.00 

-2.16 

-0.94 

2.25 

0.76 

1.17 

0.90 

— 1.33 

-1.38 

0.83 

-1.88 

1.02 

1.00 

j  0.00 

-0.84 

-1.16 

3.50 

-1.38 

1.29 

Table  2.2:  Comparison  Between  Theory  and  Simulation  for  7=0.3 


FSS  test 


SPRT 


HI 

A  ASN 

A  OC 

A  ASN 

A  OC 

A  ASN 

A  OC 

0.00 

0.00 

-0.71 

-1.92 

0.02 

-1.04 

2.09 

IliUKlI 

0.00 

-0.73 

0.32 

-1.75 

1.38 

-2.18 

0.06 

0.34 

1.59 

-0.82 

-0.71 

-1.03 

0.09 

0.00 

0.68 

-0.65 

-0.49 

0.07 

0.28 

IQHI 

0.00 

-0.24 
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2.6  TEST  EXTENSIONS 


The  analytical  expressions  for  the  ASN  and  the  OC  of  the  TST,  (2.66  and 
2.70),  can  be  used  to  determine  a  TST  with  an  optimum  mixture  of  FSS  and 
SPRT  parts.  Specifically,  the  maximum  ASN  with  respect  to  the  SNR  7'  varies 
as  a  function  of  the  mixture  constants  C\  and  Cx-  This  function  is  graphed  in 
Figure  2.10.  The  figure  indicates  that  the  optimal  TST  should  have  a  greater 
mix  of  SPRT  than  the  value  one-half  used  in  Section  2.5,  since  the  maximum 
ASN  of  smallest  value  occurs  for  smaller  values  of  the  mixture  constants  C\  and 
C-x.  The  optimal  mixture  constants  were  found  numerically  to  be  C 1  =  0.286 
and  Cx  —  0.284.  The  ASN  and  OC  of  the  optimal  TST  are  shown  in  Figures  2.11 
and  2.12.  It  is  interesting  that,  by  minimizing  the  maximum  ASN,  the  ASN  in 
the  extreme  regions  about  7'  =  0  and  7'  =  7  is  also  reduced.  This  is  believed 
to  be  a  consequence  of  the  optimal  TST  having  a  greater  SPRT  mix  than  the 
half-and-half  arbitrarily  picked  for  the  Section  2.5  simulation  and,  therefore, 
exhibiting  properties  closer  to  a  pure  SPRT.  Of  course,  if  the  first  TST  was 
specified  to  have  a  larger  SPRT  mix,  then  optimization  would  have  increased 
the  ASN  in  the  extreme  regions.  The  optimal  TST  offers  a  good  compromise 
between  the  need  to  maximize  the  ASN  performance  in  the  extreme  regions 
and  the  need  to  minimize  the  maximum  ASN. 

Another  extension  to  the  previously  described  tests  is  the  robustification  of 
the  tests  with  respect  to  the  assumed  SNR  7.  This  can  be  accomplished  by 
specifying  the  assumed  SNR  7  as  the  worst  case,  and  then  choosing  a  corre¬ 
sponding  minimum  probability  of  detection  Pq  that  is  somewhat  relaxed.  This 
procedure  produces  a  test  that  adequately  detects  over  a  broader  range  of  SNRs, 
and  it  is  a  way  to  effectively  use  the  smaller  detection  times  of  the  SPRT  and 
TST.  In  this  way,  either  a  TST  or  an  SPRT  can  be  designed  to  adequately 
detect  over  a  broader  range  of  SNRs  than  an  FSS  test  with  the  same  or  greater 
detection  time. 

The  described  tests  also  can  be  extended  to  the  slow  FH  case.  The  detector 
structure  itself  is  optimal  under  the  fast  FH  assumption  but  is  also  a  reasonable 
suboptimal  structure  for  slow  FH  signals.  This  is  especially  true  when  there  are 
a  large  number  of  hops  over  a  given  detection  time.  Even  though  the  detector 
itself  is  suboptimal  for  slow  FH,  all  the  performance  and  design  analysis  devel¬ 
oped  for  fast  FH  applies  directly.  This  is  because  all  design  and  performance 
analysis  depended  only  on  the  chip  duration  being  known  and  the  interference 
being  additive  white  Gaussian  noise. 


2.7  ASYMPTOTIC  EFFICIENCIES 

The  previous  analysis  did  not  include  the  performance  of  the  tests  when  the 
assumed  SNR  7  is  small.  This  case  will  be  examined  here.  Since  the  ASN  and 
OC  are  functions  of  both  7  and  the  actual  SNR  7',  the  ASN  and  OC  can  be 
recast  as  functions  of  7  and  the  relative  SNR  r  =  \f~t' /~t-  Test  performance 
in  the  dwindling  SNR  case  is  captured  by  the  limit  of  the  ASN  and  OC,  as  7 
diminishes  while  r  is  held  constant.  For  the  OC,  this  is  a  finite  limit,  but  the 
ASN  increases  without  bound.  Thus,  rather  than  comparing  the  ASNs  directly, 
the  limit  of  the  ASN  times  y2  is  computed.  In  other  words,  a  quantity,  identified 
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Figure  2.11:  ASN  of  Optimal  TST  versus  SPET  and  half-mixed  TST 
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Figure  2.12:  OC  of  Optimal  TST  versus  SPRT  and  half-mixed  TST 
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as  the  asymptotic  ASN  E(N/r ),  is  defined  as 

£(N/r)=  limy2£(N/r,y). 


(2.71) 


The  asymptotic  ASN  is  useful  because  it  preserves  the  relative  efficiencies  be¬ 
tween  the  ASNs  as  y  diminishes.  For  instance,  consider  the  FSS  test  ASN, 
£rss(N/r,  7),  and  the  SPRT  ASN,  ESPKT(N/r, 7)  and  write 


eFSS{N/r,y)  _  tFSS(N/r) 
7™  SSPRT{N/r,y) 


-SPRT .  .  ~ 

S  (N/r) 


(2.72) 


where  EFSS(N/r)  and  ESPRT(N/r)  are  the  asymptotic  ASNs  of  the  FSS  test 
and  SPRT,  respectively.  The  asymptotic  OC  is  simply  defined  as 

P0(r)  =  lim  Po(r,y).  (2.73) 

7-.0 


As  an  aid  in  evaluating  these  limits,  asymptotic  expressions  for  moments  of 
the  single-epoch  ALLF  derived  in  Appendix  2.B  are  defined  as  follows: 


Mr  =  Adr72  +  0(  73) 

(2.74) 

Vr  =  Vr72  +  C>(73) 

(2.75) 

where 

-  K+ 2  2  1 

Adr  =  K.  (r2  -) 

(2.76) 

Vr  -  A+2 

K  2 

(2.77) 

Throughout  this  discussion,  the  quantity  0( 7")  represents  any  function,  say 

/( 7),  such  that 

y~n  lim  /( 7)  <  00. 

7—0 

(2.78) 

The  particular  function  represented  by  0(yn)  is  determined  from  the  context 
of  the  equation  in  which  it  appears. 

To  ease  the  expression  of  the  asymptotic  ASN  and  asymptotic  OC,  the  vari¬ 
ables  L,  f,  a,  and  b  are  defined.  They  will  be  labeled  the  asymptotic  test 
parameters.  Depending  on  the  test  type,  they  have  expressions  that  corre¬ 
spond  to  that  test  type’s  parameter  equations,  where  Adr  is  replaced  with  Adr 
and,  likewise,  Vr  is  replaced  with  Vr.  For  instance,  the  FSS  asymptotic  test 
parameters  from  (2.45)  and  (2.46)  are 

-  pd)  -  v^-Hi  -  Pf)]2 

L  =  i - - - -2 - -  (2.79) 

(Adi  -  Ado) 

f  =  -  L-  ----  tof  Adi$-‘(1  ~  PF)  -  Adofl'^l  -  PD)]  . 

(Mi  —  Mq)  l  j 


(2.80) 


51 


The  Johns  Hopkins  University 

Applied  Physics  Laboratory 


By  using  the  asymptotic  expressions  (2.74)  and  (2.75),  we  have  proved  that  the 
asymptotic  ASN  and  OC  of  a  particular  test  are  exactly  those  of  a  test  with  the 
corresponding  asymptotic  test  parameters.  For  example,  this  fact  and  (2.63) 
suggest  that  for  the  FSS  test  the  asymptotic  ASN  is 

EFSS(N/r )  =  L  (2.81) 

while  the  FSS  test’s  asymptotic  OC  is 

P0(r)  =  *  ("7 Lvf )  '  (2'82) 

The  ASN  and  OC  for  the  three  different  tests  are  plotted  and  compared  in 
Figures  2.13  rind  2.14.  The  relative  relationship  among  the  tests’  asymptotic 
ASNs  is  almost  exactly  like  that  between  the  ASNs  for  7  =  1  and  7  =  0.3,  shown 
in  Figures  2.2  through  2.5.  This  indicates  that  the  three  tests  have  reached 
their  asymptotes,  even  for  7  =  1.  This  comment  also  applies  to  the  OCs.  The 
usefulness  of  this  asymptotic  analysis,  beyond  verifying  that  the  relationship 
among  tests  remains  the  same  for  diminishing  SNR,  is  that  it  simplifies  the  test 
parameter  relationships  with  respect  to  the  parameters  7  and  r.  Thus  for  each 
test,  we  could  choose  parameters  L  =  Ly~7,  r  =  f,  a  =  a,  and  6  =  6  and  have 
comparable  performance  for  any  small  7.  This  feature  simplifies  any  adaptation 
with  respect  to  7  that  might  be  added  to  this  detection  scheme. 

2.8  CONCLUSIONS 


Methods  for  the  sequential  detection  of  noncoherent  fast  FH  waveforms  have 
been  developed.  In  the  process,  the  FH  waveform  is  modeled  to  have  an  in¬ 
formation  component  that  consisted  of  a  series  of  chips  with  a  known  constant 
epoch,  where  each  chip  frequency  is  one  of  a  known  ensemble  of  frequencies. 
In  the  model,  a  particular  chip  frequency  is  independently  determined  by  a 
uniform  random  variable  on  the  frequency  ensemble.  The  F1I  waveform  is  as¬ 
sumed  to  have  an  additive  white-noise  component.  By  assuming  the  modeled 
FH  waveform  to  be  a  known  SNR,  the  optimal  detector  based  on  a  single-epoch 
observation  (SELF)  is  developed  using  likelihood-function  theory.  SELF  is  the 
sum  of  many  nearly  identical  and  nearly  independent  random  variables  and  thus 
has  nearly  Gaussian  statistics.  This  central-limit  argument  allows  a  multi-epoch 
collection  of  SELFs  to  be  considered  an  equivalent  set  of  Gaussian  i.i.d.  vari¬ 
ables.  From  these  simplified  observations,  a  log-likelihood  function  (ALLF)  is 
computed  that  is  asymptotic  to  the  exact  log-likelihood  function  as  the  number 
of  possible  hop  frequencies  becomes  large.  The  ALLF  becomes  the  test  statistic 
on  which  three  detection  tests  are  based.  The  tests  are  the  FSS  test,  SPRT, 
and  TST.  These  are  defined  to  ensure  that  detection  errors  are  below  desired 
levels.  By  modeling  the  ALLF  as  a  Wiener  process,  diffusion  theory  yields  the 
performance  of  the  three  tests  not  only  for  an  FH  waveform  of  the  assumed  SNR 
but  also  for  all  SNRs  below  the  one  assumed.  This  analysis  compares  favorably 
with  a  computer  simulation  of  the  detector  that  validates  the  analysis.  The 
analysis  also  becomes  a  tool  used  to  numerically  optimize  performance  of  the 
TST,  when  the  actual  FH  SNR  deviates  from  that  assumed.  In  order  to  study 
the  performance  of  tests  synthesized  by  assuming  an  extremely  small  FH  SNR, 
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Figure  2.13:  Asymptotic  ASN  of  FSS  Test,  SPRT,  and  TST 
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Figure  2.14:  Asymptotic  OC  Function  of  FSS  Test,  SPRT,  and  TST 
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expressions  for  the  asymptotic  test  efficiencies  are  computed.  This  asymptotic 
analysis  also  yields  simplified  test  parameter  expressions  applicable  to  the  small 
SNR  case. 

A  significant  feature  of  the  SPRT  exposed  by  the  analysis  is  that,  with  the 
same  error  probabilities,  an  FH  waveform  with  a  given  SNR  can  be  detected 
in  less  than  half  the  time  of  the  corresponding  FSS  test.  This  reduction  in 
detection  time  is  especially  significant  for  Low  Probability  of  Intercept  (LPI) 
applications  in  which  the  transmissions  are  purposely  short.  For  the  pure  SPRT, 
detection  time  increases  whenever  the  observed  SNR  differs  from  that  assumed 
in  the  test’s  synthesis.  For  SNRs  midway  between  zero  and  the  assumed  value, 
it  is  even  comparable  to  the  corresponding  FSS  test.  The  TST  significantly 
improves  this  anomaly  while  sacrificing  little  performance  over  that  of  the  purely 
sequential  test,  and  what  little  performance  is  lost  is  largely  regained  by  the 
optimal  TST.  The  decrease  in  the  detection  time  of  the  sequential  tests  can 
be  used  to  robustify  the  test  with  respect  to  the  input  SNR,  while  maintaining 
better  performance  than  that  of  the  non-robust  FSS  test.  The  simplified  test 
parameter  expressions  derived  by  asymptotic  methods  may  be  useful  for  any 
schemes  to  adapt  these  tests  for  varying  FH  SNRs.  The  three  tests  and  their 
corresponding  design  and  performance  analysis  also  apply  to  the  slow  FH  case. 
The  detector  structure  is  suboptimal  for  slow  FH  but  it  is  believed  that  the 
performance  loss  is  small,  especially  for  detection  times  that  include  a  large 
number  of  hops. 

It  is  apparent  that  other  simplifications  and  extensions  to  these  results  are 
possible.  For  instance,  it  is  assumed  that  the  starting  time  and  duration  of  the 
chip  epoch  are  known.  This  first  restriction  might  be  relaxed  by  redefining  the 
SELF  to  perform  sliding  window  integration  instead  of  the  integrate-and-dump 
operation  now  performed.  This,  of  course,  would  degrade  the  detector’s  per¬ 
formance  for  some  values  of  epoch  starting  time,  but  probably  would  exhibit 
a  better  average  performance.  There  are  also  possible  simplifications  to  the 
SELF  to  improve  its  implementability;  for  example,  removal  of  the  emphasiz¬ 
ing  function  would  make  the  detector  structure  suboptimal  but  probably  still 
asymptotically  optimal  for  small  assumed  SNRs.  Another  simplification  could 
be  coarse  subband  preselection,  where  the  total  spread-spectrum  bandwidth  is 
subdivided  into  subbands,  each  containing  a  large  number  of  chip  frequencies. 
An  algorithm  then  could  be  used  to  select  a  subset  of  the  subbands  that  most 
likely  would  contain  the  intercepted  signal.  Detailed  processing  on  these  pres¬ 
elected  bands  then  could  be  done  with  the  methods  described  in  this  chapter. 


2. A  DERIVATION  OF  SELF 


Proceeding  from  Appendix  B  of  Woodring,3  the  likelihood  function,  given 
the  carrier  phase  9  and  the  channel  Jfc,  the  conditional  likelihood  function  for 
the  i th  epoch  is 

A,  (y/M)  =  e"^e^  )T*  (2.83) 

where  E  is  the  single-epoch  energy  of  the  FH  signal,  i.e., 

E  =  25  sin^wt,!  +  9)dt.  (2.84) 

JiT„ 
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But 


E « ST/,  for  WkxTh  ^  1  (2.85) 

which,  upon  substitution  into  the  conditional  likelihood  function  (2.83)  and 
expanding  z :,(<),  yields 


A,  (y/M)  =  e-7e^(/>*5inS+Q‘C06#) 


(2.86) 


where 


( i  -t- 1  )7\ 


Pk  = 


Qk  = 


/■(<  +  ! 

-At* 


VNoTh  J,T * 

2  /‘(•+1)7\ 

iV07’h  y,r. 


y(<) cos w*.!  eft 


(2.87) 


y/No  a  /> 

Taking  expectations  with  respect  to  0  defines 


y(t)  sinw^l  dt. 


A,(y/k)  =  fe  [A, (y/M)] 


(2.88) 


which  is  the  likelihood  function  conditional  only  on  the  channel.  This  expecta¬ 
tion  can  be  evaluated  as 

£e{\ i(y/k,0)\  = 


where  =  Arg(P^  +  jQt).  Now  by  the  periodicity  of  the  integrand 

£,  [Aj(y/M)l  =  —  [2  e^Vr;+Ql™ed0  (2.91) 

In  Jo 

=  e-^o  (2.92) 

by  the  identity 

I0(a)=  [  \a^ede  (2.93) 

ilf  Jo 

where  Io  is  the  zeroth-order  modified  Bessel  function  of  the  first  kind.  Taking 
expectations  with  respect  to  the  channel  k  yields  the  single-epoch  likelihood 
function  (SELF) 


e~7 

2x 

r 2* 

1  siiiS  +  Q*  c<x»)^q 

Jo 

(2.89) 

e~7 

2  ir  , 

f  e'/2:r\/p2+Qls[n(d+'P)  (IQ 

Jo 

(2.90) 

A,(y)  =  Sk  [A,(y/fc)j 


(2.94) 


e-7 

~K~ 


) 


(2.95) 
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2.B  ASYMPTOTIC  EXPRESSIONS  OF 
THE  ALLF  MOMENTS 


We  want  to  examine  the  behavior  of  the  ALLF  moments  when  SNR  7  di¬ 
minishes  while  the  relative  SNR  r  is  held  constant.  The  asymptotic  expressions 
derived  here  encapsulate  this  behavior.  To  derive  asymptotic  expressions  for 
the  mean  and  variance  of  the  ALLF,  consider  the  mean  and  variance  M.r  and 
Vr  of  the  single-epoch  ALLF  only.  To  this  end,  it  will  be  useful  to  derive  asymp¬ 
totic  expression  of  two  functions  of  the  channel  moments:  Ht/oq  and  <r2/a2. 
Starting  with  the  first  expression  and  substituting  (2.21)  and  (2.22)  forms 


£  _  Ip&n) 

[/o(27)-1]' 


(2.96) 


We  will  need  a  partial  power  series  expansion  of  /0(x),  i.e. , 


/o(l)  =  l  +  T  +  X(Ti  +  0{x6)- 


(2.97) 


Here  and  throughout  this  discussion,  the  quantity  0(xn)  represents  any  func¬ 
tion,  say  f(x),  such  that 

x~n  lim  f(x)  <  00.  (2.98) 

z-~0 

The  particular  function  represented  by  0(xn)  is  determined  from  the  context 
of  the  equation  in  which  it  appears.  W'ith  the  above  power  series  for  / 0,  (2.96) 
becomes 


[l  +  r2  7*  +0(74)] 

72  +  b4+<?(76) 

(2.99) 

1  +2rV  +0(74) 

72  +  h4  +  ^(76) 

(2.100) 

7'2  +  (2r2  -  i)  +  C7(72) 

(2.101) 

[7-1  +  (r2  —  g)  7  +  0(j3)]2 

(2.102) 

r‘  +  ('-2-|)7  +  ^(73). 

(2.103) 

Now  let  us  evaluate  the  second  channel  moment  function.  Using  (2.21)  and 
(2.97)  plus  the  power  series  for  er,  after  carrying  out  the  multiplications  we 
get 


=  [/o(2T)  -  ir 


-f 

*  JO 


e-2^°e*/o(4rTsin|)d^-/02(2r7) 


(2.104) 
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4  =  [-r!+; i  /' 

1  —  2-,  cos <j)  +  27 2  cos2  <j>  —  57  cos3  <£3  +  |  cos4  <py4  +  <9(75)  J  x 
1  +  4r272  sin2  |  +  4r474  sin4  |  +  0(76) j 

-[l  +  rV+irV  +  0(76)]2| 

“f  =  [t2  +  ?74  +  ^(76)]  *  <  -  /  [  1  —  2  cos <?S  7 

ffo  I  *  Jo  [ 

+  ^4r2sin2  |  +  2  cos2  <pj  72  -  ^8r2co s^sin2  f  +  5  cos3  <pj  y3 
+  cos4  4>  +  8r2  cos2  <p  sin2  |  +  4r4  sin4  f  )  74  +  0(<f>,  ys)  dip 


(2.105) 


(2.106) 


-[l+2r272  +  frV  +  0(76)] 

Let  75)  be  the  particular  function  represented  by  the  symbol  0(<f>,  yh)  under 
the  integral,  then 

limj-5/  75)  =  /  \im  y~5f(<p,  y5)  <  00  (2.107) 

7-0  Jo  Jo  y—*0 

implying  that  J*  f(<j> ,  7s)  G  0( 7s).  The  interchange  of  the  limit  and  integration 
is  justified  as  follows.  The  function  f(4>,y5)  inherits  continuity  on  the  compact 
set  {<p,y  :  <p  €  [0,x]  and  7  G  (0, 1]}  from  the  integrand.  Therefore,  y~5 f(4>,y5), 
which  has  a  finite  limit  at  the  origin,  is  also  continuous  on  this  compact  set 
and  hence  is  bounded,  say  by  B,  on  this  set.  The  function  B  is  integrable 
on  <t>  G  [0,x],  from  which  the  interchange  follows  by  the  Lebesgue  Dominated 
Convergence  theorem.  The  interchange  implies  that  (2.106)  can  be  integrated 
term  wise  to  yield 

$=[7  2+h4  +  ^(76)]_1x 

{ [1  +  (2r2  +  1)72  +  2r273  +  (±  +  2r2  +  §r4)74  +  0{yb)\  (2.108) 

-  [l  +  2rV+  §rV  +  0(76)]} 

which  simplifies  to 

2 

^  =  1  +  2r27  +  2r272  +  0(y3).  (2.109) 

<ro 

With  these  asymptotic  expressions  for  /ir/<70  and  u-2/^,  we  proceed  with  the 
derivation  of  asymptotic  expressions  for  the  ALLF  moments.  The  ALLF  mean 
is  expressed  in  terms  of  the  SELF  moments  as 

Mr  =  \  ln  V  +  PVlV°}~1  {lMr  -  -  A/']2  Vo  +  [V"l  -  V<>]  Vr) 

(2.110) 
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The  last  three  terms  can  be  evaluated  as  follows: 

[2 Vi  Vo]-1  {[Mr  -  Mo]2  Vi  -  [Mr  -  M,]2  V0  +  [V!  -  V0]  Vr} 


=  te-sn— a 


-K  +  hr-i  +  2 

0q  <t0.  I  er, 


M-'}} 


(2.111) 


=  ^  [K  +  27  +  0(l2)]_1  |  [n  +  0{ 73)]2  [K  +  0(7)] 


-K[(r>-l)7  +  0(y3)Y 


+  [K  +  2r27  +  0(72)]  [27  +  272  +  0(73)]  |  (2.112) 


=  i  +  (^r!  +  IJfr)  ^  +  0iy3)- 


(2.113) 


Now  for  the  final  term.  We  will  need  the  following  power  series  expansion  for 
ln(x) 


ln(  1  +  x)  =  x  -  —  +  0(x3) 


(2.114) 


1.  (Vo 


1 ,  f,  ].! 

'o  ,n  1  “  T  +  T-! 
2  k  k  crJ 


=  -H'4 


7  +  x 72  +  0(73) 


combining  (2.113)  and  (2.117) 


Now  to  proceed  with  the  variance 


(2.115) 


(2.116) 


=  (2U7) 


Mr  =  (r2  -  5)  72  +  0(73).  (2.118) 
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The  first  term  can  be  evaluated  as 


K2  (j_  _  n  _  JL J  J  goJ 

2  \Vo  vj  2A2 


K  -  1  +  -4 
<*l 


(2.120) 


1  /  f2r-fO(T2)][A'  +  0(7)]l  (2  121) 


2  A2j  [A  +  0(7)]  J 


2rl 

A'2 


+  0(73). 


The  second  term  of  (2.119) 


(2.122) 


A-1+-* 


K 2 


-21 2 

A  -  1  + 

*oJ 


( [<r^  J  [  ffo  <^0  J  L^o  ^o'TqJJ 


(2.123) 


JL  [A  +  0(7J]  r  [2  2  2  +  Q(  3)j  [AV1  +  0(7)]  (2.124) 

A2[A+0(7)]2 

+  A[7_1  +  g7  +  ^(T3)] 

-  A  [7-1  -  h  +  0(73)]  [1  +  27  +  272  +  Cl(73)]  | 


J_  [A ±  0(7)]  (  r2A  +  2A7  +  o(72)]  +  A  [-2  -  7  +  (?(72)]  I'  (2-125) 
A2  [A  +  0(7)]2  *■  J 

=  £+0(73).  (2-126) 

K 

Combining  (2.122)  and  (2.126)  yields 

yr  =  £±V  +  0(73)  (2.127) 
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Summarizing  these  results 

m.  =  ^2  (^-i)  7*+ oft*) 

V,  =  —y’  +  Otf). 


(2.128) 

(2.129) 
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CHAPTER  3 

OPTIMAL  INTERCEPTION  VIA  A 
COMPRESSIVE  RECEIVER 


3.1  INTRODUCTION 


The  goal  of  the  intercept  receiver  is  to  detect  deceptive  electromagnetic 
sources  and  follow  up  by  extracting  waveform  features  for  use  in  the  jamming 
or  exploitation  of  that  source.  With  the  advent  of  frequency-hopped  and  other 
spread-spectrum  signals,  the  search  bandwidth  that  assures  a  reasonable  proba¬ 
bility  of  intercept  is  increased  greatly  with  a  corresponding  increase  in  the  com¬ 
plexity  of  the  intercept  problem.  For  these  cases,  wideband  interceptors  prove 
unacceptable,  while  high-performing  channelized  interceptors  prove  virtually 
unimplementable.  The  compressive  receiver,  which  simultaneously  estimates 
frequency  components  over  a  wide,  set  band,  has  promise  as  an  interceptor 
with  both  the  simplicity  of  a  wideband  device  and  the  performance  of  the  chan¬ 
nelized  device. 

The  lh  rature  is  rich1,21617'7  with  intercept  methods  for  frequency-hopped 
waveforr  There  are  also  some  analyses  of  the  detection  performance  of  the 
compressive  receiver.18,19  However,  very  little  has  been  written  on  the  appli¬ 
cation  of  the  compressive  receiver  to  the  interception  of  spread-spectrum  sig¬ 
nals  (Milstein  and  Li20  is  an  exception)  and  even  less  on  the  interception  of 
frequency-hopped  waveforms.  To  fill  the  void,  this  chapter  fully  develops  an 
optimal  and  a  simplified  subopt imal  method  for  the  detection  of  frequency- 
hopped  waveforms.  The  chapter  further  exhibits  a  structure  for  hop  frequency 
estimation. 

The  chapter  models  the  compressive-receiver  input  as  consisting  of  either 
stationary  Gaussian  noise  of  known  autocorrelation  or  noise  plus  a  hopped 
signal  of  known  hop  epoch,  unknown  phase,  and  energy  above  a  minimum  de¬ 
tectable  level.  Approximate  transfer  relationships  for  signal  and  noise  are  de- 

16A.  Polydoros,  J.K.  Holmes  and  K.T.  Woo,  “Advanced  LPI  Intercept  Detector  Research”, 
Technical  Report  No.  R8511-3,  Axiomatix,  Los  Angeles,  Calif.,  November  13,  1985. 
,7W.E.  Snelling  and  E.A.  Geraniotis,  “Sequential  Detection  of  Unknown,  Fast  Frequency- 
Hopped  Waveforms”,  IEEE  J.  Select.  Areas  Commun.,  vol.  7,  no.  4,  pp.  602-617,  May 
1989. 

,8K.D.  Breuer,  J.J  Whelehan,  and  K.  Ross,  “Compressive  Receivers  Applied  to  ESM  System 
Design”,  MSN  Microwave  Syst.  News  &  Commun.  Techno.,  vol.  16,  no.  11,  pp.  66-68, 
70,  72,  74-75,  Oct.  1986. 

19B.K.  Harms  and  D.R.  Hummels.  “Calculation  of  Detection  Probability  for  Frequency  Com¬ 
pressive  Receivers”,  IEEE  Trans.  Aerospace  and  Electronic  Systems,  vol.  AES-21,  pp. 
106-116,  January  1985. 

20K.H.  Li  and  L.B.  MiLstein,  “On  the  Use  of  a  Compressive  Receiver  for  Signal  Detection”, 
to  be  published  in  IEEE  Trans.  Commun. 
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veloped  separately  and  used  to  translate  the  detection  problem  to  an  equivalent 
one  on  the  compressive-receiver  output.  Likelihood  function  theory  is  applied 
to  the  equivalent  problem  and  yields  a  locally  optimal  (i.e.  optimal  for  small 
signal-to-noise  ratio)  detector.  The  locally  optimal  detector  has  an  unwieldy 
structure  that  defeats  the  motivation  of  using  a  compressive  receiver:  that  of 
simplicity  plus  high  performance.  Therefore,  a  time-multiplexed  detector  is  de¬ 
veloped  that,  at  the  expense  of  duty  cycle,  can  achieve  performance  as  close  to 
optimal  as  desired.  Asymptotic  statistics  of  the  detector’s  output  are  derived 
and  used  to  quantify  performance.  A  hop  frequency  estimator  is  presented  and 
its  probability  of  track  estimated. 


3.2  PRELIMINARIES 


Several  areas  need  elucidation  before  a  useful  interceptor  can  be  developed, 
in  particular,  a  precise  statistical  model  of  the  compressive-receiver  input  and  a 
precise  model  of  the  compressive  receiver  itself.  In  order  to  statistically  model 
the  compressive-receiver  output,  transfer  relationships  are  needed  for  both  the 
noise  and  signal. 

3.2.1  Input  Signal  Model 


The  signal  model  is  for  a  composite  hypothesis  problem.  Specifically,  given 
the  observation  y(t),  the  problem  is  one  of  choosing  between  // o,  which  is  the 
hypothesis  that  an  FH  waveform  is  not  present,  and  Hy,  which  is  the  hypothesis 
that  an  FH  waveform  is  present  with  an  SNR  7'  greater  than  some  minimum 
SNR  7.  Exactly,  the  model  is 


Ho  :  y{t)  =  n(t) 

versus 

Hy  :  y(t)  =  s(t)  +  n{t) 


7  <  i 


(3.1) 


where  the  frequency-hopped  signal  s(t)  is  given  by  (1.1)  and  n(l)  is  station¬ 
ary,  colored  Gaussian  noise  with  variance  <rf  and  with  autocorrelation  function 
crfRi(t).  The  hypothesized  SNR  7'  is  related  to  the  other  model  parameters  by 
7'  =  S'Th/Vi ,  while  similarly  the  minimum  SNR  7  =  STh/fff. 

Significantly,  the  signal  model  allows  for  colored  noise  and  is,  therefore, 
quite  general.  Note  that  the  model  assumes  that  all  signal  parameters  except 
amplitude  and  hop  frequency  are  known. 


3.2.2  Receiver  Model 


Figure  3.1  blocks  out  the  compressive-receiver  model.  The  compressive  re¬ 
ceiver  mixes  the  input  signal  {/;(<)  with  a  linearly  frequency-modulated  signal 

«(/)  =  cos(w0<  -  0t2)  0  <t<T,  (3.2) 

that  scans  downward  in  frequency  from  uo  to  wo  —  2 0T, .  Here  T,  is  the  scan 
time.  The  scanned  waveform  is  input  to  a  pulse  compression  filter,  hence,  the 
name  compressive  receiver.  The  filter  has  impulse  response 

h(t)  =  cos(u i0t  +  0t2)  w(t)  0  <  t  <  Tc  (3.3) 
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yi(t) 


y0(t) 


a(t) 


Figure  3.1:  Block  Diagram  of  Compressive  Receiver  Model 

where  w(t)  is  a  weighting  function  used  to  minimize  energy  spillover  between 
signals  of  different  frequencies.  The  output  of  the  compressive  receiver  now 
follows  as 

y0(t)  =  a(t  -  r)ydt  -  r)h(r)  dr  Te<t<T,  (3.4) 

=  x0(t)  +  n0(t)  (3.5) 

where 

x0(t)  i  / ^  a(t  -  r)Xi(t  -  r)h(r)  dr  Tc<t<Ta  (3.6) 

Jo 

n0(t )  ^  [TC  a(t  -  T)m(t  -  r)h(r)  dr  Tc<t<T,.  (3.7) 

Jo 

3.2.3  Output  due  to  Signal 

Using  (3.2),  (3.3),  and  the  commuted  version  of  (3.7),  the  output  of  the 
compressive  receiver  can  be  expressed  as 

x0(t)  =  f  x(r)  cos (ui0r-  0r7)  cos  [wo(<  -  r)  +  /?(<  -  t)2]  w(t  —  r)  dr  (3.8) 

Jt-Tc 

whenever  Tc  <t  <T,.  Trigonometric  manipulation  leads  to 


J  t  —  Tc 


X0(t)  -  -cos  (uiot  +  0t2)j  x(r)  cos(20tr)w(t  —  r)  dr 


+  ^sin(w0 1  +  Pt7)  f  x(r)s\n(20tT)w(t  -  r)  dr  +  c.  (3.9) 
2  Jt-Tc 


Application  of  Lemma  2  shows  that  error  term 

PwPr 


Id  < 


wo 


(3.10) 
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where  Pw  is  the  positive  variation  of  the  window  w(t)  on  t  6  [0,TC]  and  where 
Px  is  the  positive  variation  of  the  input  x(t)  on  t  £  [0,7)j.  The  definition  of 
positive  variation  appears  as  Definition  1  in  Appendix  3.B. 

The  error  bound  has  special  meaning  when  X(t)  is  a  sine  wave  of  angular 
frequency  w.  In  this  case,  Px  ss  w TJir  <  2pTaTc/ir  and  hence 

\(\  <  2PuJll~  (3.11) 

nuj0 

which  is  very  small  for  typical  values  of  u>0,  Tc ,  Ta ,  and  0. 

3.2.4  Output  due  to  Noise 


As  shown  in  Appendix  3. A  (3. 95  and  3.99),  the  normalized  autocorrelation 
(divided  by  <rj)  of  the  compressive-receiver  output  is 


=^J^Ri(ul-d)  | 


,2T 

Jim) 


2TC-|U, 


cos[(u»o  -  2 pt  +  /?U2)(ui  -  d)]  cos  [(wo  +  /?m2)«i] 


W  '  I"1  )  w  (U~  2  ~‘)  du 2  }  dui  +  c' 

The  error  term  e  is  bounded  as 


(3.12) 


I  t\<~PlBTe 


(3.13) 


where 

B_  1  1  2  2 
“  2pTc  +  2w0  -  2 0Tt  +  w0  -  2 pT,  +  -  pTc 


(3.14) 


with  Pu,  being  the  positive  variation  of  w(t)  defined  by  Definition  1  in  Appendix 
3.B. 

Under  typical  operating  constraints,  the  error  bound  can  be  simplified  fur¬ 
ther.  The  term  2P(T,  —  Tc)  represents  the  total  frequency  spanned  by  the 
compressive  filter,  which  is  very  large  (typically  on  the  order  of  megahertz). 
Additionally,  the  frequency  w0  is  usually  in  the  tens  to  hundreds  of  megahertz 
range,  hence  wo  106.  These  two  facts,  along  with  the  fact  that  the  scan  time 
is  typically  twice  the  compression  time  (i.e.  T,  =  2 Tc),  imply  that  B  <  l/pTc. 
Under  these  assumptions,  the  error  is  bounded  as 


M< 


i  Pi 

8  0  ' 


(3.15) 


Of  interest  are  special  cases  of  the  autocorrelation.  When  the  input  noise  is 
white,  meaning*  that  crfRi(t)  =  6(t)No/ 2,  then  the  output  noise  is  stationary 

‘Since  the  variance  of  a  white  noise  process  is  undefined,  arbitrarily  let  oj2  =  No  where 
No  is  the  single-sided  spectral  density  of  the  white  noise  process.  This  choice  make  the 
signal-to-noise  ratio,  "r,  consistent  with  other  definitions  in  the  literature. 
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and  has  autocorrelation 

j  r?Tc-\d\ 


RM  =  cos  [\d\(u0  +  0u2)]  W  (j2-+  — ^  W  du2 

(3.16) 

whenever  |</|  <  Tc,  otherwise  R0(d )  =  0.  If  the  window  function  w  is  rectangu¬ 
lar,  then 

_  1  sin  [0\d\(Tc  -  |d|)] 


*‘(d)=8  m 


■  cos  [|dj(wo  +  0Tc)\ 


(3.17) 


whenever  |<f|  <  Tc,  otherwise  R0(d)  =  0.  Regressing  to  the  case  of  general 
stationary  noise  but  now  considering  only  rectangular  windows, 


where 


Ro(t,d)  =  i  j  «,(«!  -  d)g(uut,d)  dux 

sin  [(0d  -  20ui)(Tc  -  |m|)] 


(3.18) 


g(ux,t,d,)  = 


(0d  —  20u\) 

x  cos  [w0d  -f  2 0t(u\  -  d)  —  2woU!  +  (0d  -  20u\)Tc ] 


+  Vn.l?d(T'  Mn  CQS  [uod  +  20t{ui  _d)  +  0dTc]  (3.19) 

0d 


3.3  LOCALLY  OPTIMAL  DETECTOR 


We  aim  to  develop  a  locally  optimal  detector  of  frequency-hopped  wave¬ 
forms  based  on  a  compressive-receiver  output.  Keying  on  the  fact  that  the 
optimal  detector  of  frequency-hopped  waveforms  integrates  coherently  over  a 
single  hop  period,17  we  conjecture  that  an  optimally  configured  compressive 
receiver  should  integrate  over  a  period  commensurate  with  the  hop  epoch  7),. 
But  because  we  are  also  interested  in  the  detector’s  performance  in  estimating 
the  hop  frequency,  we  want  to  eliminate  interference  from  adjacent  hops.  We 
thus  choose  T,  =  T),  and  assume  that  the  compressive  receiver  is  synchronized 
to  frequency  hops.  This  is  not  a  realistic  assumption  in  the  pure  detection 
problem  but  it  will  lead  to  an  optimal  detector  whose  performance  degrades 
gracefully  upon  relaxing  this  assumption. 

Because  the  interfering  noise  is  typically  of  much  larger  bandwidth  than  the 
hop  rate,  the  correlation  between  hops  is  negligible  and  so  the  optimal  multi¬ 
hop  detection  statistic  is  simply  the  sum  of  the  optimal  single-hop  detection 
statistics.  We  thus  confine  ourselves  to  the  problem  of  using  the  compressive 
receiver  to  optimally  detect,  given  an  observation  period  of  T,,  a  sine  wave  of 
unknown  amplitude  and  phase  and  whose  frequency  is  one  of  the  known  hop 
frequencies. 

Based  on  the  above  assumptions,  the  detection  problem  is  now 


Xi(t )  =  Tii(t) 

versus 

Hy  :  X{(t)  =  V2Ssin(uict  +  0)  +  n,(<) 


(3.20) 
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for  7  <  y'  and  Tc  <  t  <  T,.  The  parameters  9,  7,  7',  uik,  and  n,(t )  are  as 
defined  in  Section  3.2.1. 

By  using  the  results  in  Sections  3.2.3  and  3.2.4,  the  detection  problem  based 
on  the  output  of  the  compressive  receiver  is 

Ho  :  x0(t)  =  n0(t) 

vers  us 

Hyi  :  x0(t)  =  '/2S' cos9  yc(t,uk)  +  VliS' sinO  ya(t,tok)  +  n0(t) 


(3.21) 

for  7  <  7'  and  Tc  <  t  <  T,  and  where  n0(t)  is  stationary,  colored  Gaussian 
noise  with  autocorrelation  function  R0(t)  as  defined  by  (3.12)  and 

yc(t,Uk)  =  ^cos(w0<  +  /?<2)  f  cos(ukr)cos(2f3tr)w(t  -  r)  dr 
2  Jt—Tc 


1/ 

Jt-Tc 


+  i  sin(w0<  +  fit2)  I  cos(ojkr)s'\n(20tr)w(t  -  r)  dr 


(3.22) 


a  1  [* 

l/s(t,uk)  =  -cos  (*j0t  +  /3t7)  I  sin(ukr)  cos(2/3tr)w(t  -  r)  dr 

2  Jt-Tc 

1  f* 

-P  -  sin(u>o<  -f  0t2)  I  sin(uikr)sin(20tr)w(t  —  r)  dr. 

2  Jt-Tc 

(3.23) 

From  Van  Trees21,  the  conditional  log-likelihood  function  for  this  problem  be¬ 
comes 

\nA[x0(t)/uit,9,y]  = 

pk™9  /;  x0{r)gc(r,ojk)  dr  +  R^Sin#/r.’  *o(T)g, (r,uk)  dr- 
7  [T‘ 

2rJr  [cos0  ye(r,wk)  +  s\n0  y,(r,u>k)]lcos9  gc(r,wk)  +s\n0  g,{r,uik)]  dr 

(3.24) 

where  the  functions  gc(t,u/k)  and  g,(t,u)k)  are  respectively  defined  by  the  inte¬ 
gral  equations 

rT. 


fT • 

/  Ro 
Jtc 


T  +  t  t 

— 'r-‘j 

9c{r,ujk)  dr  = 

Vc(t,Vk) 

(3.25) 

r+t  . 

-  r  —  t 

2 

g,(r, u>k)  dr  = 

(3.26) 

J1H.L.  van  Tree*,  Detection,  Eitimation ,  and  Modulation  Theory,  Part  /,  Wiley,  New  York, 
1968. 
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for  Te  <t  <T,.  Since  we  are  interested  in  a  locally  optimal  test  (i.e.  small  7), 
we  neglect  the  last  term  of  (3.24)  and  say 

In  A  [xo(t)/u>k,0, 7]  * 


Averaging  this  approximate  likelihood  ratio  over  8  and  w*  yields 


AM0/7]  =  f> 


k=i 


(3.27) 


(3.28) 


where  /0  is  the  modified  Bessel  function  of  the  first  kind  and  zero  order  and 
where  the  complex-valued  function  G*(<)  is  defined  as 


G*(<)  =&(«,«*) +*$»(«,«*)•  (3.29) 

Consider  again  the  small  7  case  and  note  that  Io(x)  ss  1  +  x2/4  for  small  x. 
Conjure  (3.28)  into  a  locally  optimal  statistic 


r  = 


2 


*o(r)G*(r)  dr 


(3.30) 


where  the  scale  factor  l/a'-Th  is  added  for  convenience  in  future  analyses.  Fig¬ 
ure  3  2  blocks  out  (3.30).  To  complete  the  detector,  T  is  compared  against  a 
threshold  f,  whose  value  determines  the  probability  of  false  alarm  Pp.  (Sec¬ 
tion  3.5  shows  the  exact  relationship  between  t>  and  Pp.)  The  statistic  7  being 
locally  optimal  will,  for  small  signal-to-noise  ratios,  yield  the  greatest  possible 
probability  of  detection,  hence  it  is  locally  the  most  powerful  test  in  terms  of 
signal-to-noise  ratio. 


3.4  TIME-MULTIPLEXED  DETECTOR 


The  locally  optimal  detector  of  the  last  section  efficiently  detects  frequency¬ 
hopping  waveforms.  As  will  be  shown,  it  rivals  the  optimal  detector  that  directly 
observes  the  original  time  waveform.  Unfortunately,  it  also  rivals  the  optimal 
detector  in  implementation  complexity  and  thus  undermines  the  attractive  sim¬ 
plicity  of  the  compressive  receiver.  In  this  section,  we  construct  a  detector  based 
on  the  locally  optimal  detector  that  maintains  simplicity  of  implementation  for 
a  small  performance  cost. 

The  time-multiplexed  detector  depicted  in  Figure  3.3  is  both  efficient  and 
simple.  The  detector  consists  of  a  complex  filter  whose  impulse  response  G^(t) 
is  constructed  from  the  pseudo-signals  G*(t)  by  the  equation 

K 

Gs(0  =  £CjP;-t  +  (j-l)Ar].  (3.31) 

j  =  i 
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Complex  Envelope-Squared 

Matched  Filter  Detector 


Figure  3.2:  Locally  Optimal  Single-Epoch  Detector 


Complex  Envelope-Squared  Sampler/ 

Matched  Filter  Detector  Accumulator 


Figure  3  3:  Time-Multiplexed  Detector 
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The  envelope  of  the  filter’s  squared  output  is  then  sampled  at  times, 

t  =  T,  +  (k  -  1)AT  for  k-  1,  -,K  (3.32) 

then  summed  and  scaled  by  1/afTh  to  produce  the  test  statistic  I\  It  is  easily 
shown  that 


f 


2 


x0(t)Gx  [Tj  +  (fc  -  1)AT  -  r]  dr 


(3.33) 


which  by  defining 

Gk(t)  =  Gx[Ta-t  +  (k-  1)A T)]  for  Tc<t<T,  (3.34) 


yields  the  alternate  expression 


f  =  Jv  E  i  r  *<>(*)&(*)  dr 

Note  that  (3.31)  transforms  (3.34)  to 

K 

Gk(t)  =  YlGj[(j-k)AT  +  t) 

7  =  1 


(3.35) 


(3.36) 


which,  for  AT  >TS-  Tc,  implies  Gjt(<)  =  G’t(I),  since  the  GjS  are  zero  outside 
the  range  Te  <  i  <TS.  Hence  the  time-multiplexed  detector  is  equivalent  to 
locally  optimal  detector  for  this  choice  of  AT.  We  aim  to  show  that  AT  can 
be  made  significantly  smaller  for  a  small  performance  cost,  thus  concluding 
that  the  time-multiplexed  detector  captures  most  of  the  detectability  of  the 
locally  optimal  detector  but  maintains  the  simplicity  of  the  compressive-receiver 
cunfiguiation. 


3.5  DETECTOR  PERFORMANCE 
ANALYSIS 


In  both  the  locally  optimal  detector  and  the  time-multiplexed  detector,  the 
test  statistic  is  the  sum  of  squares  of  a  large  number  of  weakly  correlated  random 


variables.  Namely,  for  the  locally  optimal  detector, 

2/e 

II 

(3.37) 

where 

Ci  =  nrr  1  xo(r)hj(T)dT 

VffiTh  Jtc 

(3.38) 

and  where  hj(t)  is  defined  as  below  with  1  <  m  <  K, 

(3.39) 

m(<)  =  9j(<.Wm). 

(3.40) 
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Similarity  for  the  time-multiplexed  detector,  the  test  is 


2  K 


F  =  £</ 

(3.41) 

>=i 

where 

Cj  =  /-f—  [  x0(r)hj(T)dT 

V<*iTh  Jtc 

(3.42) 

and  for  1  <  m  <  K, 

^2m- 1(0  —  5/c(Lwm) 

(3.43) 

^2m(0  =  ga(t,Um). 

(3.44) 

Here,  gc(t,um)  =  3?  [5m(t)J  and  g,(t,ujm)  =  0  Gm(t) ]. 

In  the  analysis  to  fol- 

low,  the  hat  notation  will  be  dropped  since  the  results  apply  to  each  detector 
in  the  same  way.  In  other  words,  to  get  the  result  for  the  time-multiplexed 
detector,  add  hats  to  the  appropriate  variables. 


Because  the  test  statistic  is  the  sum  of  a  large  number  of  weakly  corre¬ 
lated  random  variables,  there  is  reason  to  believe,  despite  the  correlation,  that 
the  statistic  has  approximately  Gaussian  distribution.  We  proceed  under  this 
assumption  with  justification  to  follow  later.  To  specify  the  asymptotic  dis¬ 
tribution  of  T,  we  need  its  mean  and  variance  under  the  signal-present  and 
signal-absent  hypotheses.  For  this  purpose,  define  1  <  j  <  2A',  as 

Z2m-l(<)  =  yc(t,Wm)  (3.45) 

Z2  m(t)  =  y,(t,Wm)  (3.46) 

with  m  ranging  between  1  and  K,  while  2j(t),  for  1  <  j  <  2K ,  is  defined  as 

gc(t,um)dT  (3.47) 

T-r  05(<,wm)dr  (3.48) 

with  m  also  ranging  between  1  and  K .  For  the  same  reason  as  above,  we  define 
the  time  cross  correlations 

£m,n  =  JT  zm(T)hn{r)  dT  (3.49) 

(m,n  =  /  2m(r)hr(T)  dr  (3.50) 

7/1  Jtc 

for  1  <  m,n  <  K.  (For  the  locally  optimal  detector  case,  note  that  (j  *  =  £;,*•) 
Assume  now  that  the  signal  is  at  frequency  w;.  Then 

<*  =  Y/2ycos0£2,_,,t  +  sin  6  £2/,t  +  rjk  (3.51) 


hm(t) 


=  I  '  Ro 
Jr* 


Z2m-l(<) 


=  [  ‘  Ro 

Jtc 
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where  the  random  variable  9  is  uniformly  distributed  on  [0,2tr]  and  the  r^s  are 
zero-mean  Gaussian  with  covariances  =  £;  t-  From  (3.132)  in  Appendix 

3.C 


Pk/I  -  £{(,l/ult)  =  YtfZi-i  t  +  (21, k)  +  £*,*• 
When  averaged  over  /  and  summed  over  k,  the  mean  of  T  is 

,  2  K  2 K  2  K 

My  =  ^EE^'  +  XA*- 

1=11=1  i=l 


(3.52) 


(3.53) 


Use  (3.51)  and  (3.136)  to  construct  the  covariance  between  the  jth  and  Fth 
terms  of  T  when  the  signal  is  at  frequency  w/.  The  result  is 

PlA/l  =  2(t')2  61-1  j  £21.2  £2 l-l,k  £2 l,k  +  ~~  £?/_ij  £21-1,* 


,  il')  c2  <-2  (7')2  c2  c2  (V)2  r2  r2 

f  9  <2/,j  <2 l,k - 7)  <21-1,2  <21  ,k  7  <21, j  <2/-l,i 


+  47^21-1  ,j(.2l-\,kij,k  +  W^2l,j(2l,kij,k  +2  £?  k 

which  upon  averaging  over  /  and  summing  over  j  and  it,  becomes 
i-,'\2  2 ft  2 K  K  , 

Vi‘  =  (  2  ^2/- !  >  £21,2  £21-1,1  £21,* 


(3.54) 


2=1 *= 1 1=1 


£2/  —  lj  £21, k  -  2  ^21,2 


,  2K  2K  2 K  ,  ,  \  2 K  2K 

+  1  6  .* )  + 2  E  E  <?.»• 

2=1 *=1  /=i v  7  2=1 *=1 


(3.55) 

Of  course,  for  the  signal- not- present  case,  the  mean  and  variance  are  simply 
(3.53)  and  (3.55)  with  the  signal-to  -noise  ratio  7'  =  0. 

Since  the  test  statistic  T  has  an  approximately  Gaussian  distribution,  the 
threshold  v  and  probability  of  detection  Pd,  for  a  given  probability  of  false 
alarm,  follow  as 

v  =  \A7>-'(  1  -  Pr)  +  Mo  (3.56) 


Py  =  1  -  <t> 


>Ao*-,(1-/V)-Wv  +  A/0 


(3-57) 


where  ^(x)  is  the  distribution  function  of  the  standard  Gaussian. 

We  now  justify  the  above  use  of  the  central  limit  theorem.  First  let, 


(3.58) 


£2/1,2  K 
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which  is  the  covariance  matrix  of  the  Qs,  and 

&>/— 1,1  cos0  +  6/,i  sin0 


E$i 


{2/-1.2A' COS0  +  f2t,2K  Sin0  J 


(3.59) 


which  are  l/v/^T7  times  the  means  of  the  <y s  under  the  condition  that  the  signal 

has  phase  0  and  frequency  w/.  We  note  that,  since  £  is  nonnegative  definite 

1  XI 

and  symmetric,  there  exists  a  square-root  matrix  £2  such  that  £  2  £2  =  £. 
Consider  also  the  diagonalization  of  £~  2  £~  £  2  =  T  AT  where 


A  = 


A,  0 
0  A2 

0  0 


0 

0 

A2A 


(3.60) 


is  the  matrix  of  eigenvalues  of  £  and  T  is  an  orthogonal  matrix  of  eigenvectors. 
Use  the  above  diagonalization  and  £ 2  to  rewrite  the  test  statistic  as 


r  =  (g  +  v/2 YMB],y  A  (g  +  v/2 VM#,,) 


(3.61) 


where 


(3.62) 


with  {^}  independent,  zero  mean,  unity  variance,  and  Gaussian  and  where 


M».i  = 


=  Tt£-2E«,. 


(3.63) 


The  test  statistic  T  is  now  the  sum  of  squares  of  independent  Gaussian  variables. 
Through  application  of  the  Berry-Esseen  Theorem,22  T,  conditioned  on  0  and  l, 
is  approximately  Gaussian  distributed  with  an  error  no  more  than  4 c/cr  where 


,  487'm2#  (  +  7 

8y'm?e j  +  2 

(3.64) 

2  K 

+  2)‘ 
i—  1 

(3.65) 

If  this  error  bound  is  small,  a  fact  that  must  be  established  numerically,  then 
the  CLT  applies  uniformly  to  the  conditional  distribution  of  P.  If,  in  addition, 
the  overall  mean  A/y  and  variance  Vy  remain  essentially  constant  with  respect 

32  A.  Papoutis,  Probability,  Random  Variablea,  and  Stochaatic  Proctaaca,  Second  Edition, 
McGraw-Hill,  New  York,  1984. 
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to  /  and  9,  then  the  CLT  applies  to  the  unconditional  distribution  of  T  as  well. 
This  fact  must  also  be  established  numerically. 

The  above  analysis  using  the  asymptotic  distribution  of  T  is  supplemented 
with  upper-  and  lower-bounding  distributions.  Specifically, 


l  -Qk 


<  Pr  [r  <  a]  <  I  -  Qk 


(3.66) 

where  Amax,  Amin  are  respectively  the  maximum  and  minimum  eigenvalues  of 
the  covariance  £  and  where  emax,  emjn  are  respectively  the  maximum  and 
minimum  over  all  eigenvalues  of  the  matrices 


Bi  =  AT  £Ai,  1<1<K 


with 


At  = 


£21-1.1  £21,1 


[  £21-  1,2/e  in, ik  J 

and  finally  where  Qm  is  the  generalized  Marcum  Q-function  defined  as 

Qm(a,0)  =  f°°x  e-*(**+°a)/m_1(a*)  dx. 


(3.67) 


(3.68) 


(3.69) 


Of  interest  is  that  the  upper  bound  equals  the  lower  bound  only  when  the 
channel  outputs  are  i.i.d.  and  the  sum  of  the  square  magnitude  of  the  signal 
component  across  the  channels  is  independent  of  signal  phase.  In  a  sense,  the 
bounds  give  an  indication  of  how  well  the  detector  fits  the  i.i.d.  assumption, 
since  the  upper  bound  corresponds  to  the  detector  distribution  under  the  i.i.d. 
and  phase  independence  assumptions,  but  with  an  increased  noise  level,  while 
the  lower  bound  has  the  same  interpretation,  but  with  a  decreased  noise  level. 
These  bounds,  when  averaged,  approximate  the  detector  distribution,  the  use¬ 
fulness  of  which  will  be  studied  and  compared  with  the  asymptotic  distribution 
in  Section  3.7. 


3.6  HOP  FREQUENCY  ESTIMATOR 


The  time-multiplexed  detector  can  be  modified  to  estimate  hop  frequency 
as  depicted  in  Figure  3.4.  Here,  instead  of  accumulating  samples  of  the  square 
envelope,  the  maximum  sample  is  found  as 

|<7.|2  =  max  |Gi|2.  (3.70) 

'  1  <1<K 

Then  u >tj  is  declared  the  hop  frequency  estimate. 

We  need  the  performance  of  this  maximum  likelihood  estimator.  For  this 
purpose,  the  probability  of  track  Pr  is  defined  as  the  probability  that,  given 
the  presence  of  a  signal,  the  hop  estimate  matches  the  actual  hop  frequency. 
Because  of  the  lack  of  symmetry  and  independence  between  the  squared  enve¬ 
lope  samples,  an  exact  expression  of  Pr  that  is  also  computable  is  very  difficult 
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X0(t) 


Complex  Envelope-Squared  Hop 

Matched  Filter  Detector  Estimation 

f - A - N  f - " - \  f - A - h 


Figure  3.4:  Hop  Frequency  Estimator 


to  find,  although  it  can  be  expressed  as  a  A'-dimensional  integral.  We  will  be 
satisfied  to  tightly  bound  Pt  when  it  is  large.  The  results  will  be  accurate 
under  the  useful  operating  conditions  of  the  estimator. 

Proceeding  in  this  fashion,  we  will  need  to  define  the  distribution  function 
of  each  sample 

Fk(z)  =  Pr  (|G*|2  <  x)  (3.71) 


and  note  that  its  density  /t(x)  exists  and  will  be  computed  exactly  later.  Now 
suppose  the  signal  is  at  hop  /;  then  the  probability  that  |(7/|2  is  between  x 
and  x  +  Ax  is  approximately  ft(x)Ax  and,  given  that  it  is  between  these  two 
numbers,  the  probability  that  it  is  the  maximum,  hence  correctly  chosen,  is 
approximately 


T,(x)  =  Pr 


fl  (lG*!2<*) 


*  =  < 
l*  *  i 


(3.72) 


Therefore,  upon  letting  Ax  — >  0  and  averaging  over  /, 


M*)T,(x) 


dx. 


(3.73) 


In  order  to  produce  the  promised  bounds  on  Pr,  we  will  need  to  produce 
bounds  on  7}(x).  Using  De  Morgan’s  law, 


T,(x) 


1  -  Pr 


U  (I^I2>-) 


\k  =  l 
Lit  1 1 


(3.74) 


=  1-  £  Pi(*)+  £  W  +  ■  •  •  (3-75) 

singles  pairs 

~  /  -I 

where  Pi(x)  is  the  probability  that  a  single  sample  exceeds  x  and,  similarly, 
P2(x)  is  the  probability  that  both  of  a  given  pair  of  samples  exceed  x.  The 
above  expression  is  useful,  because,  for  the  case  of  large  x,  the  events  that 
samples  other  than  the  /th  exceed  x  are  approximately  disjoint.  In  this  case, 
the  first  two  terms  of  the  series  accurately  approximate  7)(x).  With  enough 
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terms,  this  series  can  determine  T;(x)  to  any  desired  degree  of  accuracy  and, 
being  alternating,  implies  bounds  on  7/(z )  since  the  series  out  to  positive  terms 
is  above  the  actual  probability  and,  similarly,  the  series  out  to  negative  terms 
is  below.  We  will  use  the  first  three  terms  of  the  series  in  just  this  way  along 
with  (3.71)  to  write 

b,(x)  <  T,(x)  <  b,(x)  +  t,(x)  (3.76) 

where 

K 

6, (*)  =  1  -  Y,  [1  -  **(*)!  (3-77) 

i  =  1 

K  K 

d(x)  =  53  5Z  {1-m®*  [**(*),  fj(*)]} 

i  =  i  j  =  l 
*  t  I  j  t  l,k 

Equations  (3.73)  and  (3.76)  now  determine 

Pt  >  Jr  T"  f  fi(x)m&x[b(x),0]  dx 

K  ,=  i  Jo 

Pt  -  7^53  /  /*(*) min  [&(z)  +  f(x),  1]  dx  (3.80) 

h  i=i  Jo 

which,  along  with  explicit  expressions  for  fi(x)  and  F/(z)  below,  complete  the 
performance  analysis  of  the  hop  estimator: 


(3.78) 


(3.79) 


fiW=  '  /^  r  rTe-«z-Mf s-\z-M)Mdedr 

\E\*  Jo  Jo  Jo 


where 


and 


and 


S  = 


^2/-  1 ,2fc—  1  £2f-l,2* 

£2I,2*-1  £21,2  k 


z  = 


r  cos  <f> 
r  sin  <t> 


M  =  v/277 


£21-1,2*-!  co»0  +  £21.2*  —  1  sinf? 
£21-1,2*  cos  9  d-  £21,2*  sin# 

following  upon  differentiation  is 

r2ir  r2w 


1M  =  Jrr  /  7  c-«v-Mrr-'(V-M)  ^ 
<112/ 1 2  Jo  Jo 


where 


Y  = 


i/z  cos  <t> 
yfx  sin  <t> 


(3.81) 

(3.82) 

(3.83) 

(3.84) 

(3.85) 

(3.86) 
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3.7  PERFORMANCE  COMPARISONS 


This  section  graphically  compares  the  performance  of  the  locally  optimal 
detector,  based  on  the  compressive-receiver  output,  to  that  of  the  optimal  de¬ 
tector,  based  on  receiver  input.  Also  evaluated  is  the  performance  of  the  time- 
multiplexed  detector  as  a  function  of  the  sampling  AT. 

The  parameters  chosen  to  make  comparisons  are:  Tc  =  50  ps,  T,  =  1000 
ps,  vo  =  2tt  X  40  Mhz,  the  minimum  hop  frequency  equals  2  Mhz  and  the 
maximum  hop  frequency  equals  4  Mhz,  100  hop  frequencies,  and  a  hop  rate  of 
20  Khops/sec. 

Figure  3.5  shows  how  the  locally  optimal  detector  compares  with  the  opti¬ 
mal  detector  based  on  the  original  observations.  As  expected,  for  low  SNRs  the 
locally  optimal  and  the  optimal  compare  favorably.  On  the  other  hand,  there 
is  about  a  3  db  difference  between  the  performances  in  the  high  region.  Two 
factors  are  responsible.  One  is  that  the  analysis  of  the  optimal  receiver  used  the 
CLT,  a  poor  model  when  one  channel  dominates  as  in  the  high  SNR  case.  It 
is,  in  fact,  optimistic.  Another,  as  explained  later,  is  a  modeling  phenomenon 
that  can  account  for  up  to  3-db  error  between  predicted  and  actual  perfor¬ 
mances.  Also  note  how  the  distribution  bounds  are  pessimistic  relative  to  the 
CLT  analysis.  This  also  may  be  due  to  the  inappropriateness  of  the  CLT  when 
one  channel  dominates. 

Figure  3.6  shows  how  performance  degrades  immediately  as  AT  becomes 
less  than  the  compression  time  Tc.  This  indicates  that  any  tradeoff  between 
performance  and  duty  cycle  would  not  be  worth  the  degradation. 

Figure  3.7  differs  from  the  previous  case  in  that  the  noise  is  bandpass  instead 
of  white.  It  has  the  same  general  character  as  far  as  the  degradation  with  respect 
to  AT  but  is  approximately  3  db  better.  This  is  not  an  actual  performance 
difference  but  a  modeling  phenomenon.  As  far  as  noise  analysis  is  concerned, 
the  compressive  receiver  is  a  mixer  followed  by  a  narrow  bandpass  filter.  If  the 
input  noise  is  white,  there  will  be  uncorrelated  noise  contributions  from  the 
sum  and  difference  frequencies  produced  by  the  mixer,  hence  the  3  db.  A  more 
realistic  scenario  has  noise  of  uniform  spectral  density  across  the  analysis  band 
of  the  compressive  receiver.  In  this  case,  the  sum  contribution  is  filtered  out. 
This  reasoning  also  explains  the  difference  between  the  optimal  and  the  locally 
optimal  detectors. 
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Figure  3.5:  Performance  of  Locally  Optimal,  Compressive  Receiver- Based  De¬ 
tector  versus  Optimal  Detector  Based  on  Direct  Observations 
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Signal  to  Noise  Ratio  (db) 


Figure  3.6:  Performance  of  Locally  Optimal  Detector  Compared  with  Time 
Multiplexed  Detector,  White-Noise  Case 
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3.8  CONCLUSIONS 


Presented  were  two  detectors  of  frequency-hopped  waveforms  based  on  the 
compressive  receiver.  The  first  was  developed  by  applying  the  likelihood-ratio 
theory  to  the  observed  compressive-receiver  output  and  yielded  a  locally  op¬ 
timal  (low-SNR)  detector.  The  second,  motivated  by  the  simplicity  of  imple¬ 
mentation,  was  a  time-multiplexed  version  of  the  first  that,  through  the  choice 
of  a  parameter,  either  could,  at  the  expense  of  a  low  duty  cycle,  achieve  the 
detectability  of  the  first  or,  at  the  expense  of  degraded  performance,  achieve 
higher  duty  cycles.  The  second  detector  was  modified  into  a  maximum  likeli¬ 
hood  estimator  of  hop  frequency.  Performance  of  both  detectors  and  the  hop 
frequency  estimator  were  analyzed  and  compared. 

The  compressive  receiver  fulfills  its  promise  as  a  simple,  yet  high-performing 
interceptor.  The  performance  of  the  locally  optimal  detector  shows  that  rela¬ 
tively  little  detectability  is  lost  in  the  compressive-receiver  processing.  Most  of 
the  discrepancy  is  due  to  the  difference  in  coherent  integration  time  (one-half  for 
the  parameters  used).  Furthermore,  for  a  small  performance  cost,  the  simplicity 
of  the  compressive-receiver  approach  can  be  retained  by  the  time-multiplexed 
detector.  The  hop  frequency  estimator  again  compares  favorably  with  the  cor¬ 
responding  device  that  used  raw  input  instead  of  compressive-receiver  output. 


3. A  DERIVATION  OF 

COMPRESSIVE-RECEIVER 

AUTOCORRELATION 


The  normalized  autocorrelation  of  the  compressive-receiver  output  is  defined 
to  be 

1 


Ro(t,d)=  -*E 


under  the  restrictions  that 


n°  0 + 2)  n° 

Tc<t<T, 


d 

1  2 


(3.87) 

(3.88) 


and 

^<min  (t-Tc,T,-t).  (3.89) 

Substitute  the  expression  for  the  output  noise,  (3.7),  interchange  expectations 
and  integration,  and  use  the  definition  of  the  normalized  input  correlation; 
R,(t)  =  S[n,(t)rti(t  +  r)]/af,  to  get 


1  /Tc  fTc 

Ro(t,d)=-j  \ 

ai  J a  Jo 


R>{t\  -  r2-d) 


COS 


cos 


cos  [won  +  (dry]  w(ti  ) 


cos  [w0r2  +  0t2]  w(r2)  dr2  dry. 
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Table  3.1:  Coefficients  of  (3.94) 


m 

U)j 

m 

i 

u0t  -  2 (3l2  -  ^d2  +  /3dm 

2  fit 

0 

2 

—uiq  i  +  2(3t2  +  — d2  —  (3du\  +  (3u 2 

<5. 

l 

o 

3 

3 

-(wo  ~  20t)(«i  -  d)  +  | u2 

w'o  +  0d  —  0%i\ 

4 

—u>at  +  2  f3t2  +  ^ d 2  —  (uiq  +  pd)ui  +  ^u2 

w0  -  20t  -  0u\ 

0 

2 

5 

(u0i  -  2/3l)(m  -d)+  ^ u 2 

w0  -  (3d  +  /?«! 

0 

2 

6 

—u>ot  +  20t2  +  —  d2  +  (wo  —  (3d)u\  + 

w0  -  2/3t  +  /?«! 

0 

2 

< 

wo d  +  2(3t{ui  -  d) 

(3d 

0 

8 

(uig  —  20t)d  4-  (— 2wq  +  2/3t)ui 

fid  —  2fixi\ 

0 

(3.90) 

To  exploit  the  stationarity  of  the  input  noise,  first  transform  the  above  integral 
with 


c 

II 

J 

1 

J 

(3.91) 

u2  =  n  +  r2 

(3.92) 

and  then  reduce  the  cosine  components  with  multiple  applications  of  the  identity 
cos(A)  cos(B)  =  l/2cos(A  +  B)  +  1/2  cos(A  —  B)  (3.93) 


to  form 

R0(t,d) 


j  cos(gj  -t-u )jU2  +  (3:ul)w  ^U2-^  —  J  w  ~2~^)  du2}dui 


(3.94) 

where  0j,  Wj,  and  (3j  are  given  by  Table  3.1.  Use  more  trigonometry  and  apply 
Lemma  2  to  terms  3-8  to  rewrite  the  autocorrelation  as 
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Rc{t,d)  =  J/J  *(„-*{ 

/■2T,-|u,| 

/  cos  [(w0  -  2/Ji  +  /?u2)(ui  -  d)]  cos  [(w0  +  0u2)ui] 

J  l«.l 

W  (  ~)  w  (  du2  |  dui  +  f  (3.95) 


with  error  term 


<  1  p2  V  f  C  l^(ui  ~  d)j 

_  8  w  jLj  J_Tc  (uij  +  20 j  |ui  |) 


dui 


(3.96) 


where  is  the  positive  variation  of  w(t)  defined  by  Definition  1  of  Appendix 
3.B.  Simplify  the  error  bound  by  minimizing  each  term,  Wj  +  20j\ui\,  with 
respect  to  u i,  d,  and  t,  while  noting  the  restrictions  (3.88)  and  (3.89).  The 
result  is 

rT< 

II  1  n  • 

k  < 


|/^(«i  -d)|du, 


where 


B 


1 


+ 


1 


+ 


2 0TC  2u>o  -  20T,  w0  -  2/37)  w0  -  3TC ' 

The  relation  |/?(<)|  <  1  further  simplifies  the  bound  as 

kl  <\pIbtc. 


(3.97) 

(3.98) 

(3.99) 


3.B  BOUNDS  ON  INTEGRALS  OF 
LINEARLY 

FREQUENCY-MODULATED 

SINUSOIDS 


The  first  bound,  Lemma  1,  is  a  tool  used  only  to  prove  the  second  bound, 
Lemma  2,  which  is  used  for  derivations  concerning  the  output  of  the  compressive 
receiver:  namely,  the  derivation  of  the  noise  autocorrelation  and  the  derivation 
of  a  simplified  expression  for  the  output  signal  component. 


Lemma  1  Under  the  restrictions  that  b  >  a,  w  >  0,  0  >  0,  and  w  +  2 0a  >  0, 


cos(9  +  w<  +  0t 2)  dt 


2 

~  w  +  20a 


(3.100) 


Proof 

Define  the  function 

T (6,u,0)=  [  sin((?  +  w<  +  0l2)  dt  (3.101) 

Jo 
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where 


«.  = 


-1  +  0+4  0(x-6) 

20 

x  —  6 


u> 


0>  0 
0  =  0 


(3.102) 


6  =  6  mod  x. 


(3.103) 


This  is  simply  the  integral  of  sin(0  +  ut  +  0t 2)  from  t  =  0  to  its  first  zero 
crossing. 

Preliminarily,  three  facts  need  to  be  proven:  first,  that  T (6,u,0)  decreases 
with  respect  to  0\  second,  that  it  decreases  with  respect  to  w;  and  third,  that 
with  0  =  0  it  decreases  with  respect  to  6. 

Beginning  with  the  first  fact  we  will  show  that  T (6,u>,0)  decreases  with 
respect  to  0  by  proving  its  partial  derivative  to  be  negative.  Employ  the  chain 
rule  to  get 

(6,U,0) 


60 


=  /'' 


V  cos (6  +  cut  +  0t2)  dt. 


Let  x  =  t/tx  and  observe  that  0t 2  —  x  —  6  —  wtT;  then 


ST(fl,u,,/?) 

d0 


-«/' 


x 2  cos[#  +  u)trx  +  (x  —  6  —  utw)x 2]  dx. 


(3.104) 


(3.105) 


To  tightly  bound  the  above  integral,  find  its  supremum  by  observing  that  from 
wt*  >  0  follows  0  +  u)t„x  +  (v  —  6  —  w t„ )x 2  >  jtx2  on  x  €  [0, 1],  from  which 


l  f' 


6  +  uit^x  +  (x  —  6  —  wtw)x 


1  dt<L 


x2  cos(?rx2)  dx  (3.106) 


upon  noting  that  the  cosine  argument  is  always  within  the  region  [0,rr],  a  region 
where  the  cosine  decreases.  Substitution  u  =  x2  yields 


j: 


x2  cos(  7TX2  )  dx  =  2 


< 


Jo 

I  cos( 
Jo 


3  cos(xu)  du 


iru)  du  =  0 


(3.107) 


(3.108) 


after  observing  that  cos(xu)  <  cos (xu)  on  u  £  [0,1].  The  function 

T(6,u/,0)  is  decreasing  with  respect  to  0  >  0,  since  (3.105),  (3.106),  and  (3.108) 
imply  that  the  partial  derivative  of  T(0,u>,0)  with  respect  to  0  is  negative. 

The  second  preliminary  fact,  that  T (6,wy0)  decreases  with  respect  to  ui, 
will  be  proven  similarly  by  applying  the  chain  rule  to  compute 


6T(0,u>,0) 

6ui 


-S" 


t  cos (6  +  uit  +  0t2)  dt. 


(3.109) 


Again  let  x  =  </<„,  observe  that  0t2„  =  x-6-u)tr,  and  apply  the  same  reasoning 
leading  to  (3.106);  then 

*T(fl,w,/?)  <ja  [' 


du) 


t2  /  x  cos( xx2)  dx  =  0.  (3.110) 

Jo 
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I 

I 


The  partial  dei'v.-tive  ST(9,ui,  0)/Sui  <  0  implies  the  promised  result. 

The  third  preliminary  fact,  that  T (ui,9,0)  decreases  with  respect  to  9 ,  fol¬ 
lows  because  9  =  n  —  uitz  implies 

9  —  cos  0  4*  1 

T(0,w,O)  =  /  sin(0  +  uit)  dt  = - .  (3.111) 

Jo  w 

The  fact  that  9  £  [0,  ir),  a  region  on  which  the  cosine  decreases,  clearly  demon¬ 
strates  that  T (9,u>,0)  also  decreases. 

Now  with  the  preliminary  facts  established,  consider  f\os(9  +  uit  +  0t2)dt 
and  let  {rfr}"=1  be,  in  order,  the  zeros  of  the  integrand  on  (a,  6).  In  other  words, 
a  <  <  j]i  <  6  and,  whenever  t  £  (a,  6),  cos (9  -\-ujt  +  (it2)  =  0,  if  and  only  if 

t  —  rji  for  some  i  (if  there  are  no  zeros  then  set  r/i  =  6).  Decompose  the  integral 
into  subintegrals  between  the  zeros  and  get 

j  f  cos (9  +  u it  +  0t2)  dt 

j  J  a 

where 

'  =  T  (9  -  |,w  +  20a,/?)  for  i  =  0 

'  =  T(0,w  +  2 0^,0)  for  0  <  i  <  n  (3.113) 

.  <  T  (0, ui  +  20r]n,0)  for  i  =  n. 


1=0 


(3.112) 


Since  the  t^s  are  an  alternating  sequence  of  elements  whose  magnitude,  after 
the  first  element,  decreases, 

<  max(eo, ei)-  (3.114) 


f 


cos (9  +  u/t  +  0t2)  dt 


Use  the  fact  that  T (9,ui,0)  decreases  with  respect  to  0  to  maximize  eo,ei 
by  putting  0  =  0.  Next,  maximize  with  respect  to  the  other  arguments  to  show 

max(e0,  ei)  <  T(0,u/  +  2 0a,  0)  =  - — .  (3.115) 

ui  +  10a 

The  conclusion  of  the  lemma  now  follows  from  (3.115)  and  (3.114). 

Definition  1  Given  the  partition  V  =  [a  =  <0  <  t\  ■  ■  tn_i,(n  =  6],  the  positive 
variation  of  x(t)  on  [a,  6]  is 

n 

Pz  =  supy'[ar(t1)  -  x(f,_,)]+  +  x(a)+  +  x(b)+  (3.116) 

V 

where  r+  has  the  value  r,  if  r  >  0,  and  the  value  zero,  otherwise. 

Definition  2  Given  the  partition  V  =  [a  =  t0  <  t\  ■■  ■<„_!,<„  =  6),  the  nega¬ 
tive  variation  of  x(t)  on  [a,  6]  is 

n 

Nz  =  inf  ][>«.)  -  x(t,_1)]~  +  x(a ) ”  +  x(6)“  (3.117) 

y  1=1 

where  r~  has  the  value  r,  if  r  <  0,  and  the  value  zero,  otherwise. 
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A  function  is  said  to  be  of  bounded  variation  if  both  its  positive  and  negative 
variations  are  finite. 


Under  the  restrictions  that  b  >  a,  w  >  0,  u>  +  2 0a  >  0,  and  that  x(t)  and 
y(t)  are  piecewise  continuous  and  of  bounded  variation  on  [a,  6], 

fb  IP  P 

/  cos(0  +  ud  +  t)t2)x(t)y(t)  dt  <  (3.118) 

Ja  w  +  2/Ja 

where  PT  and  Py  are  the  pc  it'-  nations  of  x(t)  and  y{t). 


Since  the  function  x(t),  '••-'ded  variation,  it  is  integrable  implying 

that,  for  arbitrary  c  >  0,  there  (  ..  ..cp  function  x,(t)  =  x(f,_i 

with  corresponding  part  it  ion  [a  —  <  t\  •••<„_  i ,  tn  =  6),  such  that 


-  x,(/)|  dt  <  e. 


(3.119) 


'I'he  step  function  represented  above  is  the  sum  of  nonoverlapping  steps.  We 
want  to  reconstruct  it  as  the  sum  of  overlapping  steps  with  the  property  that 
the  accumulated  absolute  amplitudes  of  the  steps  are  minimal.  We  do  this  in 
an  iterated  fashion  by  ordering  the  step  amplitudes  p0  =  x(Lu)  <  ■  p,  _  i  = 
x(til_l ),  }i]  =  0,  ;>.  +  1  =  xff,  )  •  •  •  <  pn  —  x(t, n)  with  the  zero  amplitude 
included  and  defining  the  increments  I-  =  [f’j-\<P])-  Starting  with  the  Urth 
step,  an  increment  l}  is  considered  '  ■  -i”  if  l}  £  [0,  x(  t, )]  and  “closed” 

otherwise.  Whenever  an  increment  l}  transitions  from  open  to  closed,  define 
a  step  of  amplitude  rk  }  =  C(/;)  and  of  duration  d*  y  ranging  from  f*  to  the 
time  when  the  increment  was  last  opened.  Proceeding  in  this  manner,  the  step 
function  now-  has  the  form 


.1  (3.120) 

J=0  *■  =  0 

upon  setting  rt }  =  I)  for  previously  undefined  values 

At  each  stage  in  the  iteration,  notice  that  the  sum  of  the  lengths  of  incre¬ 
ments  either  opened  or  closed  is  equal  to  the  variation  of  the  step  function  at 
that  point .  imply  ing  t  hat  rt  ;  =  x(t}  )  —  x(t}  _  |  ).  Not  ice  also  t  hat  an  incre¬ 

ment  opened  by  an  im  Tease/ decrease  in  tin-  step  function  can  be  closed  only  by 
a  future  opposite  derrease/mcrease.  meaning  every  step  is  uniquely  associated 
with  a  point  of  increase  1  In-se  facts  mean  that 

fl  »* 

rr  ^,W/,  )-/■(/._,)]♦  4  [xf/„))+  (3.121) 

;  ~0  *  -  0  i-l 

r  PT  (3122) 

w  here  t  lie  last  relat  ion  follows  from  lie  deli  ml  ion  of  posit  ive  variat  ion  Properly 
define  c/  arid  T\  m  terms  of  rt  1 1  and  </;-  ,  to  write 

m 

'Jr, 

I  .  I 


-  i  T-J 


m) 
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with  <  Px  and  fa  |x(<)  —  <  c.  A  similar  step  function  exists 

for  y(t),  namely, 


y»(t)  ~  b< /|u,-|,ti|] 


(3.124) 


i=t 


with  Ei=i  b‘  <  py  and  /afc  l»(0  -  y>(t)W  < <■ 

With  these  two  step  functions  in  hand,  compute 

f  sin (9+ut+0t2)x(t)y(t)  dt  -  f  sir \(6  +wt+  fit2  )x  ,(t)ys{t)  dt+e(t)  (3.125) 
J  a  J  a 

where  |e(<)|  <  ((A/^i  +  MM)  +  c2,  M\x\  =  sup,  |x(<)|  <  oo,  and  M\y\  - 
sun,  |y(<)|  <  oo.  Putting  (3.123)  and  (3.124)  into  (3.125)  yields 

m  n  b 

'%2J2bkC‘  sin(^  +  w/+/?/2)/[Tl_,,r,jn[Uk_,,uk]  dt.  (3.126) 

k  =  l  t=l  Ja 

Applying  Lemma  1  and  maximizing  the  bound  by  replacing  the  starting  time 
for  each  step  with  the  worst  case  a  forms 


[  cos (0  +  tut  +  l3l2)x(t)y(l)  dt  < 
J  a 


9P  P 
u  +  2/fa 


(3.127) 


after  noting  that  c/  <  Pr  and  Yi’k=  l  bk  <  Py  bet  i  — *  0,  then  |e(f)|  —  0 

and  the  lemma  is  proved. 


3.C  MOMENTS  BETWEEN  SQUARES  OF 
CORRELATED  GAUSSIAN  RANDOM 
VARIABLES  WITH  RANDOM  PHASE 
COMPONENT 


We  have  in  this  section  two  random  variables 

P  =  o  cost?  +  /?sin  0  +  v  (3.128) 


Q  =  7  cos  9  4-  6  sin  9  +  y 


(3.129) 


where  the  random  variable  6  is  uniformly  distributed  on  the  on  [0,2t]  and  the 
Gaussian  random  variables  v  and  t)  are  zero-mean  with  covariances  a2,  <r2,  and 
<r2n.  We  want  to  compute  the  mean  and  variance  of  P2,  the  mean  and  variance 
of  P2,  and  the  covariance  between  P2  and  Q2 . 

Our  calculations  will  be  assisted  by  a  formula  of  the  general  fourth  moment 
between  the  Gaussian  random  variables  xq,  i),  rn.  and  *3  with  means  m,  = 
£(x, )  and  covariances  <r2  =  £[(x,  —  m,)(x}  —  m; )].  The  formula  is 


£(xnxix2x3) 


^,4,  +  4-  o2yAc\2  4  nfn  m-jrri:,  4-  <722m,m3 


+ 


<T23r«irn2  +  (T22m„»ri3  +  <rf:,rnofn-j  4  tr2:im„mi 


4-  rrii/tn  1  m2  rn-. 


(3  130) 
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We  will  now  compute  the  mean  of  P2  and  Q7 .  Equation  (3.128)  IliiptiV/O 
£[P2/9]  =  a2  cos2  6  +  02  sin2  8  +  2o/?cos  8 sin  8  +  <J2  (3.131) 

which,  upon  averaging  over  9,  becomes 

e[P7]  =  y  +  y  +  ^2-  (3.132) 

Similarly, 

£[<?2]=y  +  y  +  <  (3-133) 

Onward  to  the  covariances.  Equations  (3.128),  (3.129),  and  (3.130)  imply 
that 

£[P2Q2/8)  =  <x2<r2  +  2<74r?  +  <t2(7  cos  0  +  6  sin  8)2 

+  4aln(ct  cosd  +  /?sin  8)(j  cos  8  +  6  sin  9)  +  <r2(a  cos8  +  /?sin  0)2 

+  (a  cos  9  +  /?  sin  9)7(~(  cos  9  +  6  sin  0)2  (3.134) 

which,  upon  averaging  over  0 ,  becomes 

£[P2Q 2]  =  ^2  +  2^  +  ^(72  +  ^)+^(«2  +  /?2)  +  2^VaT  +  W 

+  ^(q272  +  07b2)  +  l(7202  +  4a0yb  +  a2^2).  (3.135) 

o  o 

With  (3.132),  (3.133),  and  (3.135),  we  conclude  that 
cov[P2,Q2]  = 

2<  +  2a2„(a7  +  fib)  +  ~(a272  +  02b2  -  y202  -  n7b7)  +  ^a/?7*(3.136) 
which  specializes  to 

var[P2]  =  2^  +  2(T2(o2  +  /?2)+^(a4  +  /?4)+  -a207  (3.137) 

o  4 

var[Q2]  =  2<r4  +  2rr2(72  + /i2)  +  i(74 -(■  A4)  +  i7262.  (3.138) 
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3.D  DERIVATION  OF  UPPER  AND 
LOWER  BOUNDING 
DISTRIBUTIONS  FOR  THE  SUM  OF 
SQUARES  OF  CORRELATED 
GAUSSIAN  RANDOM  VARIABLES 
WITH  RANDOM  PHASE 
COMPONENT 


Theorem  1  Define  the  n-dtmensional  vector  Pr  =  [pi,  •••,?«]  to  have  com¬ 
ponents 

Pi  ;=  ai  cos  0  -f  b,  sin  0  -f  Wi,  1  <  i  <  n  (3.139) 

where  each  ai  and  bi  is  a  constant,  6  is  a  uniformly  distributed  random  variable 
on  [0,2n],  and  { t  is  a  sequence  of  zero-mean  Gaussian  random  variables 
with  an  invertible  covariance  matrix  £ .  Then,  with  P  defined  as, 


l -Q* 

\J ^max  j  \) 

/  * 

\ 

<  Pr\pTP  <  it]  <  1-Q* 

\/€ min  >  \j 

1  * 

\ 

V 

''max 

L  J 

V 

/'min 

(3.140) 


where  Amax,  Amin  are  respectively  the  maximum  and  minimum  eigenvalues  of  £ 
and  where  emax,  emm  are  respectively  the  maximum  and  minimum  eigenvalues 
of  the  matrix 


D  =  AT£A 


(3.141) 


with 


(3.142) 


and,  finally,  where  Qm  is  the  generalized  Marcum  Q-function  defined  as 


QmW)* 


(3.143) 


Proof 

The  proof  consists  of  three  parts.  In  the  first  part,  the  conditional  proba¬ 
bility  Pr[PTP  <  k/0]  is  expressed  as  an  integral  of  a  multidimensional  Gaus¬ 
sian  density  over  a  spheroid  centered  at  the  origin  Upon  transformation  with 
a  decorrelating  matrix,  the  region  of  integration  becomes  ellipsoidal  and  the 
Gaussian  density  becomes  independent  with  each  of  its  marginal  densities  hav¬ 
ing  unity  variance.  Then  through  eigenvalue  analysis,  the  ellipsoidal  region  of 
integration  is  inscribed  and  circumscribed  with  spheroids  yielding  corresponding 
bounds  on  the  integral.  The  second  part  of  the  proof  shows  that  the  integral  of 
an  independent  Gaussian  distribution  over  an  arbit  rary  spheroid  depends  only 
on  the  magnitude  of  the  mean  vector  and  decreases  with  respect  to  it.  This  fact 
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will  enable  further  bounding  in  the  third  part  after  computing  the  minimum 
and  maximum  of  the  mean  as  a  function  of  9.  An  aftereffect  will  be  the  removal 
of  9  dependence  in  the  bounds  allowing  their  direct  application  to  the  uncondi¬ 
tional  probability  Pr [PT P  <  £].  Next  the  bounds,  which  are  still  expressed  as 
integrals,  are  evaluated  in  closed  form  via  the  generalized  Marcum  Q-function. 


Part  I 

By  applying  the  expression  for  a  multivariate  Gaussian  density,  the  condi¬ 
tional  probability 


P r[P?P  <  k/S]  =  — ‘  /  e-i(/,-AC]Tr-[P-AC]  dp 

2x1171*  Jli 


2w|27|*  J[PTP<k\ 


where 


C  = 


cos  9 
sin  6 


(3.144) 

(3.145) 


The  matrix  £~  *  with  the  property  £~  ^  £~  *  =  X-1  is  guaranteed  to  exist, 
since  X  is  an  invertible  covariance  matrix.  Furthermore,  for  the  same  reason, 
there  exists  an  orthogonal  matrix  T,  such  that  £  =  Tt AT ,  where  A  is  a 
diagonal  matrix  of  eigenvalues  of  £.  We  can  now  define  the  transformation 
X  —  T£~?P,  from  which  follows 


Pr [PTP<k/0]=~f  dX  (3.146) 

2  *J{XTAX<k\ 

where 

Mg  =  T£~*  AC.  (3.147) 

Now,  since  £  is  an  invertible  covariance  matrix,  each  entry  of  A  (i.e.  eigenvalues 
of  £)  is  positive.  Hence 

[AmaxXTX  <  Jfc]  c  \XTAX  <  jfc]  C  [An„nXTX  <  Jfc]  (3.148) 

from  which 

Pr  [PTP<k/0}>—[  e-\[X-Me\T[X-Me}  dx  (3  149) 

2 *  J[XTX<kixm.,} 

and 

Pr[PrP  <  kJO]  <  /  e-| [X-Me]T[X-Me]  dX  (3  150) 

~  1  ~  ^  J[XTX<k,xmn] 


Part  II 

We  aim  to  show  that,  for  a  given  r, 


rti-L  [  e-hlx-M)r{X-M]dx 

2»  J  [  X  7  X  <  r  | 


(3.151) 


depends  only  on  the  magnitude  of  M  and  decreases  as  \M\  gets  larger. 
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There  exists  an  orthogonal  matrix  U,  such  that  UM  -  [jAfj.C,  ,0]T. 
The  matrix  U  is  simply  a  change  of  orthonormal  basis  to  one  that  includes 
M/\M\  as  its  first  member.  Now  Y  —  UX  transforms  (3.151)  to 


-i/r 


[■ YTY<r ] 


(3.152) 


As  promised  for  a  given  r,  (3.152)  depends  only  on  the  magnitude  of  M ,  hence 
the  notation  /?(|Af|). 

We  now  show  that  R(\M\)  decreases  with  respect  to  the  magnitude  of  M 
by  showing  that,  for  any  positive  increment  AjAT  j,  the  corresponding  difference 

AftflAf |)  =  tf(|M|  +  A|Af|)  -  R(\M\)  is  negative. 

Make  the  respective  substitutions  zy  =  yy  —  |Af  |,  z;  —  yy,  for  2  <  i  <  n,  to 
R(\M\)  and  Zj  =  yy  -  \M |  -  A|M|,  z,  =  y.,  for  2  <  i  <  n,  to  ft(|M|  +  A|M(). 
Then 


AR(\M\)=~  [  e  -z  z  dz  -  2-  f  e-*zTz  dZ 

2ir  JG  2ir  JH 


where  the  sets  G  and  //  are 


G  = 


Z  :  (z,  +  |M|  +  A|M |)2  +  jS?  <  r 


i  =  2 


(3.153) 


(3.154) 


H  = 


Z:(z,  4  |Mj)2  4  ]T  *1  <  r  ■ 

i  =  2 

Cancel  out  the  common  points  of  G  and  //  then 

ar(\m\)  =  ^-[  e-?zTzdz-^-[  e-^7,  dZ. 
2*  Ja~n  2tt  Jh~g 


(3.155) 


(3.156) 


Let  ui  =  — zi  —  2|M|  —  A|M|,  vt  =  z, ,  for  2  <  i  <  n,  in  the  second  integral 
above;  then  //  ~  G  is  mapped  to  G  ~  //  and  the  integrals  can  be  combined  to 
yield 


A/e(|M|)  = 


(3.157) 

The  coordinate  z i  is  in  G  ~  //,  if  and  only  if 

n 

(z,  +  \M  |)2  4  2A|Af|(zi  4  |M|)  +  (A|Af  |)2  4  £  z(  <  r  (3.158) 

*  =  2 

and 

n 

(z,  +|M|)2  +  £z.  >r.  (3.159) 

1  =  2 
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Upon  subtracting  both  relations,  Z\  must  satisfy 

z\  <  —\M\  —  — i^— -I -  (3.160) 

Using  this  relationship  in  (3.157)  implies  that  A/7(|A/|)  <  0,  meaning  that 
f?()M|)  decreases  with  increasing  magnitude  of  M. 

Part  III 

Further  bound  (3.149)  and  (3.150)  by  respectively  maximizing  and  mini¬ 
mizing  magnitude  of  the  mean  over  the  random  phase  6.  From  (3.147),  the 
magnitude  of  the  mean  is 


|M«|2  =  Me1  Ms  =  CtAtS~1AC 


(3.161) 


following  from  the  fact  that  T,  being  orthogonal,  satisfies  TtT  =  I.  The  matrix 
At A,  being  symmetric,  ensures  that  it  can  be  diagonalized  making 


|Af#|2  —  [cos#,  sin#]!/7 


0 


U 


cos  6 
sin  0 


(3.162) 


where  U  is  an  orthogonal  matrix  and  emax  and  emm  arc  the  eigenvalues  of 
At E~'A.  Since  U  is  orthogonal,  it  rotates  the  plane  by  some  angle  <p.  This 
means  that 


U 


cos  0 

cos(#  —  <j>) 

sin  0 

sin(0  -  <f>) 

and  hence  that 


\Me\7  =  emax  cos J(0  -</>)  +  emm  sin2(#  -  <{>). 


(3.163) 


(3.164) 


The  eigenvalues  ernax  and  cmm  are  nonnegative  since  AT 1 7  1 A  has  a  square 
root,  namely,  27~5  A.  This  fact,  along  with  (3.164),  implies 


<  \M,\  ^  \f  e  m  ax .  (3.165) 

Use  (3.165)  and  (3.152)  to  deduce  from  (3.150)  and  (3.149)  the  bounds 

g-  Vi  n/' '■*! )2  JY 


V'[rTp < n  >  if 


YTY<k/xm.,  | 


Pt[PtP  <  k]  < 


hi 


YTY  <k/xm,„ 


e~  5  v?  - 


(3.166) 

dY. 


(3.167) 

The  integrals  in  the  above  bounds  are  simplified  by  showing  that  they  are 
the  distribution  functions,  evaluated  respectively  at  kf Amax  and  kj Am,n  of  the 
sum  of  n  non-central  \2  random  variables  with  noncentrality  cmax,  emm.  An 
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explicit  expression  for  this  distribution  is  given  in  Torrieri2  and  leads  to  the 
conclusion  of  the  theorem: 


i-o* 

max  »  \  l 

1  k 
\ 

<  Pr  \PTP  <  jfcl  <  1  -Q* 

y{ prmn  j  \j 

/  * 

\ 

V 

''max 

L  J  2 

V 

''umn 

(3.168) 
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CHAPTER  4 

DETECTION  VIA  THE  AMPLITUDE 
DISTRIBUTION  FUNCTION 

4.1  BACKGROUND  AND  MOTIVATION 

The  unfriendly  detection  and  interception  of  secure  communications  is  a 
topic  of  much  current  research.  Secure  communications  usually  involve  some 
variety  of  spread-spectrum  modulation,  whose  purpose  is  to  add  ambiguity  or 
“randomness”  to  the  communication  waveform  as  a  measure  against  unintended 
detection  or  interception.  The  usual  procedure  for  randomizing  the  waveform  is 
pseudo-random  variation  of  transmission  times  (time  hopping),  phases  (direct 
sequence),  or  frequencies  (frequency  hopping).  The  development  of  a  method 
to  detect  frequency-hopped  waveforms  is  the  subject  of  this  chapter. 

The  use  of  the  Amplitude  Distribution  Function  (ADF)  for  detection  is 
a  new  idea  with  potentially  many  diverse  applications.  However,  although 
this  work  focuses  solely  on  the  detection  of  frequency-hopped  waveforms,  the 
general  ADF  technique  can  be  applied  to  related  areas  such  as  radar  or  sonar 
detection.  The  central  idea  of  the  technique  is  that  the  ADF  of  an  observed 
signal  in  additive  noise  is  the  convolution  between  the  individual  ADFs  of  signal 
and  noise.  We  have  shown  through  the  use  of  deconvolution  techniques  that 
the  signal  component  can  be  separated  and  thus  detected  even  for  small  signal 
levels. 

There  are  previous  works23,24'25'26  that,  in  essence,  use  the  ADF  but  none 
has  given  a  precise  definition  and  mathematical  development  like  those  offered 
here.  Moreover,  to  our  knowledge  there  is  no  reference  that  directly  uses  the 
ADF  idea  for  detection. 

The  ADF  indicates  the  time  fraction  that  a  waveform  is  below  a  given  am¬ 
plitude,  much  like  a  probability  distribution  function  measures  the  probability 
that  a  random  variable  is  below  a  given  value.  Previous  researchers  have  used 
this  concept  but  failed  to  give  a  precise  definition  of  the  ADF  as  it  applies  to 

23R.F.  Paw  ula,  S  O.  Rice,  and  J.H.  Roberts,  “Distribution  of  the  Phase  Angie  Between 
Two  Vectors  Perturbed  by  Oaussian  Noise”,  IEEE  Trans.  C'ommun.,  vol.  COM-30,  pp. 
1828- 1840  August  1982. 

2*J.H.  Higbie,  “Adaptive  Nonlinear  Suppression  of  Interference”,  MILCOM  *88  ConJertn.ee 
Proceedings ,  23.3.1  21.3.9,  October  1988. 

2I,S.O.  Rice,  “Statistical  Properties  of  a  Sine  Wave  Plus  Random  Noise,"  Bell  Syst.  Tech. 
J..  pp.  109  157,  January  1948. 

26 .1.  Salz,  and  S.  Stein,  '‘Distribution  of  Instantaneous  Frequency  for  Signal  Plus  Noise", 
IEEE  Trans.  Injorm.  Theory,  vol.  IT-10,  pp.  272  274,  October  1904. 
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Figure  4.1:  Definition  of  ADF 


both  deterministic  and  stochastic  signals.  For  signal  X(t),  our  definition  is 

Fx(a)  =  lini  ~S{C[t  :  X(t)  <  a,  0  <  t  <  T]}  (4.1) 

where  C  is  set  function  givmg  length.  Figure  4.1  illustrates  that  the  ADF  is 
simply  the  time  fraction  that  X(<)  is  below  a  given  threshold  a.  With  this  def¬ 
inition,  we  have  proved  that,  under  very  general  conditions,  the  ADF  of  signal 
plus  noise  is  the  convolution  of  the  signal  ADF  and  the  noise  ADF  individually. 
This  result  would  not  have  been  possible  without  a  definition  that  applied  to 
both  deterministic  and  stochastic  signals. 

The  adf  or  amplitude  density  function  is  the  density,  if  it  exists,  implied 
by  the  ADF.  There  is,  of  course,  a  corresponding  convolutional  relationship 
between  the  adf  of  signal  and  noise  and  the  individual  adf’s  of  signal  and  noise. 
Some  examples  exemplify  this  convolutional  relationship  and  hint  at  the  poten¬ 
tial  of  the  ADF  in  signal  detection.  Figure  4.2  shows  the  ADF  of  a  modulated 
sine  wave  (the  signal),  the  ADF  of  noise,  and  the  ADF  of  signal  plus  noise. 
The  main  point  here  is  that  the  ADF  of  a  sine  wave  is  invariant,  under  most 
phase  and  frequency  modulations,  but  these  are  exactly  the  modulations  used 
to  thwart  a  potential  interceptor.  Therefore,  the  most  typical  spreading  modu¬ 
lations  will  not  degrade  the  performance  of  an  ADF  detector.  Also  of  interest 
is  the  complex  ADF,  that  is,  the  two-dimensional  distribution  in  amplitude  of 
the  signal’s  in-phase  (I)  and  quadrature  phase  (Q)  components  (Figure  4.3). 

Both  of  the  last  two  examples  illustrate  the  convolutional  spreading  of  the 
ADF  due  to  additive  noise.  Like  a  photograph  taken  while  the  camera  is  out  of 
focus,  noise  smears  the  signal  part  of  the  ADF.  By  a  technique  borrowed  from 
image  processing,  the  picture  (i.e.,  ADF)  can  be  refocused  to  reveal  the  un¬ 
derlying  picture  detail  (i.e.,  signal).  Deconvolution,27  as  this  process  is  called, 

J7H.C.  Andrew*  and  B.R.  Hunt,  Digital  Image  Retloration,  Prentice-Hall,  Englewood  Cliff*, 
New  Jersey,  1977. 
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involves  convolving  the  picture  with  a  kernel  function  k(x)  that  has  been  spe¬ 
cially  constructed  from  noise  details.  Figure  4.4  shows  the  result  of  the  process, 
an  out-of-focus  picture  before  and  after  deconvolution. 

Our  original  idea  allows  the  direct  application  of  the  deconvolution  tech¬ 
niques  to  signal  detection.  Suppose  a  noisy  signal  X(<),  observed  over  the 
interval  T ,  is  transformed  into  the  function  Fx(x)  by 

Fx(x)=^j\[x-X(i)}dt  (4.2) 

where  k(x )  is  the  deconvolution  kernel.  We  have  shown  that,  for  large  T, 
Fx(x)  converges  directly  to  the  deconvolved  ADF.  For  this  reason,  we  call 
this  the  “Deconvolution  Statistic”.  Furthermore,  we  have  shown  that  samples 
of  the  Deconvolution  Statistic  are  approximately  jointly  Gaussian,  to  which 
well-known  optimal  detection  techniques  apply. 

In  summary,  our  approach  consisted  of  precisely  defining  the  ADF  and  prov¬ 
ing  the  existence  of  an  intuitive  relationship  between  signal  and  noise,  that  of 
convolution.  Borrowing  techniques  from  image  processing,  we  showed  that  the 
effects  of  the  noise  could  be  separated  from  the  signal  by  a  process  called  de- 
convolution.  By  transforming  the  observed  waveform,  we  generated  a  random 
process  that  converges  directly  to  the  deconvolved  ADF  and  upon  which  stan¬ 
dard  detection  techniques  apply. 

4.2  MATHEMATICAL  TOOLS  FOR  THE 
ADF 


The  ADF,  as  defined  here,  is  an  original  concept  and  thus  needs  a  firm 
mathematical  foundation.  This  section  precisely  defines  the  ADF  in  a  way  that 
applies  equally  well  to  deterministic  and  stochastic  signals.  This  basic  definition 
is  extended  to  include  the  concepts  of  a  joint  ADF  and  the  notion  of  amplitude 
independence,  a  notion  analogous  to  independence  in  probability.  By  way  of  a 
sequence  of  lemmas  and  theorems,  two  significant  results  are  established.  The 
first  is  the  already  promised  result  that  the  ADF  of  signal  plus  noise  is  the 
convolution  of  the  signal  ADF  with  the  noise  ADF.  This  is  proved  under  the 
very  general  constraint  that  the  second  derivative  of  the  noise  autocorrelation 
exists  and  is  finite  at  time  difference  zero.  The  remaining  result  of  significance 
is  a  linkage  between  the  ADF  and  the  instantaneous  probability  distributions 
of  the  signal  plus  noise.  We  begin  with  a  precise  definition  of  the  ADF. 

Definition  3  The  Amplitude  Distribution  Function  (ADF),  written  Fx(a)  for 
a  stochastic  process  X(t),  is 

Fx(a)=  lim  )-€  {C  [t  :  A'(<)  <  a,  0  <  t  <  T]}  (4.3) 

T  — *  oo  1 

where  C  is  a  set  function  giving  length.  Additionally ,  the  limit  must  exist  for 
all  a. 

If  the  signal  is  deterministic,  then  the  definition  of  the  ADF  reduces  to 

FX{a)  =  lim  ^  {£[<  :  A'(f)  <  «,  0  <  t  <  7’]}  (4.4) 

7  ’  —  'lo  J 
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The  ADF  is  not  a  distribution  function  in  the  strict  sense  because  its  extreme 
values  are  not  necessarily  one  or  zero  and  it  may  not  be  right-continuous.  For 
example,  the  function 


for  0  <  t  mod  1  <  1/2 
for  other  t 


(4.5) 


has  ADF  identically  equal  to  1/2.  Consider  also  the  function 


Its  ADF  is 


for  0  <  t  mod  1  <  1/2 
for  other  t. 


Fs(a) 


0  for  oo  <  a  <  —  1 
1/2  for  —  1  <  a  <  1 
1  for  1  <  a 


(4.6) 


(4.7) 


which  is  right-continuous  at  —1  and  left-continuous  at  1.  The  ADF,  not  neces¬ 
sarily  being  a  true  distribution,  creates  problems  in  situations  which  require  an 
ADF-induced  measure,  for  instance,  the  Lebesgue-Stieltjes  integral.  In  these 
cases,  we  use  the  right-continuous  extension  of  the  ADF,  defined  as 


F*(a)  =  lim  F(r). 

r— a  + 


(4.8) 


Analogous  to  the  joint  probability  function  of  random  variables,  there  exists 
a  joint  ADF  of  two  different  stochastic  processes,  defined  as  follows. 

Definition  4  The  joint  ADF,  written  FXy  (a,  b)  for  stochastic  processes  X{t) 
and  Y(t),  is 


Fx  v(a,  b)  —  lim  ~£  {£  [t  :  X(t)  <  a  and  Y(t)  <  6,  0  <  t  <  71} .  (4.9) 

This  definition  will  be  used  to  define  the  following  concept  of  amplitude  inde¬ 
pendence,  analogous  to  that  of  independence  between  random  variables. 

Definition  5  Two  stochastic  processes  X(t)  and  Y(t)  are  amplitude  indepen¬ 
dent,  if  their  joint  ADF  is  the  product  of  the  ADF  for  each  process.  In  other 
words, 

Fx.y(a,b)  =  Fx(a)Fy(b).  (4.10) 

We  now  will  show  a  relationship  between  the  ADF  of  stationary  Gaussian 
noise  and  its  instantaneous  probability  distribution,  thus  enabling  the  estab¬ 
lishment  of  more  directly  applicable  results. 


Lemma  3  Let  Y(t)  be  a  stationary,  zero-mean,  Gaussian  process  with  auto¬ 
correlation  R(t),  such  that  —R"( 0)  <  oo;  then,  for  any  measurable  set  A, 


£{C[t  :  Y{t)  <  a,  t  £  A))  =  $ 


CA 


(4.11) 


where  $  is  the  distribution  function  of  a  standard  Gaussian  random  variable 
and  (To  =  y/Wi 
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The  proof  for  Lemma  3  is  included  in  Appendix  4. A. 

In  words,  the  above  lemma  means  that  the  average  time  that  the  noise 
process  is  below  the  threshold  a  on  the  set  A  is  equal  to  the  percent  of  time 
that  the  noise  process  is  below  a  at  any  single  point  times  the  length  of  the  set 
A.  This  result  implies  that  the  ADF  of  stationary  Gaussian  noise  is  identical 
to  its  instantaneous  distribution. 

With  the  help  of  the  previous  results,  we  now  can  prove  amplitude  inde¬ 
pendence  between  a  deterministic  signal  and  stationary  Gaussian  noise,  whose 
autocorrelation  has  finite  second  derivative  at  time  difference  zero.  This  result 
will  be  necessary  to  prove  the  convolutional  relationship  between  signal  and 
additive  noise. 

Theorem  2  Lei  S(t)  be  a  deterministic  signal  and  let  N(t)  be  a  stationary, 
zero-mean,  Gausstan  process  with  autocorrelation  R(t)  such  that  -/?"(())  <  oo; 
then  S(t)  and  N(t)  are  amplitude  independent.  Stated  symbolically, 

FsAa,b)  =  Fs(a)Fjvlb).  (4.12) 


Proof 

By  Lemma  3, 

£{£[<:  S(t)  <  a  and  N(t)  <  b,  0  <  t  <  T\) 


=  £{C[t:N(t)<b,t£S;ln[0,T)]} 

(4.13) 

=  Fs(b)C  [S”‘  n[0,T)] 

(4.14) 

where  5"1  =  [<  :  S(<)  <  «]•  Hence  the  joint  ADF  of  S(t )  and  N(l)  is 

Fs,N(a,b)  =  Fff(b)  lim  [S" 1  D  [0, T)] 

T  —♦  oo  i 

(4.15) 

=  Fs(a)FN(b) 

(4.16) 

implying  the  amplitude  independence  of  S(t)  and  N(l). 

We  now  can  prow  the  most  important  result  of  this  section,  that  of  convolu¬ 
tion  between  the  ADFs  of  signal  and  additive  noise.  The  idea  of  deconvolution 
and  the  deconvolution  statistic  rest  firmly  on  this  result. 

Theorem  3  Let  S(t)  and  N(t)  be  amplitude  independent  and  let  either  Fs  or 
F$  be  continuous;  then  the  ADF  o/Y(t)  =  S(t)  +  N{t)  is 

Fy(a)  =  r  Fs(a-s)dFZ(s) 

J -oo 

(4.17) 

=  f°°  Fs(a-n)dFs(n) 

J -OO 

(4.18) 

where  Fg  and  F ^  are  the  right-continuous  extensions  of  Fs  and  F s- 

The  proof  for  Theorem  3  is  given  in  Appendix  4.B. 

The  next  theorem  is  a  general  statement  about  the  ADF  of  signal  with 
additive  noise  and  its  instantaneous  distribution.  Alternatively,  this  theorem 
could  have  served  as  the  definition  of  the  ADF,  but  then  there  would  be  technical 
difficulties  in  determining  the  ADFs  of  purely  deterministic  signals. 


102 


The  Johns  Hopkins  University 

Applied  Physics  Laboratory 


Theorem  4  Let  y(<)  =  S(t)  +  X(t)  and  let  Fs  be  nght-continuous;  then,  if 
either  Fpj  or  Fs  continuous,  the  ADF  ofY(t)  »s 


Fy(a)  =  lim 

T— ►co 


FN[a  -  5(1)]  }  dt. 

The  proof  for  Theorem  4  is  included  in  Appendix  4.C. 


j;e{j%N[a-s(t))^ 


(4.19) 


4.3  DECONVOLUTION  STATISTIC 


As  shown  in  the  previous  section,  the  ADF  of  signal  plus  noise  is  the  ADF 
of  the  signal  convolved  with  the  ADF  of  the  noise.  What  was  exactly  shown  is 
that 


for  X(t)  =  S(t)  +  N(t)  where  N(t)  is  a  stationary  Gaussian  process 
with  autocorrelation  R(t)  satisfying  —R"( 0)  <  oo  and  where  S(t)  is 
a  deterministic  signal  with  defined  ADF;  then  F\  =  Fg  *  F^  where 
Fx ,  F$,  and  Fn  are  the  respective  ADFs  of  X(t),  S(t),  and  N(t). 

To  apply  this  result  in  the  construction  of  a  detector,  we  will  assume  from  herein 
that  the  above  restrictions  are  met  and  that  the  noise,  signal,  and  observations 

have  densities  defined  as /x  =  dFx(x)/di,  fs  =  dfs /da,  and  fs  =  dF^/da. 
These  densities  will  be  called  the  amplitude  density  functions  (adf).  We  will, 
for  reasons  explained  later,  make  the  restriction  that  the  noise  autocorrelation 
is  zero  after  some  duration  (i.e.  R(t)  =  0,  for  t  greater  than  some  T\ ).  This  is 
a  sufficient  but  probably  not  necessary  condition  for  the  asymptotic  statistical 
characterization  of  the  detection  statistic. 

If  the  adf  of  the  observed  signal  could  somehow  be  measured  or  estimated, 
then  by  deconvolution  the  signal  component  could  be  separated  from  the  noise 
component  and  hence  detected.  An  understanding  of  deconvolution  is  prerequi¬ 
site  to  understanding  of  how  the  deconvolution  statistic  effects  this  separation. 

In  general  to  deconvolve,  take  a  function  that  is  the  convolution  of  two 
different  functions  and  convolve  it  again  with  a  kernel  function  whose  net  effect 
is  to  undo  the  first  convolution.  As  such,  convolution  of  the  kernel  function 
with  the  original  convolving  function  should  result  in  a  delta  function.  The 
construction  of  a  kernel  function  with  this  property  is  usually  not  possible 
exactly  and  so  some  approximation  must  be  made. 

To  use  the  idea  of  deconvolution  in  a  detection  scheme,  the  adf  could  be 
measured  using  standard  density  estimation  techniques  and  then  convolved  with 
a  kernel  function  to  separate  signal  from  noise.  But  there  is  a  more  direct 
approach,  that  of  the  deconvolution  statistic 

fx (x)  =  k{x  -  A'(0]  dt  (4.20) 

where  fc(z)  is  the  deconvolution  kernel. 

The  usefulness  of  this  statistic  is  that,  as  T  gets  large,  it  converges  uniformly 
in  probability  to  the  desired  convolution  of  the  kernel  function  with  the  adf  of 
the  observations.  In  a  sense,  the  statistic  maps  the  observations  into  another 
domain,  the  amplitude  domain.  The  problem  of  detection  is  now  one  on  this 
new  dcr,.,»ii.  and  from  ii<-run,  it  is  approached  in  a  classical  manner. 
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4.3.1  Deconvolution  Kernel 

We  want  to  construct  a  kernel  with  the  property  that  when  it  is  convolved 
with  the  noise  adf,  a  function  closely  approximating  the  delta  function  results. 
The  approach  is  to  solve  the  equivalent  problem  after  Fourier  transformation, 
that  is,  to  find  a  function  that  when  multiplied  by  the  Fourier  transform  of 
the  noise  adf  yields  a  constant.  It  is  shown  below  that  this  can  be  done  only 
approximately. 

Proceeding  by  Theorems  3  and  4,  the  noise  has  adf 

Mi)=7ke'^-  (4-2,) 

where  <Tq  =  72(0) .  The  Fourier  transform  of  the  noise  adf  is 

/OO 

fN(x)&7xuz  dx  =  e~  *<,»<  (4.22) 

-OO 

From  the  convolutional  relationship  between  noise  and  signal,  the  transform 
of  the  observation  adf  is  Cx(w)  =  Cn(u)Cs(u).  From  this  relationship,  we 
might  be  led  to  construct  a  kernel  with  characteristic  function  1  /Cn(ui),  but 
this  proves  fruitless  since  the  inverse  Fourier  transform  does  not  exist.  Instead, 
if  before  Fourier  inversion,  we  multiply  1  /Cn(u)  with  window  function 

cos2  (y)  ,  ~<w<  J,  d  >  0  (4  23) 

we  then  get  a  family  of  kernels  kd(x)  indexed  by  d  with  the  property  that 
kd  *  Jn(z)  —  Hx)  as  d  — ' •  0,  hence  convolution  with  kd(x)  can  approximate,  to 
an  arbitrary  accuracy,  perfect  deconvolution. 

Application  of  the  above  strategy  with  equations  (4.22)  and  (4.23)  produces 

kd{x)  =  i  J  eHj)COs2  cos  dd.  (4  24) 

The  quantity  d,  which  controls  the  amount  of  deconvolution,  is  called  the  de- 
convolution  index. 

4.3.2  Statistical  Characterization 

In  order  to  set  up  the  detection  problem,  we  need  at  least  an  asymptotic 
statistical  characterization  of  the  deconvolution  statistic.  Specifically,  we  want 
to  show  that  samples  of  the  deconvolution  statistic  are  asymptotically  jointly 
Gaussian  for  large  T,  and  we  want  its  asymptotic  mean  and  variance. 

As  for  the  jointly  Gaussian  property,  consider  samples  of  the  deconvolution 
statistic 

«  =  fj0  k[*i-X[t)]<U  (4.25) 

for  some  finite  sequence  {x;}"=i  •  1°  order  to  prove  that  the  Zj’s  are  jointly  Gaus¬ 
sian,  it  is  sufficient  to  prove  that  c^z*  is  Gaussian  for  arbitrary  constants  C{. 
Rewrite 

n  l  m 

(4.26) 

.  j= l  j=i 
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where 


(4.27) 


(4.28) 

(4.29) 


It  is  easily  seen  that  the  /2;’s  are  independent  and  the  Iy-i's  are  independent. 
Hence,  for  large  T,  each  sum,  h j  and  HjLi  hj-i< 's  individually  asymp¬ 

totically  Gaussian.  Furthermore,  even  though  the  sums  are  correlated,  the 
overall  sum  is  approximately  Gaussian,  since  each  component  sum  is  Gaussian. 
Hence,  samples  of  the  deconvolution  statistic  are  jointly  Gaussian. 

As  for  the  mean, 


fi(x) 


(4.30) 


*[*-*(<)] 


(4.31) 


%  k  *  fx(x)  for  large  T 

with  the  last  step  following  from  Theorems  3  and  4. 

As  for  the  variance, 


<r2(r,y)  =  £{[/*(*)  -  p(x)][fx (y)  -  p(y)]} 


(4.32) 


(4.33) 


k(x  —  v)k(y  —  w)h(v,  w,  T)  dvdw  —  p(x)p{y) 


(4.34) 

where 

1  [T  fT  1 

=  75 1  i  2^Ae  *  isdi  <«5> 


A  =  [l-p2(s-<)]  (4.30) 


M  = 

v  -  5(s) 
w  —  S(t) 

(4.37) 

S  = 

[  1  -p(s  - 1) 

L  -M*  ->)  1 

(4  38) 

p(s-l)  = 

R(s  - 1) 

(4.39) 
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Figure  4.5:  ADF-Based  Detector 


We  want  to  show  that  the  variance  goes  to  zero  as  the  integration  time  T 
goes  to  infinity  from  whence  it  follows  that  the  detection  statistic  converges,  in 
probability,  to  its  mean.  To  this  end,  assume  T  >  T\  and  rewrite 


V'W'T)  =  A® 


.-A-Mt27M 


dsdt 


+ 


ML 


1  --Jr-MrSM 

— "  *1*  dsdt 


27T«ToA 


(4.40) 


(4.41) 


where  the  regions  within  the  square,  Q  =  [s,  t  :  0  <  s  <  T,  0  <  t  <  T\,  are 

Di  =  [s,t  :  s-t  <Ti)nQ  (4.42) 


D-i  -  (s,i:s-(>r,]nQ. 


(4.43) 


The  area  of  D\  is  less  than  2T7j  and,  for  iTq  >  0,  the  integrand  is  bounded  by 
l/2jr<To-  Furthermore  A  =  1  in  region  D2  Hence, 


*2T 

h{v,w,T)  <  -~  + 


2nalT2 


L 


ds 


l 


- ^7  fw/  —  S(  t  J]  ^ 

e  “I  dt. 


(4.44) 


Now  zissuming  k{x)  is  integrable,  we  have  from  (4.44)  and  (4.32)  that  cr2(x,  y)  — > 
0  for  large  T,  as  conjectured. 

The  fact  that  samples  of  the  deconvolution  statistic  are  asymptotically 
jointly  Gaussian,  and  the  expressions  of  its  mean  and  variance  will  be  necessary 
for  the  detector  development,  below. 


4.4  ADF-BASED  DETECTOR 


Using  the  fact  that  samples  of  detection  statistic  are  asymptotically  jointly 
Gaussian,  a  classical  detector  can  be  constructed  that  observes  the  amplitude 
domain.  We  assume  the  original  observed  waveform  to  be  of  the  form  X(t)  = 
n/25'  sin  wo<  +  n(i) ,  for  0  <  t  <  T,  where  n(t)  is  white  noise  of  spectral  density 
N0/2,  S'  is  the  average  signal  energy,  and  wo/2tr  is  unknown  frequency  in  the 
band  [fc  -  W/2,fc  +  W/2],  In  this  setting,  the  detection  problem  is  one  of 
choosing  between  Hq  (signal  absent)  and  H\  (signal  present  with  SNR,  7'  = 
S'T/No  >  7  =  ST/Nq).  Before  we  apply  our  detector,  we  filter  the  observations 
with  a  bandpass  tdter  to  produce  the  waveform  Y(t)  (see  Figure  4.5).  This  filter 
has  unity  gain  within  the  bandwidth  W  and  has  center  frequency  fc.  It  further 
has  response  such  the  noise  has  autocorrelation  R(t)  with  R(t)  =  0  for  /  greater 
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than  some  7)  .*  We  next  transform  the  output  of  the  filter  into  the  amplitude 
domain  via  the  deconvolution  statistic.  The  new  detection  problem  becomes, 
after  sampling,  one  of  deciding  between  the  presence  and  absence  of  a  signal 


given  {zi  =  /y(xj)}"=1  As  shown  earlier,  the  x,’s  are  jointly  Gaussian  and  have 
means  (4.32), 

.  ,  f  P(*i),  signal  absent 

/v(*0  =  {  ,  ,  ,  .  (4-45) 

1  p*r{Xi),  signal  present 

where 

P(*) 

1  i 

—  p  J  77^W 

y/2*N0W 

(4.46) 

r(x) 

=  —  \/25'  —  x2,  x2  <  25' 

it 

(4  47) 

and  covariances 

<V(x;,Xj) 

defined  by  (4.34)  with 

m  = 

Arolf/(l  -  Wt),  — 1  <  Wt  <  1 

(4.48) 

5(<)  = 

V2 S’  sin  uit. 

(4.49) 

The  R(t)  used  here  is  a  first  approximation  to  the  one  described  earlier.  Be¬ 
cause  of  this,  there  will  be  a  small  amplitude  variation  in  the  signal  with  respect 
to  frequency  that  will  be  ignored  in  the  analysis  to  follow.  The  above  analy¬ 
sis  completely  characterizes  the  asymptotic  joint  statistics  of  samples  of  the 
deconvolution  statistic. 

We  now  use  classical  methods  to  design  a  detector  on  the  amplitude  domain. 
The  philosophy  assumes  that  both  the  noise  and  signal  amplitudes  are  known. 
Then  under  this  condition,  an  optimum  detector  is  synthesized.  To  relax  this 
unrealistic  assumption,  we  will  assume  that  the  noise  level  is  known  or  mea¬ 
sured  and  that  the  signal  level  is  above  that  used  in  the  detector’s  synthesis. 
This  will  be  suboptimum  in  general  but  in  the  important  low-signal-level  case, 
performance  will  approach  the  optimum. 

We  can  now,  via  the  likelihood  ratio,21  define  the  optimal  test  statistic  as 

L  =  eIr^'Eq-  (4.50) 


where 


R-.  = 


^7(^l,*l)  •••  <r7(xi,x„) 


j  ^1  )  J 


(4.51) 


and 


E7  = 


-i  Pti1  l) 


(4.52) 


L 

In  order  to  evaluate  the  performance  of  this  detector,  we  note  that,  since 
Hy  is  nonnegative-definite  and  symmetric,  there  exists  a  matrix  such  that 

‘This  filter  is  not  realizable  exactly  but  can  be  approximated  to  an  arbitrary  degree. 
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I  I  _  _  T  __  _  I 

Ry  Ry  =  Ry.  Use  this  fact  and  the  diagonalization  of  Ry  7  R  Ry  2  = 
Ty  7Ayi7Tyi7  to  rewrite  the  test  statistic  as 

L  =  (G  +  My,o)TAy,o(G  +  My,0)-(G  +  My,7)TAy,7(G+ Afy,7)  (4.53) 
where 


G  = 


9 1 


L  92  J 


(4.54) 


with  {</,}  independent,  zero  mean,  and  unity  variance  and  where 


My.,y  i 


*^i,V 


■~n  _  X 

—  T'  /?  2 

—  ±  y,7  fly' 


Hy>(Xn)  -  Hy{X„)  J 


(4.55) 


We  note  that  (G+Afy  7)rAyi7(G+Afyi7)  is  asum  of  squares  of  independent 
Gaussian  variables.  We  now  have,  through  application  of  the  Berry- Esseen 
Theorem,22  that  this  term  is  approximately  Gaussian  distributed  with  an  error 
of  no  more  than  4 cja  where 


,  24m.\y,7  +  7 

c  -  max  A,  y  7  — — 5 - — 

*  ,y  ,y  4mf  y,  y  +  2 


(4.56) 


a2  =  ^  A?7,t7(4m?y  7  +2). 


(4.57) 


»=i 


For  each  particular  detection  problem  and  for  each  value  of  SNR,  this  error 
bound  determines  the  validity  of  the  CLT  argument.  Assuming  the  bound  is 
small,  then  L  itself  must  be  approximately  Gaussian,  being  the  sum  of  two 
Gaussian  distributed  random  variables.  The  distribution  and  hence  perfor¬ 
mance  of  L  is  determined  by  its  mean  and  variance  as  computed  below.  The 
mean  is 


My,y  =  Y,  [Vv.oKy.o  +  1)  -  A,,y.7(m2y7  +  1)]  (4.58) 

*=1 

while  the  variance  is 

n 

Vy,7  =  )  1  [2  (A?,y,o  —  ^i,f',y)  "b  4  (m',7',oAil7',o  -  mi>7',7  Aji7'i7)  j. 


1=1 


(4.59) 

Since  the  test  statistic  L  has  an  approximately  Gaussian  distribution,  the 
threshold  v  and  probability  of  detection  Pd,  for  a  given  probability  of  false 
alarm,  follow  as 

0=  y/V^*-'(l-  Pr)+Mo,y  (4.60) 

and 

pD  •  1 .  ♦  (4.6D 

where  <t»(z)  is  the  distribution  function  of  the  standard  Gaussian. 
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4.5  CONCLUSIONS 

The  ADF  was  introduced  and  shown  to  be  roughly  the  average  probability 
distribution  of  a  random  process.  Because  the  ADF  is  original,  a  mathematical 
foundation  was  laid  consisting  of  a  sequence  of  definitions,  lemmas,  and  theo¬ 
rems,  the  most  significant  of  which  was  the  fact  that  the  ADF  of  signal  plus 
noise  is  the  convolution  of  the  ADF  of  signal  and  the  ADF  of  noise  taken  sepa¬ 
rately.  The  technique  of  deconvolution  used  in  image  processing  was  the  germ 
from  whence  emerged  the  deconvolution  statistic,  a  statistical  transform  that 
mapped  a  stochastic  process  into  the  amplitude  domain.  It  was  shown  that 
by  proper  choice  of  a  kernel  any  degree  of  separation  between  the  noise  and 
signal  components  could  be  obtained.  For  the  particular  problem  of  detecting 
a  modulated  sinusoid  in  stationary  Gaussian  noise,  a  detector  was  developed 
around  the  detection  statistic.  The  detector’s  performance  was  analyzed  and 
compared  with  that  of  a  radiometer. 

The  use  of  the  ADF  in  detection  has  room  for  many  new  developments  in 
just  its  mathematical  development,  not  to  mention  specific  applications  in  de¬ 
tection  theory.  For  instance,  the  basic  results  possibly  could  be  extended  to 
non-Gaussian  and  nonstationary  noise.  In  terms  of  the  deconvolution  statistic, 
the  class  of  kernels  used  was  very  narrow.  The  investigation  of  various  kernels, 
especially  those  used  in  inv>ge  processing,  would  be  in  itself  a  worthwhile  un¬ 
dertaking.  Within  the  class  of  kernels  presented,  the  choice  of  an  optimal  kernel 
for  various  classes  of  detection  problems  would  be  a  possible  topic  to  pursue. 
Finally,  the  most  fertile  ground  for  the  application  of  the  ADF  to  detection 
could  be  the  detection  of  noise  in  noise. 

4. A  PROOF  OF  LEMMA  3 

We  first  consider  the  sets  A,  which  are  finite  half  open  intervals.  Of  these, 
it  is  necessary  to  consider  only  the  interval  [0,T),  since  Y(t)  is  stationary. 
Partition  A  =  [0,T)  into  n  subintervals 

B{  —  (i  —  I)  — ,  i—  j  for  r  =  1, . . . ,  n  (4.62) 

n  n ) 

and  define  the  set  Y~ 1  =  [t  :  Y{l)  <  a].  Observing  that  the  length  of  the  set 
y-1  O  A  is  the  sum  of  the  lengths  of  the  sets  V'a_I  0  B,,  we  can  write 


n 


e[C(Y-lDA)}  =  ££[£(>--' OB,)] 

I  =  1 

(4.63) 

=  ^[qy-'ntf,)] 

(4.64) 

since  V(<)  is  stationary.  Now  define  the  following  three  events: 

C  is  the  event  that  Y{t)  <  a,  for  some  t  G  B i 

D  is  the  event  that  Y(t)  <  a,  for  all  t  E  B\ 

E  is  the  event  that  Y(t)  crosses  a  on  B\. 

Notice  that 

C(Y-'nBt)>  -ID 
n 

(4.65) 
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where  Id  is  the  indicator  of  the  event  D.  To  understand  this  relationship, 
consider  the  case  that  the  sample  path  V(f)  is  in  D,  meaning  that  it  is  not 
above  a  during  the  entire  interval  fl, .  It  follows  that  the  amount  of  time  that  it 
is  not  above  a,  i.e.  C(Y~l  fl  B\ ),  equals  the  length  of  B, ,  which  is  Tfn.  Upon 
taking  expectations  of  this  relationship, 

£  [£(V7‘  D  fl,)]  >  Zeno)  (4.66) 

=  — Pr(  D)  (4.67) 

n 

=  -  {PrfV'(O)  <  a]  -  Pr[fl  and  (V'(0)  <  a)]}  (4.68) 

n 

since  the  probability  that  Y(t)  is  not  above  level  a  over  the  interval  fl,  is  exactly 
the  probability  that  V'(0)  is  not  above  a  and,  under  this  condition,  Y(t)  does 
not  cross  a.  Now,  since  Pr[£  and  (V'(0)  <  a)]  is  less  than  or  equal  to  Pr(£7), 
we  have 

£  [C(Y~l  nfl,)]  >  T-  {Pr[V(0)  <  a]  -  Pr (£)}  .  (4.69) 

In  an  analogous  manner,  we  can  produce  a  complementary  inequality  by 
observing  that 

C(Y~‘  n  Bi)  <  -fc-  (4.70) 

n 

This  inequality  follows  by  considering  two  cases.  When  the  sample  path  Y(t)  is 
in  C,  meaning  that  it  is  below  a  sometime  during  the  interval  B j,  the  amount 
of  time  that  it  is  not  above  a  is  not  greater  than  the  length  of  B ;  or  Tjrx. 
Alternatively,  whenever  Y(t)  is  not  in  C,  it  is  not  above  a  for  zero  time.  Taking 
expectations  of  this  inequality  yields 

£  [£07*  ns,)]  <  ^£(lc)  (4.71) 

=  -Pr(C')  (4.72) 

n 

=  -  {Pr[T(0)  <  a]  -  Pr[£7  and  (V(0)  >  a)]}  (4.73) 

n 

since  the  probability  that  Y{t)  is  below  level  a  for  some  time  during  the  interval 
B i  is  exactly  the  probability  that  y(0)  is  not  below  a  but  crosses  a  during  the 
interval  B\  plus  the  probability  that  V(0)  is  below  a  initially.  Now,  since 
Pr[£  and  V(0)  >  a]  is  less  than  or  equal  to  Prff),  we  have 

£  [C(Y~l  n  fl,)]  <  ^  {Pr[V (0)  <  a]  -  Pr(fl)}  .  (4.74) 

Equations  (4.69)  and  (4.74)  applied  to  (4.64)  imply 

T  (Pr[V'(0)  <  a]  —  Pr(£)}  <  £  [£  (Ta' 1  n  fl, )]  <  T  {Pr[V(0)  <  a]  +  Pr(£)} 

(4.75) 
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thus 


€  [C  (V--1  n  /?!>] 


Pr[V(0)  <  a)T 


(4.76) 


(4.77) 


if  linin^oo  Pr(E)  0. 

In  order  to  prove  that  l;mrl_00  Pr(£)  — *  0,  we  define  the  counting  process 
Na(t)  as  the  number  of  crossings  of  the  threshold  a  by  the  process  Y'(l)  on  the 
interval  [0,<).  By  Chebyshev’s  inequality, 


Pr(£7) 


but  from  Karlin  and  Taylor. 28 


e 


T  a 

- e 


a  2 
2<r0 


n  (T0 


(4.78) 

(4.79) 


(4.80) 


where  a2  =  -  fi( ())".  Upon  letting  n  — *  oc,  the  last  two  equations  imply 
Pr(/T)  —  0.  Now  that  the  result 


£  (avr'n.d)]  =  *(j^jca  (4.81) 

has  been  proven  for  A ,  an  interval,  it  can  be  extended  to  any  finite  set  as  follows. 
Let  A  be  a  set  of  finite  length;  then  for  any  c  >  0  there  exists  a  finite  set  of 
intervals  {/,,}j  =  l  such  that 


and 


■X  ~  U  L' 


i=  1 


( 

<  - 
~  2 


r 


(see  Royden29).  Hence, 

i  i 

Y^C(Y-'  <  C(Y~lnA)  <  ^£(v;lru1)+^. 

i=t  ■=! 


(4.82) 


(4.83) 


(4.84) 


Taking  expectations  and  applying  the  result  for  intervals 


<  €{C(Y~lnA)] 


<  «I> 


CLi  + 


e 

2 


(4.85) 


28S.  Karlin  and  H.M.  Taylor,  A  First  Course  in  Stochastic  Processes ,  Second  Edition,  Aca¬ 
demic  Press,  San  Diego,  Calif.,  1975. 
a9H.L.  Royden,  Real  Analysis,  MacMillan,  New  York,  19f»8. 
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But 

i 

LA~\  -  HCL' s  CA  + 1 

i=l 

implying 


$ 


CA  —  ( 


<  £[C(Y~lCiA) j  <  <J>  (  — 'j  CA  +  r. 

V^o  / 


(4.86) 


(4  87) 


But,  since  c  is  an  arbitrary  positive  number,  we  have  for  a  set  A  with  finite 
length, 


€[C(Y-'nA)\  =  *(±'jcA.  (4.88) 

To  extend  to  the  case  in  which  ,4  is  not  of  finite  length,  write  A  =  U“0.4, , 
where  the  /L's  are  disjoint  and  of  finite  length;  then 

OO 

s  [£(>Vl  n  A)]  =  ^^  [£(>;-' nT,)]  (4.89) 

i=0 

OO  /  \ 

=  <«»> 

=  ,'-9” 


4.B  PROOF  OF  THEOREM  3 


Without  loss  of  generality,  assume  Fs  is  continuous.  To  prove  the  result, 
we  wish  to  compute 


Fy(a)  =  Yim^  ~  {£  [Y~l  n  [0,7’)] }  (4.92) 

where 

V-1  =(<:  S(t)+ N(l)<a].  (4.93) 

Begin  by  selecting  an  integer  m  and  constructing  ?  partition  of  the  real  line. 
(s_ ma  ^  ^  ^mJ)>  where 


—  OO 

for  «  =  —  rn2 

i/m 

for  —  m2  <  i  <  rn 

OO 

for  i  -  m2. 

Continue  by  defining  the  sets 


(4.94) 


Ai  =  [t  :  Pu(t)  <a  —  Si  and  s*_i  <  S(t)  <  s*]  (4.95) 

A?  =  (<  :  N(t)  <  a  -  s,_i  and  s,_i  <  5(1)  <  s,]  (4.96) 

for  i  =  —  m2  +  1,  m2.  Observe  that  A~  C  Y„~'<  because  S(l)  +  JV(1)  <  a 

whenever  1V(1)  <  a  -  s,  and  <  5(1)  <  s,  Hence, 


UA'cl;1.  (4.97) 

»=1 
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Furthermore, 

n 

Y-lc{jAf  (4.98) 

i=l 

because,  for  any  t  where  N(t)  +  5(<)  <  a,  there  exists  an  i  in  the  range 
— m2  <  i  <  m2,  such  that  «i_i  <  S(t)  <  s,  implying  N(t)  <  a  —  s<.  Be¬ 
cause  {A"}”3_ma+1  and  {At }^_m3+1  are  disjoint,  (4.97)  and  (4.98)  imply 
that 

1  m  i 

f€{C[Y~ln[0,T)]}  >  £  j^{C[A~n[0,T)}}  (4.99) 

»  =  —  m3  + 1 

■  m3  i 

^{rfy-^io.r)]}  <  £  f£  ^  E^,+  n  [0,T)] (4100) 

»=:  —  m3  +  1 

Now  notice  that 

lim  {C  [F"1  n[0,T)]}  ~  Fy(a).  (4.101) 

T—  oo  1 

Since  S(t)  and  ,V(<)  are  amplitude  independent,  Fs  is  continuous,  and  Fs(+ oo) 
=  lim^oo  Fs(s)  and  Fs(-oo)  =  lim,—.,*  Fs(s),  it  follows  that 

lim  U{C[A-n[Q,T)}}  =  FN(a-s,)[Fs(si)-Fs(si-l)].  (4.102) 

T—ooJ 

Similarly, 

{£[A,+  n(0,T)]}  =FA(a-si.,)[Fs(s,)-Fs(s<_1)].  (4.103) 

By  defining  two  particular  step  functions,  the  above  results  with  equation 
(4.100)  and  (4.99)  form  a  relationship  between  the  ADF  of  S(t)  +  N(t)  and 
the  integrals  of  the  two  step  functions.  The  step  functions  are 

Ffr  ( a  —  s )  =  FN(a-Si)  (4.104) 


F+{a~s)  =  Fjv(« -*<_,)  (4.105) 

whenever  s,_!  <  s  <  s,.  Upon  passing  T  to  oo,  (4.100)  and  (4.99)  become  with 
the  aid  of  the  above  definitions  and  (4.101).  (4.102),  and  (4.103), 


/  Fy(a-s)dFs(s) 

J  —  OG 

[OO 

<  Fy(a)  <  /  F+(a-s)dFs(s). 

J  — OO 

(4.106) 

Proceed  by  enlarging  m 
which  coupled  with  (4.106), 

and  find  that  FN  and  Fy  converge  weakly  to 
implies  that 

Fy(a) 

=  f  Fs(a  -  s)  dFs(s) 

J  -  OG 

(4.107) 

=  f  FN(a  -  s)  dF's(s) 

J  — OO 

(4.108) 
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where  the  last  equation,  following  from  the  continuity  of  F$,  proves  the  first 
convolution.  Note  that  the  interchange  between  limit  and  integration  is  justified 
in  the  last  operation,  because  the  integrals  have  essentially  all  mass  within  a 
bounded  domain,  upon  which  the  integrand  is  bounded  between  zero  and  one. 
Note  also  that,  if  F n  is  not  continuous,  then  Ffi  or  Fj y  may  converge  to  some 
function  that  differs  from  F/v  at  a  countable  number  of  points,  but  luckily  the 
hypothesis  that  F5  is  continuous  makes  (4.106)  invariant  to  the  limit  value  of 
Fn  and  Fjv  at  these  problem  points. 

The  complementary  convolution  is  obtained  by  integrating  the  product  mea¬ 
sure  dFg(s)  x  dFtf(n)  over  the  half  plane  //  =  [s,  n  :  s  +  n  <  a}.  Proceeding, 
we  obtain  with  the  help  of  Fubini’s  theorem 


/  dF's(s)  x  dF-N(n)  =  r\r  * dF'N(n ) 

J H  oo  U-oq 

]  dF's(s) 

(4.109) 

=  f  F'N(a  -  s)dF's( 

J  —  OO 

s ) 

(4.110) 

or,  alternatively, 

/  dF-s(s)  x  dF'N(n)  =  ^  f  f  "  dFg(s) 

J H  J — oo  LJ-oo 

drN(n) 

(4.111) 

=  j00  Fs(a-n)dF'N( 

J  —  OO 

n) 

(4.112) 

implying 

r  m*-»)dm»)=  r  Fs(a-n)dF^n). 

(4.113) 

—  oc  J  —  OO 

The  continuity  of  Fs  and  (4.108)  yields 

Fy(a)=  ^  Ftf(a  —  s)  dFg(s) 

(4.114) 

which  with  (4113) 

Fv(a)  =  f  Fs(a-n)dF-„(n) 

J  —  OO 

(4.115) 

=  f  Fs(o  -  n)  dFff(n) 

(4.116) 

-OO 


by  the  continuity  of  Fs .  Equation  (4.116),  the  remaining  convolutional  rela¬ 
tionship,  is  now  proven. 

4.C  PROOF  OF  THEOREM  4 


Begin  by  selecting  an  integer  rn  and  constructing  a  partition  of  the  real  line 


(s_mi  <  s_mj+i  <  •  •  •  <  »mj),  where 


for  i  =  —m2 

for  — m2  <  i  <  m2  (4.117) 

for  i  =  m2. 
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Continue  by  defining  the  sets 

=  [<  :  Si_i  <  5(0  <  «i] 


(4.118) 


for  i  =  -m2  +  l,  -,m2.  The  sets  {Si}  partition  the  real  line  into  disjoint 
subsets,  therefore, 

ma  . 

/  FN[a-S(t)]dt  =  V  /  FN[a-S(t)]dt.  (4.119) 

J(0,T]  te_»*+i  ■'S.n(°.Tl 


Furthermore,  since  F/v  is  increasing, 


m3 

f  FN[a-S(t)]dt  >  T  I  FN(a-8i)dt  (4. 

J(0,T]  iJ=M0.T) 


120) 


and 


.  m3  . 

/  FN[a-  S(t)]dt  <  Y]  /  Fs(a 

J(0,T)  i  =  -m^  +  l'/H'n(°'Tl 


-Si_i)dt.  (4.121) 


Divide  (4.120)  by  T,  take  expectations,  and  enlarge  T ;  then,  noting  that  Fs(+oo) 
=  lim,_oo  Fs(s)  and  Fs(-oc)  =  lim,  —  Fs(s),  implies 


lim 

T-oo 


Hi 


F*[a  -  5(0]  dt 

Operate  similarly  on  (4.121)  and 


1  m3 

}>-  E 

)  i  =  -m3  +  l 


FN(a  -  s,)  [Fs(si)  -  Fs(si_j)]  • 

(4.122) 


lim  \-C  {  /  Fat[o  -  5(0]  dt  1  <  Y  F/v(a-Si_i)  [Fs(si)  -  Fs(si_i)] . 
T~°°  1  [Jo  )  i=-m3+[ 

(4.123) 

In  order  to  express  (4.122)  and  (4.123)  in  convolutional  form,  define  the  two 
step  functions 


FN(a-s)  -  FN(a-Si) 
F^(a-s)  =  Fn(q  —  Si_i) 


(4.124) 

(4.125) 


whenever  s;_  i  <  s  <  s,.  Since  Fs  is  assumed  to  be  at  least  right-continuous,  it 
follows  from  (4.122)  and  (4.123)  that 

f  Ftf(a  -  s)  dF|(s)  <  ^5  jjT  FN[a  -  5(0]  d<  j  (4.126) 


and  hence 


/OO  FOO 

F^(a-s)dF's(s)  <  /  Ffi (a  —  s)  dFg(s).  (4.127) 

•OO  j  —  OO 


115 


The  Johns  Hopkins  University 

Applied  Physics  Laboratory 


Finish  the  proof  by  enlarging  rn  and  find  that  and  Ffi  converge  weakly 
to  Fjv  which,  coupled  with  (4.127),  implies 

FN[a~S{t)]  =  /  F"(a  ~ S)  dF's(s)  (4128) 

=  Fy(a)  (4.129) 

where  the  final  step  follows  from  Theorem  3.  Note  that  the  interchange  between 
limit  and  integration  is  justified  in  the  last  operation,  because  the  integrals 
have  essentially  all  mass  within  a  bounded  domain  upon  which  the  integrand  is 
bounded  between  zero  and  one.  Note  also  that,  if  F^ r  is  not  continuous,  then 
Fj y  or  Fjj  may  converge  to  some  function  that  differs  from  Fn  at  a  countable 
number  of  points;  but  for  this  case,  F?  is  hypothesized  to  be  continuous,  making 
(4.128)  invariant  to  the  limit  value  of  and  Fjy  at  these  problem  points. 
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CHAPTER  5 

CONCLUDING  REMARKS  AND 
EXTENSIONS 

5.1  SUMMARY  AND  REMARKS 

In  the  first  chapter,  the  general  setting  for  the  interception  problem  was  de¬ 
scribed.  A  brief  exposition  on  the  function  of  the  intercept  receiver  was  given, 
emphasizing  the  functions  of  initial  detection  and  feature  detection  that  are  per¬ 
tinent  to  this  work.  A  brief  description  of  existing  detectors  was  presented  with 
a  qualitative,  and  in  some  cases  quantitative,  performance  evaluation.  These 
were  the  radiometer,  the  optimal  channelized  detector  with  various  suboptimal 
simplifications,  and  the  autocorrelation  detector.  Of  existing  feature  detectors 
of  hop  frequency,  a  maximum  likelihood  receiver  was  described  along  with  a 
related  receiver  employing  course  subband  selection.  The  detectors  presented 
in  this  work  were  then  briefly  described  and  related  to  existing  detectors.  They 
were  a  sequential  detector,  two  detectors  and  a  hop  frequency  estimator  based 
on  the  compressive  receiver,  and  another  detector  based  on  the  new  concept  of 
the  amplitude  distribution  function. 

In  the  second  chapter,  methods  for  the  sequential  detection  of  noncoher¬ 
ent  fast  FII  waveforms  were  developed.  In  the  process,  the  FH  waveform 
was  modeled  to  have  an  information  component,  which  consisted  of  a  series 
of  chips  with  a  known  constant  epoch  where  each  chip  frequency  was  one  of 
a  known  ensemble  of  frequencies.  In  the  model,  a  particular  chip  frequency 
was  independently  determined  by  a  uniform  random  variable  on  the  frequency 
ensemble.  The  FH  waveform  also  was  assumed  to  have  an  additive  white-noise 
component.  By  assuming  the  modeled  FII  waveform  was  of  a  known  SNR,  the 
optimal  detector  based  on  a  single-epoch  observation  (SELF)  was  developed 
using  likelihood-function  theory,  SELF  was  the  sum  of  many  nearly  identical 
and  nearly  independent  random  variables  and  thus  had  nearly  Gaussian  statis¬ 
tics.  This  central-limit  argument  allowed  a  multi-epoch  collection  of  SELFs  to 
be  considered  an  equivalent  set  of  Gaussian  i.i.d.  variables.  From  these  sim¬ 
plified  observations,  a  log-likelihood  function  (ALLF)  was  computed  that  was 
asymptotic  to  the  exact  log-likelihood  function,  as  the  number  of  possible  hop 
frequencies  became  large.  The  ALLF  became  the  test  statistic  on  which  three 
detection  tests  were  based.  The  tests  were  the  FSS  test,  the  SPRT,  and  the 
TST.  These  were  defined  to  ensure  that  detection  errors  were  below  desired 
levels.  By  modeling  the  ALLF  as  a  Wiener  process,  diffusion  theory  yielded 
the  performance  of  the  three  tests  not  only  for  an  FH  waveform  of  the  assumed 
SNR,  but  also  the  test  performance  for  all  SNRs  below  the  one  assumed.  This 
analysis  compared  favorably  with  a  computer  simulation  of  the  detector  and 
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thus  validated  the  analysis.  The  analysis  also  became  a  tool  used  to  optimize 
numerically  the  performance  of  the  TST  when  the  actual  FH  SNR  deviated 
from  the  one  assumed.  In  order  to  study  the  performance  of  tests  synthesized 
by  assuming  an  extremely  small  FH  SNR,  expressions  for  the  asymptotic  test 
efficiencies  were  computed.  This  asymptotic  analysis  also  yielded  simplified  test 
parameter  expressions  applicable  to  the  small-SNR  case. 

A  significant  feature  of  the  SPRT  exposed  by  the  analysis  is  that,  with  the 
same  error  probabilities,  an  FH  waveform  with  a  given  SNR  can  be  detected 
in  less  than  half  the  time  of  the  corresponding  FSS  test.  This  reduction  in 
detection  time  is  especially  significant  for  Low  Probability  of  Intercept  (LP1) 
applications,  where  the  transmissions  are  purposely  short.  For  the  pure  SPRT, 
detection  time  increased  whenever  the  observed  SNR  differed  from  that  assumed 
in  the  test’s  synthesis.  And  for  SNRs  midway  between  zero  and  the  assumed 
value,  it  was  even  comparable  to  the  corresponding  FSS  test.  The  TST  signif¬ 
icantly  improves  this  anomaly,  while  sacrificing  little  performance  over  that  of 
the  purely  sequential  test;  whatever  little  performance  is  lost,  the  optimal  TST 
largely  regains.  The  decrease  in  the  detection  time  of  the  sequential  tests  can 
be  used  to  robustify  the  test  with  respect  to  the  input  SNR  while  maintaining 
better  performance  than  that  of  the  non-robust  FSS  test.  The  simplified  test 
parameter  expressions  derived  by  asymptotic  methods  may  be  useful  for  any 
schemes  to  adapt  these  tests  for  varying  FH  SNRs.  The  three  tests  and  their 
corresponding  design  and  performance  analysis  also  apply  to  the  slow-FH  case. 
The  detector  structure  is  suboptimal  for  slow  FH,  but  it  is  believed  that  the 
performance  loss  is  small,  especially  for  detection  times  that  include  a  large 
number  of  hops. 

In  the  third  chapter,  two  detectors  of  frequency-hopped  waveforms  based  on 
the  compressive  receiver  were  presented.  The  first  was  developed  by  applying 
likelihood  ratio  theory  to  the  observed  compressive-receiver  output  and  yielded 
a  locally  optimal  (low-SNR)  detector.  The  second,  motivated  by  simplicity  of 
implementation,  was  a  time-multiplexed  version  of  the  first  that,  through  the 
choice  of  a  parameter,  could  either,  at  the  expense  of  a  low  duty  cycle,  achieve 
the  detectability  of  the  first  or  could,  at  the  expense  of  degraded  performance, 
achieve  higher  duty  cycles.  The  second  detector  was  modified  into  a  maximum 
likelihood  estimator  of  hop  frequency.  Both  detectors  and  the  hop  frequency 
estimator  were  performance  analyzed  and  compared. 

The  compressive  receiver  fulfilled  its  promise  as  a  simple,  yet  high-performing 
interceptor.  The  performance  of  the  locally  optimal  detector  shows  that  rela¬ 
tively  little  detectability  is  lost  by  the  processing  of  the  compressive  receiver. 
Most  of  the  discrepancy  is  due  to  the  difference  in  coherent  integration  time 
(one  half  for  the  parameters  used).  Furthermore,  for  a  small  performance 
cost,  the  simplicity  of  the  compressive-receiver  approach  can  be  retained  by 
the  time-multiplexed  detector.  Also,  the  hop  frequency  estimator  again  com¬ 
pares  favorably  with  the  corresponding  device  that  used  raw  input  instead  of 
compressive-receiver  output. 

In  the  fourth  chapter,  a  new  idea  in  detection,  the  Amplitude  Distribution 
Function  (ADF),  was  introduced.  The  ADF  is  roughly  the  average  probability 
distribution  of  a  random  process.  Because  the  ADF  is  original,  a  mathematical 
foundation  was  laid  consisting  of  a  sequence  of  definitions,  lemmas,  and  theo¬ 
rems,  the  most  significant  of  which  was  the  fact  that  the  ADF  of  signal  plus 
noise  is  the  convolution  of  the  ADF  of  signal  and  the  ADF  of  noise  taken  sepa¬ 
rately.  The  technique  of  deconvolution  used  in  image  processing  was  the  germ 
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from  which  emerged  the  Deconvolution  Statistic,  a  statistical  transform  that 
mapped  a  stochastic  process  into  the  amplitude  domain.  It  was  shown  that 
by  proper  choice  of  a  kernel  any  degree  of  separation  between  the  noise  and 
signal  components  could  be  obtained.  For  the  particular  problem  of  detecting 
a  modulated  sinusoid  in  stationary  Gaussian  noise,  a  detector  was  developed 
around  the  Detection  Statistic.  The  detector’s  performance  was  analyzed  and 
compared  with  a  radiometer. 


5.2  EXTENSIONS — DIRECTIONS  FOR 
FUTURE  WORK 


In  applying  sequential  detection  to  interception,  other  simplifications  and 
extensions  are  possible.  For  instance,  it  was  assumed  that  the  starting  time  and 
duration  of  the  chip  epoch  were  known.  This  first  restriction  might  be  relaxed 
by  redefining  the  SELF  to  perform  sliding  window  integration  instead  of  the 
integrate-and-dump  operation  now  performed.  This,  of  course,  would  degrade 
the  detector’s  performance  for  some  values  of  epoch  starting  time,  but  it  would 
probably  exhibit  a  better  average  performance.  There  are  also  possible  simpli¬ 
fications  to  the  SELF  to  improve  its  implementability.  Among  these  could  be 
the  removal  of  the  emphasizing  function,  which  would  make  the  detector  struc¬ 
ture  suboptimal  but  it  probably  still  would  be  asymptotically  optimal  for  small 
assumed  SNRs.  Another  simplification  could  be  coarse  subband  preselection, 
where  the  total  spread-spectrum  bandwidth  is  subdivided  into  subbands,  each 
containing  a  large  number  of  chip  frequencies.  An  algorithm  could  be  used  to 
select  a  subset  of  the  subbands  most  likely  to  contain  the  intercepted  signal. 
Detailed  processing  on  these  preselected  bands  then  could  be  done  with  the 
methods  described  in  this  work. 

There  are  many  remaining  avenues  to  be  traveled  in  the  use  of  the  compres¬ 
sive  receiver  to  interception.  We  provided  many  results  applying  to  frequency- 
hopped  waveforms,  but  the  essence  of  these  ideas  can  apply  to  other  spread- 
spectrum  modulations  as  well.  But  even  in  realm  of  frequency-hopped  wave¬ 
forms,  much  work  remains  to  be  done.  For  instance,  a  locally  optimal  detector 
was  derived  consisting  of  a  bank  of  filters  operating  on  the  compressive-receiver 
output.  From  the  filter  responses  in  this  configuration  was  formed  the  filter 
response  for  the  time-multiplexed  detector.  The  direct  derivation  of  this  fil¬ 
ter  response  using  some  optimality  criterion  over  the  class  of  expected  signals 
might  be  a  promising  endeavor.  The  extension  of  the  detectors  to  a  multihop 
observation  period  was  largely  ignored.  We  simply  assumed  that  results  of  the 
individual  detectors  would  be  combined,  as  if  they  were  independent.  Here 
the  issue  of  performance  versus  overlap  between  data  windows  is  one  to  be  po¬ 
tentially  explored.  We  showed  how  the  compressive  receiver  could  be  used  to 
estimate  hop  frequency.  Other  types  of  feature  detectors  such  as  the  hop-rate 
detectors  or  the  carrier-frequency  detectors,  also  could  be  pursued. 

The  use  of  the  ADF  in  detection,  being  a  new  idea,  has  room  for  many 
new  developments  in  just  its  mathematical  development,  not  to  mention  spe¬ 
cific  applications  in  detection  theory.  For  instance,  the  basic  results  could  be 
extended  to  non-Gaussian  and  nonstationary  noise.  In  terms  of  the  Deconvo¬ 
lution  Statistic,  the  class  of  kernels  used  was  very  narrow.  The  investigation  of 
various  kernels,  especially  those  used  in  image  processing,  would  be,  by  itself, 
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a  worthwhile  undertaking.  Within  the  class  of  kernels  presented,  the  choice  of 
an  optimal  kernel  for  various  classes  of  detection  problems  is  a  possible  area 
to  pursue.  Finally,  the  most  fertile  ground  for  the  application  of  the  ADF  to 
detection  could  be  the  detection  of  noise  in  noise. 


5.3  FINAL  REMARKS 


This  work  represents  a  significant  step  in  the  development  of  new  detec¬ 
tion  and  interception  techniques  for  spread-spectrum  waveforms  and  frequency- 
hopped  waveforms  in  particular.  With  the  three  new  detectors  and  the  hop 
frequency  estimator  presented,  many  new  trails  were  blazed.  Armed  with  the 
idea  of  sequential  detection,  an  existing  optimal  detector  using  a  fixed  number 
of  samples  was  extended,  with  improved  performance  resulting.  This  useful  and 
eclectic  endeavor  brought  previous  theoretical  and  practical  results  on  sequen¬ 
tial  detection  into  the  context  of  interception.  The  second  detector  and  hop 
frequency  estimator  solved  the  problem  of  how  to  apply  effectively  the  com¬ 
pressive  receiver  to  interception.  The  firm  mathematical  development  starkly 
contrasted  with  previous  ad  hoc  attacks  on  the  problem.  Out  of  this  work  also 
came  some  mathematical  results  of  general  interest;  among  these  were  bound¬ 
ing  distributions  on  the  sum  of  squares  of  Gaussian  random  variables  and  an 
extension  of  the  Riemann-Lebesgue  Lemma  to  integrals  of  linearly  frequency- 
modulated  sinusoids.  A  brand  new  idea  in  detection  was  conceived,  yielding  yet 
another  new  detector.  The  original  idea  melded  the  image  processing  technique 
of  deconvolution  to  those  of  density  estimation.  The  detector  thus  developed 
was  quite  independent  of  the  details  of  signal  modulation. 

While  there  is  much  presented  here,  there  is  still  much  to  be  done.  There  are 
other  interceptors  to  be  analyzed  and  new  detection  techniques  to  be  developed. 
As  such,  the  area  of  spread-spectrum  interception  will  yield  new  results  for  many 
years  to  come. 


120 


The  Johns  Hopkins  University 

Applied  Physics  Laboratory 


REFERENCES 


*M.K.  Simon,  J.K.  Omura,  R.A.  Scholtz,  and  B.K.  Levitt,  Spread  Spectrum 
Communications,  voi.  Ill,  Computer  Sciences  Press,  Rockville,  Md.,  1985. 

2D.  Torrieri,  Principles  of  Secure  Communication  Systems.  Artech  House, 
Dedham,  Mass.,  1985. 

3D.G.  Woodring,  “Performance  of  Optimum  and  Suboptimum  Detectors  for 
Spread  Spectrum  Waveforms” ,  Naval  Research  Laboratory,  Washington,  D.C., 
Technical  Report  No.  8432,  December  1980. 

4A.  Polydoros  and  J.K.  Holmes,  “Autocorrelation  Techniques  for  Wideband 
Detection  of  FH/DS  Waveforms  in  Random  Tone  Interference”,  MILCOM 
’83  Conference  Proceedings,  Boston,  Mass.,  pp.  781-785,  October  17-20, 
1983. 

5A.  Polydoros  and  K.T.  Woo,  “LPI  Detection  of  Frequency-Hopping  Signals 
using  Autocorrelation  Techniques”,  IEEE  J.  Select.  Areas  Commun.,  vol. 
SAC-3,  no.  5,  September  1985. 

6N.C.  Beaulieu,  W.L.  Hopkins,  and  P.J.  McLane,  “Interception  of  Frequency 
Hopped  Spread  Spectrum  Signals”,  to  be  published  in  IEEE  Trans.  Commun 

7W.W.  Short  and  R.D.  Chapman,  “Adaptively  Configured  Channelized  Re¬ 
ceiver  for  Frequency  Hopped  Signal  Detection  and  Tracking”,  IEEE  Interna¬ 
tional  Conference  on  Communications,  1985,  pp.  832-838,  1985. 

8J.D.  Edell,  “Wideband,  Noncoherent,  Frequency-Hopped  Waveforms  and 
their  Hybrids  in  Low-Probability-of-Intercept  Communications”,  Naval  Re¬ 
search  Laboratory,  Washington,  D.C.,  Technical  Report  No.  8025,  November 
1976. 

9A.  Wald,  Sequential  Analysis,  Wiley,  New  York,  1947. 

10N.  Shiryayev,  Optimal  Stopping  Rules,  Springer- Verlag,  New  York,  1977. 

nD.A.  Darling  and  A.J.F.  Siegert,  “The  First  Passage  Problem  for  a  Contin¬ 
uous  Markov  Process”,  Ann.  Math.  Stat.,  vol.  24,  pp.  624-639,  1953. 

12T.W.  Anderson,  “A  Modification  of  the  Sequential  Probability  Ratio  Test  to 
Reduce  the  Sample  Size”,  Ann.  Math.  Stat.,  vol.  31,  pp.  165-197,  1960. 

13S.  Tantaratana  and  H.V.  Poor,  “Asymptotic  Efficiencies  of  Truncated  Sequen¬ 
tial  Tests”,  IEEE  Trans.  Inform.  Theory,  vol.  IT-28,  no.  6,  pp.  911-923, 
November  1982. 

14S.  Tantaratana  and  J.B.  Thomas,  “Truncated  Sequential  Probability  Ratio 
Test”,  Inform.  Set.,  vol.  13,  pp.  283-300,  1977. 

15G.N.  Watson,  A  Treatise  on  the  Theory  of  Bessel  Functions,  Cambridge  Uni¬ 
versity  Press,  New  York,  1980. 

l6A.  Polydoros,  J.K.  Holmes  and  K.T.  Woo,  “Advanced  LPI  Intercept  Detector 
Research”,  Technical  Report  No.  R8511-3,  Axiomatix,  Los  Angeles,  Calif., 
November  13,  1985. 

17W.E.  Snelling  and  E.A.  Geraniotis,  “Sequential  Detection  of  Unknown,  Fast 
Frequency- Hopped  Waveforms”,  IEEE  J.  Select.  Areas  Commun.,  vol.  7,  no. 
4,  pp.  602  617,  May  1989. 

18K.D.  Breuer,  J.J  Whelehan,  and  K.  Ross,  “Compressive  Receivers  Applied 
to  ESM  System  Design”,  MSN  Microwave  Syst.  News  &  Commun.  Techno., 
vol.  16,  no.  11,  pp.  66-68,  70,  72,  74-75,  Oct.  1986. 


121 


The  Johns  Hopkins  University 

Applied  Physics  Laboratory 


19B.K.  Harms  and  D.R.  Hummels,  “Calculation  of  Detection  Probability  for 
Frequency  Compressive  Receivers”,  IEEE  Trans.  Aerospace  and  Electronic 
Systems,  vol.  AES-21,  pp.  106-116,  January  1985. 

20K.H.  Li  and  L.B.  Milstein,  “On  the  Use  of  a  Compressive  Receiver  for  Signal 
Detection”,  to  be  published  in  IEEE  Trans.  Commun. 

2lH.L.  van  TVees,  Detection,  Estimation,  and  Modulation  Theory,  Part  I,  Wi¬ 
ley,  New  York,  1968. 

22A.  Papoulis,  Probability,  Random  Variables,  and  Stochastic  Processes,  Second 
Edition,  McGraw-Hill,  New  York,  1984. 

23R.F.  Pawula,  S.O.  Rice,  and  J.H.  Roberts,  “Distribution  of  the  Phase  Angle 
Between  Two  Vectors  Perturbed  by  Gaussian  Noise”,  IEEE  Trans.  Com¬ 
mun.,  vol.  COM-30,  pp.  1828-1840,  August  1982. 

24J.H.  Higbie,  “Adaptive  Nonlinear  Suppression  of  Interference”,  MILCOM  '88 
Conference  Proceedings,  23.3.1-23.3.9,  October  1988. 

25S.O.  Rice,  “Statistical  Properties  of  a  Sine  Wave  Plus  Random  Noise,”  Bell 
Syst.  Tech.  J.,  pp.  109-157,  January  1948. 

26J.  Salz,  and  S.  Stein,  “Distribution  of  Instantaneous  Frequency  for  Signal 
Plus  Noise”,  IEEE  Trans.  Inform.  Theory,  vol.  IT-10,  pp.  272-274,  October 
1964. 

27H.C.  Andrews  and  B.R.  Hunt,  Digital  Image  Restoration,  Prentice-Hall,  En¬ 
glewood  Cliffs,  New  Jersey,  1977. 

2SS.  Karlin  and  H.M.  Taylor,  A  First  Course  in  Stochastic  Processes,  Second 
Edition,  Academic  Press,  San  Diego,  Calif.,  1975. 

29H.L.  Royden,  Real  Analysis,  MacMillan,  New  York,  1968. 


122 


The  Johns  Hopkins  University 

Applied  Physics  Laboratory 


BIBLIOGRAPHY 


'T.W.  Bristol,  “Review  of  Spectrum  Analysis  with  SAW  Chirp  Transforms 
and  Filter  Banks”,  Proceedings  of  International  Specialist  Seminar  on  Case 
Studies  in  Advanced  Signal  Processing,  IEE,  London,  Engl.,  pp.  226-231, 
Sept.  1979 

2P.  Das  and  L.B.  Milstein,  “Adaptive  Spread  Spectrum  Receiver  Using  Acous¬ 
tic  Surface  Wave  Technology”,  Report  No.  ARO  17570-8-EL,  May  31  1984, 
U.S.  Army  Research  Office,  Post  Office  Box  12211,  Research  Triangle  Park, 
N.C.  27709. 

3G.M.  Dillard,  “A  Moving-Window  Detector  for  Binary  Integration”,  IEEE 
Trans.  Inform.  Theory,  IT-13,  no.  1,  pp.  2-6,  January  1967. 

4R.A.  Dillard,  “Detectability  of  Spread  Spectrum  Signals”,  IEEE  Trans.  Aero¬ 
space  and  Electronic  Systems,  vol.  AES-15,  no.  5,  pp.  526-537,  July  1979. 

5R.  Fry  and  P.  Yansouni,  “Detection  and  Direction  Finding  of  Frequency 
Hoppers  by  Amplitude  Comparison  Techniques”,  Defense  Research  Estab¬ 
lishment,  Ottawa,  Ontario  Canada,  Report  No.  940,  December  1986. 

6W.A.  Gardner,  “Signal  Interception:  A  Unifying  Theoretical  Framework  for 
Feature  Detection”,  IEEE  Trans.  Commun.,  vol.  COM-36,  pp.  897-906, 
1988. 

7A.B.  Glenn,  “Low  Probability  of  Intercept”,  IEEE  Commun.  Mag.,  vol.  21, 
no.  4,  pp.  26-33,  July  1983. 

8L.  Kleinrock,  “Detection  of  Energy  Peak  of  an  Arbitrary  Signal”,  MIT  Lincoln 
Laboratory,  Lexington,  Mass.,  Technical  Report  No.  325,  August  23,  1963. 

9N.F.  Krasner,  “Optimal  Detection  of  Digitally  Modulated  Signals”,  IEEE 
Trans.  Commun.,  vol.  COM-30,  no.  5,  pp.  885-895,  May  1982. 

10N.F.  Krasner,  “Maximum  Likelihood  Parameter  Estimation  for  LPI  Signals”, 
MILCOM  '82  Conference  Proceedings,  Boston,  Mass.,  pp.  2. 3-1-2. 3-4,  Octo¬ 
ber  17-20,  1982. 

nR.J.  Mammone,  “Frequency  Resolution  Enhancement  of  a  Compressive  Re¬ 
ceiver  by  Spectral  Estimation”,  Proceedings  of  the  IEEE  Military  Communi¬ 
cations  Conference,  vol.  3,  pp.  713-19,  1983. 

12G.L.  Moule,  “SAW  Compressive  Receivers  for  Radar  Intercept”,  IEE  Pro¬ 
ceedings,  Part  F  —  Communications,  Radar  and  Signal  Processing,  vol.  129, 
pt.  F,  no.  3,  pp.  180-186,  June  1982. 

13D.L.  Nicholson,  “Spread  Spectrum  Signal  Design” ,  Computer  Science,  Rockville, 
Md.,  1988. 

l4J.E.  Ohlson,  “Efficiency  of  Radiometers  Using  Digital  Integration”,  Radio 
Science,  vol.  6,  no.  3,  pp.  341-345,  March  1971. 

I5A.  Polydoros  and  C.L.  Weber,  “Optimal  Detection  Considerations  for  Low 
Probability  of  Intercept”,  MILCOM  '82  Conference  Proceedings,  Boston, 
Mass.,  pp.  2. 1-1-2. 1  5,  October  17-20,  1982. 

I6A.  Polydoros,  J.K.  Holmes,  and  K.T.  Woo,  “Advanced  LPI  Intercept  Detector 
Research”,  Technical  Report  No.  R8403-1,  Axiomatix,  Los  Angeles,  Calif., 
March  26,  1984. 

17A.  Polydoros  and  C.L.  Nikias,  “Detection  of  Unknown  Frequency  Sinusoids  in 
Noise:  Spectral  versus  Correlation  Domain”,  IEEE  Transactions  on  Acous¬ 
tics,  Speech,  and  Signal  Processing ,  vol.  ASSP-35,  no.  6,  June  1987. 


123 


The  Johns  Hopkins  University 

Applied  Physics  Laboratory 


18D.C.  Schleher,  Introduction  to  Electronic  Warfare ,  Artech  House,  Inc.,  Ded¬ 
ham,  Mass.,  1986. 

19W.E.  Sneiiing  and  E.A.  Geraniotis,  “The  Optimal  Interception  of  Frequency- 
Hopped  Waveforms  via  a  Compressive  Receiver”,  to  be  published  in  IEEE 
Trans.  Commun.. 


124 


INITIAL  DISTRIBUTION  EXTERNAL  TO  THE  APPLIED  PHYSICS  LABORATORY* 


The  work  reported  in  TG  1378  was  done  under  Navy  Contract  N00039-89-C-5301  and  is  related  to  Task  BKBR9SLX  F2B,  supported  by  NAVSEA 


ORGANIZATION 

LOCATION 

ATTENTION 

No.  o< 
Copies 

DEPARTMENT  OF  DEFENSE 

— 

Office  of  the  Under  Secretary  of  Defense, 

Washington.  DC  20301 

Accessions 

Research  and  Engineering 

Defense  Technical  Information  Center 

Alexandria.  VA  23314 

Accessions 

12 

Defense  Advanced  Research  Projects  Agency 

Arlington,  VA  20331 

Accessions 

DEPARTMENT  OF  THE  NAVY 

Chief  of  Naval  Operations 

Washington,  DC  20350 

Accessions 

2 

COMOPTEVFOR 

Norfolk.  VA  2351  1-6388 

Library 

1 

Naval  Air  Development  Center 

Warminster,  PA  18974-5000 

Library 

1 

R.  J.  Klump 

1 

A.  C.  Knobloch 

1 

FCDSSA,  Dam  Neck 

Virginia  Beach,  VA  23461-5300 

Library 

1 

Office  of  the  Assistant  Secretary 

of  the  Navy 

Washington.  DC  20350 

Accessions 

2 

Office  of  Naval  Research 

Arlington.  VA  22217 

Library 

2 

K.  G  Dial 

1 

Office  of  the  Naval  Technology 

Arlington.  VA  22209 

Library 

2 

Naval  Air  Systems  Command 

Washington,  DC  22202 

Library.  Air  7226 

2 

R.  H  Ekslrom 

1 

Naval  Electronic  Systems  Command 

Arlington.  VA  20360 

Library 

2 

Naval  Postgraduate  School 

Monterey.  CA  92940 

Library 

2 

Naval  Research  Laboratory 

Washington.  DC  20375 

Library 

2 

J .  W.  Crider 

1 

J.  R.  Davis 

1 

D.  !.  Himes 

1 

Space  &  Naval  Warfare  Systems  Command 

Washington.  DC  20360 

Library 

1 

Nava!  Sea  Systems  Command 

Washington,  DC  22202 

Library,  Sea  9961 

2 

PM5400 

2 

PM5400B 

2 

PM5400B5 

1 

PM5400B51 

1 

PM5400B52 

1 

SEA  06 

1 

SEA  06AP 

1 

J.  F.  Blincoe 

1 

R.  M.  Cassidy 

1 

A  L.  Cetel 

1 

H.  W.  Howard 

1 

R.  G.  Meinig 

1 

R  E.  Muir 

1 

J,  S.  Polk 

1 

A  Rutherford 

1 

J  R.  Whalen 

1 

Naval  Ship  Weapon  Systems 

Engineering  Station 

Fort  Hueneme,  CA  93043-5007 

Library 

2 

S.  S  Chima 

1 

Naval  Avionics  Center 

Indianapolis,  IN  46219-2189 

Library 

2 

Naval  Surface  Weapons  Center 

Dahlgren,  VA  22448 

Library 

2 

T.  C  Henderson 

1 

Naval  Weapons  Center 

China  Lake,  CA  93555 

Library 

2 

Naval  Surface  Warfare  Center 

White  Oak,  MD  20903 

Library 

2 

Pacific  Missile  Test  Center 

PI.  Mugu,  CA  93042 

Library 

2 

U  S.  Naval  Academy 

Annapolis,  MD  21402 

Dirctor  of  Research 

2 

NAVPRO 

Laurei.  MD  20723-6099 

2 

SPAWAR 

Arlington,  VA  20363-5100 

Library 

2 

J  D.  Campbell 

1 

D.  E.  Cowles 

_ 

1 

Requests  lor  coploi  ol  this  report  from  DoD  activities  and  contractor*  ahould  ba  dlrectad  to  DTIC. 
Camaron  Station,  Alaxandrla,  Virginia  22314  ualng  OTIC  Form  1  and,  l(  necessary,  OTIC  Form  55. 


'tnttlal  distribution  of  this  document  within  the  Applied  Physics  Laboratory  tie*  been  made  In  accordance  with  a  Hat  on  Ilia  In  the  API  Technical  Publication*  Q  oup 


INITIAL  DISTRIBUTION  EXTERNAL  TO  THE  APPLIED  PHYSICS  LABORATORY* 


ORGANIZATION 

LOCATION 

ATTENTION 

No.  ol 
Copies 

DEPARTMENT  OF  THE  NAVY  (continued) 

Naval  Weapon  Support  Center 

Crane,  IN  47522 

Library 

ML  Dean 

D.  A.  Winkler 

2 

1 

1 

Naval  Ocean  Systems  Center 

San  Diego,  CA  92152 

Library 

H.  F.  Wong 

R.  B.  Lowry 

1 

1 

1 

DEPARTMENT  OF  THE  ARMY 

Ballistic  Missile  Defense  Advanced 

Technology  Center 

Huntsville,  AL  35807 

Library 

1 

Missile  Command 

Huntsville,  AL  35898 

Library 

i 

DEPARTMENT  OF  THE  AIR  FORCE 

Atr  Force  Avionics  Laboratory 

Aeronautical  Systems  Division 

Wright-Patterson,  AFB,  OH  45433-6503 

Library 

2 

CONTRACTORS 

Advanced  Technology,  Inc. 

E-Systems/EC! 

12005  Sunrise  alley  Dr. 

Reston,  VA  22091 

1501  72nd  St.  N. 

St.  Petersburg,  FL  33733-2248 

R.  F.  Gregorits 

K.  M.  Allison 

M.  A.  Larson 

T.  R.  Marchant 

M.  A.  Mayor 

R  A  Saraydar 

J.  F.  Slosic 

W  L.  Dillon 

1 

1 

1 

1 

1 

1 

1 

1 

RCA  GE 

Borton  Landing  Road 

Moorestown,  NJ  08057 

P.  D.  Edgar 

D.  J.  Herman 

S.  B.  Jo.ce 

H.  C.  Kelly 

1 

1 

1 

1 

Requests  lor  coplea  ol  thla  report  Irom  DoD  activities  end  contractors  should  be  directed  to  OTIC. 

Cameron  Station,  Alexandria,  Virginia  22314  using  DTIC  Form  1  snd,  II  necssssry,  OTIC  Form  55. 

*  Initial  distribution  el  title  document  within  ttte  Applied  Physics  Laboratory  has  been  made  In  accordance  with  a  llal  on  Ilia  In  the  Aft.  Technical  Publications  Group. 


